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Introduction 



The Standard Model of particle physics (SM) is a renormalisable quantum field theory, based 
on the SU{3) x SU{2) x U{1) gauge group. This model is very successful, since it can explain 
almost all observations. However, a few points arc unexplained. For example, among many other 
questions, we do not understand why there arc three (and only three) generations of fcrmions, 
why the fermions exhibit such a mass pattern, why there is such a flavour structure among quarks 
or leptons. We do not know why the electrical charge is quantified, what is the origin of CP 
violation and the anomally cancellation in the SM looks rather miraculous. If the Higgs mech- 
anism indeed provides masses to the fermions, the stability of the Higgs mass under radiative 
corrections is a problem. There are also a few observations that possibly provide evidence for 
physics beyond the SM. Among the latter, is the origin of neutrino masses. About 3/4 of our 
Universe is constituted of Dark Energ}^ the exact nature of it we do not grasp. Only 1/4 of 
our Universe is made of matter, but only 15% of this matter is known, the remaining part, the 
Dark Matter, being yet unknown. Moreover, wc do not know why is there matter at all, if the 
Universe was initially matter-antimatter symmetric. 

In this thesis, we focus on two of these observations, namely the matter-antimatter asymmetry 
of the Universe and the fact that neutrinos are massive. 

The dominance of matter over antimatter is an obvious and happy constatation that we do every 
day. Observations tell us however that matter dominates, but the domination is really small. 
This tiny number, the baryon asymmetry of the Universe, is unexplained in the SM. In chapter 1 
we will discuss the most important observations, and present some of the possible explanations. 
The other failure of the SM we are interested in stems from the observation of a deficit in solar 
and atmospheric neutrino fluxes. As wc will explain in chapter 2, these deficits imply that neu- 
trinos are massive particles, and this really constitutes a first experimental evidence of physics 
beyond the SM. 

Among the different mechanisms advocated to explain neutrino masses, we consider in more 
detail the class of seesaw models, where some heavy fields which couple to fight neutrinos are 
added to the SM particle content. 

The seesaw model constitutes so far the preferred solution to the non-zero neutrino masses, 
mainly for two reasons. First, these heavy fields naturally emerge in larger gauge groups which 
attempt to unify the SM interactions. The second reason is that these fields provide an expla- 
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nation for the observed mattcr-antimattcr as^'mmctr)', through the leptogenesis mechanism. 
Chapters 3, 4, 5 and 6 arc devoted to the study of leptogenesis, in which the heavy fields in- 
troduced can generate during their decay an excess of leptons over anti-leptons. The lepton 
asymmetry is then partly reprocessed into a baryon asymmetry. Chapter 3 can be viewed as an 
introduction to the work I have achieved during my thesis, which is discussed in chapters 4, 5 
and 6. 

In chapter 3, we discuss in detail leptogenesis, in particular the constraints under which does this 
mechanism accounts for the observed baryon asymmetry. Chapter 3 relies on the old wisdom 
based on the one-flavour approximation. This picture, even if yielding good estimates of the 
lepton asymmetry, is nevertheless incorrect. 

Indeed, as explained in chapter 4, since the interactions involving charged leptons are usually 
in- equilibrium for the temperatures relevant in leptogenesis, lepton flavours have to be consid- 
ered. We show in this chapter under which conditions the flavours must be included, and we 
show what are the good objects to study in this framework. 

In chapter 5, we study the effects of lepton flavours in leptogenesis in the type I seesaw, where 
fermion singlets are added to the SM particle content. We show in this chapter how flavours af- 
fect leptogenesis. This chapter constitutes the main results obtained during my thesis, which are 
important results in all leptogenesis models, since the inclusion of flavours profoundly modifies 
the constraints derived in the one-flavour approximation. 

Finally, in chapter 6 we investigate the effects of flavours in the type II seesaw, in which scalar 
triplets are added to the SM. The framework of this study is a supersymmetric grand-unified 
(SUSY-GUT) model, in which scalar triplets as well as fermion singlets naturally appear. In 
this chapter we study the influence of lepton fiavours, as well as the effect of including second 
lightest fermion singlet in leptogenesis. We further correct the GUT relation between charged 
lepton and down-t3^pe quark masses. These inclusions and corrections are found to be fruitfull, 
since they render this model of SUSY-GUT leptogenesis a viable mechanism to account for the 
observed baryon asymmetry. 



Chapter 1 

Baryon Asymmetry of the Universe 



It is an every-day observation that matter dominates over anti-matter. On Earth, matter is 
present everywhere around us, whereas antimatter is only significantly produced in accelerators. 
The Solar system or our galaxie as well, are clearly matter-antimatter asymmetric. One could 
think that the Universe contains as much matter as anti-matter, each kind living apart in distinct 
regions. However, if it was so, one would expect annihilations X + X ^ I'y to occur on matter- 
antimatter borders, and therefore high-energy photons would be produced with a clear signature. 
As we have not observed such annihilation, we must conclude that our Universe is globally matter- 
antimatter asymmetric. 

More precisely, two distinct cosmological observations, based on two very different physics, give 
us similar probes: in a volume which contains 10^° photons, the number of particles exceeds that 
of antiparticles by 6: 

7,5 = ^^^^^:^6x10-10 . (1.1) 
This quantity rjB is defined as the Baryon Asymmetry of the Universe (BAU) 

1.1 Evidence for a baryon asymmetry 

Let us briefly discuss the experimental evidences of the BAU. 
1.1,1 Big-Bang Nucleosynthesis 

Big-Bang Nucleosynthesis (BBN) |2] is a remarkably trustful window on the early Universe, 
since it is based on well-known Standard Model physics. It predicts the abundance of light 
elements D, ^He, ^He, '^He as a function of a small set of parameters, among which the matter- 
antimatter asymmetry 7]b- 

In the early Universe, at T ~ 20MeV, baryon-antibaryon pairs annihilate, leaving any baryon 
excess frozen-out and constant in a comoving volume. This amount of baryons consists of protons 
and neutrons, which are in thermal equilibrium together with and neutrinos. For T > 1 MeV, 
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these species undergo weak interactions: 



n + Ve ^ p + e 



p + Ve ^ n + e 



which maintain the neutron-to-proton ratio at its equilibrium value: 



n 



eq 



n 



eq 



_ mn-mp _ 1.293 McV 



At high temperatures, protons and neutrons are most likely free particles, as light nuclei are 
disfavoured, due to the high energetic gamma rays of the thermal bath. 

Since the rate of weak interactions, r„^p ~ G'jpT^, decreases faster than the Hubble expansion 
rate, H{T) ~ 1.66 a/^TGjv J"^, the electroweak interactions will freeze as the temperature of 
the Universe cools down. This occurs at Tfr ^ 0.8 MeV and at this temperature, according to 



eq.(1.2|, the neutron-to-proton ratio is about Un/up ~ 1/6. 



During the cooling process, light nuclei become energetically favoured, as the temperature drops 
below their different binding energies. '^He, the more stable of them, is likely to be formed. 
However, direct production, 2p+2n — > '^He+j, is highly inefficient due to the small cross-section 
of the electroweak processes, and also due to the relative low density of species. Therefore '^He 



is more likely produced through a series of chain reactions, which are depicted in fig (1.1 




Figure (1.1): Schematic representation of the nucleosynthesis chain reaction. Picture taken from |3]. 



The deuterium production is a prior to the chain reaction, through the process 

p + n — > D + ^ . 

However, given its small binding energy, Eb — 2.2 MeV, as long as the gamma rays of the thermal 
bath are energetic enough to photo-dissociate D via 'y + D ^ p + n, no nuclei can be efficiently 
produced. This is the so-called "deuterium bottleneck". 
The chain reaction begins when the deuterium to baryon ratio is ~ 1: 

riD fTnucV^"^ ^ . 

— — Vb e^""c ~ 1 . 

rib \fnp J 



1.1.1 Big-Bang Nucleosynthesis 
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This happens for T ~ 0.06 MeV. Once deuterium production starts, the chain reaction rapidly 
occurs, until ^He is synthethised. All neutrons are to be converted in ^He, which mass fraction 
is then roughly found to be: 

Yp = - 1/4. 

Precise computations of this mass fraction show a tiny dependance on the baryon asymmetry 
A fit of Yp, accurate over the value rjB — G x 10^^°, gives 

Yp ~ 0.2485 ± 0.0006 + 0.0016(7yio - 6) , (1.3) 

where r]iQ = IO^^tje, is the normalised baryon-to-photon ratio. 

In fact, the ^He mass fraction is not a good indicator of the baryon-to-photon ratio. Rather, 
it strongly depends on the expansion rate, and so it is a good indicator of the time at which 
nucleosynthesis occured. Indeed, since in a first approximation all neutrons are captured in this 
nuclei, the later this capture happens, the longer time neutrons have to decay, and hence a lower 
mass fraction is obtained. 

The prefered "baryometer" turns out to be deuterium, since its abundance directly results from 
the competition between production and destruction processes. The primordial D to H ratio is 
accurately fitted around its central value by [S]: 

= 2.64 (1± 0.03) f—") , (1.4) 

where yo = W^{D/H). This ratio clearly denotes a stronger dependence on the baryon asym- 
metry. 

Moreover, the post-BBN evolution of the deuterium is easy to track, the latter being only de- 
stroyed in stars. Therefore, observations of D in "old" regions of the Universe yield the value of 
its primordial abundance: 

yD = 2.68t^,il. (1.5) 

Using this observation, the baryon-to-photon ratio is found to be 

r/io = 6.0 ± 0.4 , (1.6) 

which is in excellent agreement with the value predicted by the cosmic microwave background 

observation, as we will see in the next section. 

Comparatively, the ^He primordial mass fraction is inferred to be 

yp = 0.240 ± 0.006 , (1.7) 

which, together with the D abundance, provides the value 



r/io = 5.7±0.4. 



(1.8) 
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Similarly to deuterium and defining the number abundance of various nuclei relative to the 
Hydrogen abundance, 2/3 = \{)^(^He/H) and yu = Li/H), the mass fraction of ^He and 

Li are found to be [3]: 



^ 3.1(1 ± 0.01) 

VLi ^ 



!?io 

8.5 



From these fits, and the corresponding primordial abundances inferred by observation, 7710 can 
be deduced. However, the determinations of the primordial abundance for ^He and Li are more 
complex and model dependent, hence the inferred rjiQ should be regarded with less confidence. 



Fig, (1.2) shows how the analytical fits for the various mass fractions, together with the primordial 



abundances inferred from observation, can result in a determination of the baryon asymmetry. 



Baryon density ^b'^^ 



0.02 0.03 




^Li/Hlp 



Baryon-to-photon ratio rj x 10 



Figure (1.2): Abundance of light elements. Figure taken from |6]. 



In conclusion, BBN tells us that the baryon-to-photon ratio is non-zero, but small: 



■qs = 5.7 ±0.4 X 10^ 



-10 



(1.9) 



1.1,2 Cosmic Microwave Background 

The observation of the CMB from the WMAP satellite gives a very precise determination of 
the baryon-to-photon ratio, which is consistent with the BBN value given above. However, while 
BBN started around T ~ 0.08 MeV, i ~ 3 min and was almost complete at t ~ 20 min, the CMB 
photons come from matter-radiation decoupling, which occured at recombination, T 1 eV or 



1.1.2 Cosmic Microwave Background 
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equivalently t ~ 4 x 10^ years [H]. 

As the temperature drops below T ~ 13.6 eV, photons are no longer energetic enough to photo- 
dissociate Hydrogen, hence the reaction 

e~+p^H + -i (1.10) 

deviates from equilibrium, resulting in the absorption of free electrons by protons and ionized 
atoms. The photons, which were in equilibrium with electrons, via scatterings and annihilations 
processes: 

e — h 7 + 7 , e~ + 7 <-> e~ + 7 , (l-H) 

then decouple and become free streaming particles. The observation of these photons provides 
us a snapshot of the Universe at that time. 

Before decoupling, the photons and the baryons are strongly coupled, and form a fluid which 
is liable to interactions with opposite effects. On the one hand, in a gravitational well, gravity 
pushes this fluid down the well, whereas on the other hand the radiation pressure tends to push 
it out. The fluid thus goes through accoustic oscillations, which last until the baryons and the 
photons decouple. 

After their last scattering, photons freely propagate, and the surface of this last scattering appears 
to us as a sphere of homogeneous temperature Tq ~ 2.273K corresponding to a decoupling tem- 
perature Tdecoupiing — 0.26 eV. Howcver, the temperature distribution is not fully homogeneous, 
and the study of the anisotropics provides an accurate determination of a set of cosmological 
parameters. More specifically, the angular distributions of rescaled anisotropics 

Ar(g,0) ^ T{9,^)-Tmean 
m pan J- m.p.an 



are decomposed on spherical harmonics 

AT{e,ct)) 



J2 E '^imYiie,<P), (1.13) 



i=l m= 



so that the study of the anisotropics is translated into the study of the angular power spectrum 
Ci = (|a^^mP)- The cosmological parameters are extracted after a fit of the power spectrum, over 
a set of priors, which are theoretical assumptions [3] . By modifying the values of the differents 



parameters, the relative height and space between the different peaks vary. For example, fig (1.3) 
show the different power spectra obtained by varying tjb around its central value, increasing or 
lowering the baryon asymmetry up to 50%. 

We clearly see how the first accoustic peak, located at ^ ~ 200, is influenced by the value of rjB- 
This first peak is precisely the one which provides the constraint on the baryon asymmetry. 
Using WMAP 5 years only, the baryon asymmetry is found to be [5] 



6.225 ± 0.17 X 10"^°. (1.14) 
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Figure (1.3): Temperature angular power spectrum. Figure taken from j3] 



The combination of different observations leads [S] to the value 



ryB = 6.11^°|6 xlO-i° 



(1.15) 



at 95% confidence level. This is the value we will use throughout this thesis. 



1.2 Generating the Baryon number of the Universe 



The observed baryon asymmetry calls for explanation. One could argue that this asymmetry 
was an initial condition of the Universe, or that is has been produced in the very early stages of 
its evolution. Assuming that initially the Universe was asymmetric goes against a naturalness 
principle, and is not satisfactory. Furthermore, an inflationary period is believed to have occurred 
in the Standard Model of Cosmology, (see, e.g. ^0]), during which the scale factor ait) of the 
Universe is increased by a factor ~ 10'^'^. Consequently, any asymmetry existing prior to inflation 
would have been extremely diluted, by a factor ~ 10^^°. Therefore, one has to suppose that the 
Universe was initially matter-symmetric, and find a way to produce 7?_b 7^ after the inflation 
period. We know that the annihilation of baryons and antibaryons is not fully efficient, and so 
one can wonder whether such an asymmetry could be produced at the freezing of annihilations. 
Following Appendix B we have that for 1 MeV < T < 1 GeV, baryons and antibaryons are in 
thermal equilibrium, with 



n 



eg 



n 



eq 



m 



rrir 



3/2 



-rrip/T 



(1.16) 



Baryons and anti-baryons annihilate and their density decreases until the reactions freeze-out. 
Since the annihilation rate is Ta ^ ""^"^ '-^'"''^ ,.ruu „ ■i-t^^v.^^iu, „,,„v.„„^,. v.^-*-^ /^i„.i\ 1 /^2 



~ 71^'^ {a\v\), with a thermally averages rate {a\v\) 



^Cf. also chapter 5 and 6 of [T]. 
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the decoupHng occurs when T annihilation — H{T), at a temperature Tq ~ 20MeV. At this 
temperature, the baryon-to-photon ratio is 

eq „e9 

!^ ^ ~ 10-18 « 10-10 , (1.17) 

which is far below the observed matter-antimatter number density: the Standard Model of 
Cosmology cannot explain the baryon asymmetry. Therefore, a model is needed, in which this 
number is generated after inflation and before the freeze-out of baryon-antibaryon annihilations. 

Sakharov's conditions 

It has been shown that three conditions, known as Sakharov's conditions [H], are necessary to 
produce a baryon asymmetry in the Universe, starting from a symmetric initial state: 

• baryon number must be violated, 

• C (Charge conjugation) and CP (charge x parity) must not be conserved, 

• the system needs to undergo an out-of-equilibrium period. 

The first condition is somewhat obvious, if starting from B = one aims at generating B ^ 0. 
The second condition requires both C and CP to be violated. Indeed, for a given process X ^ i 
in which the particle i bears a non-zero baryon number Nb^i and X is any particle, if C were 
conserved, one would have M{X i) = M-{X i) and therefore the global baryon number 
variation IS.B would be approximately given by the difi'erence between baryon number creation 
A6 and antibaryon number creation A6: 

= A6 - A6 oc Nb,^ ^ \M{X ^ + N^-i ^ {M^X ^ i)|2, (1.18) 

X X 

As Nh^i = —Nfji, one trivially observes that no baryon number can be generated. 

If CP is conserved, then M{X ^ i) = M-{i ^ X) = M{X ^ i), the last equality coming 

from CPT invariance. Again, it turns out that no excess of baryons over antibaryons can be 

generated. 

Finally, a departure from thermal equilibrium is needed. Indeed, looking at the Boltzmann 
equation for the number densities, neglecting the expansion term, one has 

^ = Y.J d'^xdHi {nx{l ± n,) r{X ^ i) - (1 ± nx) r(i ^ b)) , (1.19) 

where (1 it n) are Pauli blocking or stimulating factors, for fermions (-I-) or bosons (— ). If the 
difi'erent species i and X were in thermal equilibrium all along their evolution, we would have 

nx{l±ni) = n,{l±nx): (1.20) 

which, combined with the unitary condition: 

Y;^\M{X^a)\ = Y,\M{a^X)\ (1.21) 

X X 
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would ensure that drii/dt = 0. 

Another way to see how these three conditions arise is to consider how the baryon number 
behaves under C, P and T transformations. Being odd under C, we have 



and so C conservation implies (B) = 0. 

Similarly, B is odd under CP so that the same argument holds. 

Finally, in thermal equilibrium the system is described by a density matrix p{t) = e~f^^ , so that 
one has 



Given that initially -B(O) = 0, no baryon asymmetry can be generated if the system is in thermal 
equilibrium. 

These three conditions are the basic requirements that a model should fulfil in order to produce 
a non- vanishing baryon asymmetry. 

1.3 Models of Baryogenesis 

Sakharov's conditions can be satisfied in many different ways. 

The first attempts to explain the BAU were based on Grand-Unified Theories (GUT) [T^. Indeed, 
C and CP violation are easily fulfilled, at least qualitatively if the Yukawa couplings are complex 
and the departure from thermal equilibrium is found to be naturally satisfied during the freeze- 
out of a heavy particle. Moreover, B violation, which was not known to occur in the SM, is a 
common by-product of any GUT theory. 

1.3.1 GUT baryogenesis 

The early realisations of GUT baryogenesis were based on the SU{5) gauge group, in which a su- 
perheavy field, whose mass is related to the various breaking scales (usually 10^^ GeV < Mx < 
Mgut — 10^® GeV), decays when the temperature drops below its mass. This implies a Boltz- 
mann suppression of inverse-reactions, and thus out-of-equilibrium decays. Under the condition 
that the X couplings to its decay-products are complex, these decays can produce a non-zero 
Baryon number, which was supposed to remain constant during the subsequent cooling of the 
Universe. 

However, as already stated by Dimopoulos and Suskind in [12], "a quantum mechanical source 
of baryon-number violation (...) is possible (...) to seriously alter the results of this [paper]". 
Indeed, in 1976, 't Hooft [T3| showed that non-perturbative effects, which he called instantons, 
and which emerge from the SM gauge structure, are a possible source of B violation. This 
violation occurs due to the tunneling over an energy-barrier of the instantons. However, this 
instant on- effect was calculated in zero-temperature field theory, and the probability for such a 




(1.22) 



{B{t)) = tr{e-^^B{t)) = tr{e-^^ e-'^'B{{)) e'^') = (5(0)) . 



(1.23) 



1.3.2 Electroweak baryogenesis 
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B violation to occur was shown to be exponentially suppressed. 

Taking temperature into account, Kuzmin, Rubakov and Shapovnikov showed in 1985 [Qj that 
while the tunnelling through the energy barrier was exponentially suppressed, a non-static gauge 
field configuration, lying on the top of the energy barrier, allows B violating processes to quanti- 
tatively occur. More precisely, they showed that this field configurations, called sphalerons, are 
in equilibrium for 10^ GeV < T < 10^^ GeV. These processes violate B + L, the excess of which 
they rapidly wash-out, whilst conserving the orthogonal B — L. We refer the reader to appendix 
A, where sphalerons are discussed in more detail. 

The GUT early attempts were B — L conserving, and given that any B + L number created before 
sphalerons come into equilibrium is completely erased, such GUT-baryogenesis were ruled-out. 
It is customary now to effectively add another condition to Sakharov's three, which could be 
formulated as "Prevent B number creation from being erased by sphalerons". There are many 
ways to do so. For example, GUT models in which h B — L number is created are not affected 
by sphalerons erasure, as we will see in chapter 6. 

1,3.2 Electroweak baryogenesis 

Another class of models consists in the electroweak-baryogenesis mechanism [T^. While GUT- 
baryogenesis was based on the out-of-equilibrium decays of a superheavy field, electroweak- 
baryogenesis relies on the possibility of having a strong departure from thermal equilibrium 
during the electroweak (EW) phase transition, when the Higgs field(s) acquires a non-zero vaccum 
expectation value (vev). The mechanism can be roughly described as follows. Assuming the 
electroweak phase-transition to be of first order, then, during the transition, two degenerate vacua 
will coexist, one corresponding to the broken phase, with (0) / 0, while the other corresponds 
to the symmetric phase {(p) = 0. As the temperature drops down, the regions of broken phase 
grow until they fill all the space, a phenomena called "bubble nucleation". 

In the SM, the only source of CP violation resides in the 6km phase of the C/sTiWj^ mixing 
matrix for quarks, i.e. originates from the non-degeneracy of up and down quark type masses. 
As fermion mass generation is closely related to the Higgs mechanism, it is natural to think 
that CP will be violated during the EW phase-transition, when going from the unbroken phase 
T > Tew into a bubble of true vacuum T < Tew- 

Since sphalerons are in equilibrium for T > Tew, but highly suppressed for lower temperatures, 
the creation of a net B number in the unbroken phase, close to the bubble frontier, can result in 
a non-zero B number in the broken phase if the phase transition is fast enough compared to the 
characteristic time of the sphaleron erasure. 

Actually, it has been shown that electroweak-baryogenesis does not work in the SM. 
The first reason lies in not having enough sources of CP violation to account for the observed 
baryon asymmetry [17]. The second reason comes from the fact that the strong first-order phase- 
transition needed for this mechanism to occur requires the Higgs mass to be small, uih ^ 40 GeV 



^Cabibbo-Kobayashi-Maskawa. 
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in the SM. Given the LEP result, mfl^ > 114 GeV [IH], the electroweak-scale SM baryogenesis 
is clearly ruled out. 

Supersymmetric extensions of the SM could solve this problem [11], by increasing the number of 
CP violating sources, and providing a sufficiently first- order phase-transition. For the minimal 
supersymmetric (MSSM) extensions to work, some fine-tuning is required [20] • Electroweak- 
baryogenesis in singlet extensions of the MSSM is found to be viable 
We do not consider this class of mechanism in this thesis. 

1.3.3 Affleck-Dine baryogenesis 

This class of models relies on fiat directions of the supersymmetric scalar potentials [52]. During 
infiation, fermion condensates can form along fiat directions of the superpotential and can develop 
large vevs. After infiation these condensates coherently oscillate around their minima, and baryon 
as well as lepton number can be stored in these oscillations. 

Albeit being an interesting possibility, we will also not study Affleck-Dine baryogenesis here. 

1.3.4 Baryogenesis via leptogenesis 

There is another class of models for baryogenesis, which has received an increasing amount of 
attention since the original paper of Fukugita and Yanagida [2B]: the leptogenesis mechanism. 
SU{2) sphalerons violate B + L via the effective interaction term (cf. appendix A) 

= llqiqlqiii. (1.24) 

i 

On the other hand, the B — L direction is conserved by these sphalerons. Hence the proposal 
of leptogenesis [21], [2S]: a lepton number could be generated from the out- of- equilibrium decays 
of a heavy field which couples to leptons. Then this lepton number is converted, at least partly, 
into a baryon number. 

This mechanism has the interesting feature of relating the observed baryon asymmetry of the 
Universe to the yet unexplained fact that the neutrinos are massive particles. 
Indeed, in seesaw models, as we will see in the following chapter, the smallness of the light 
neutrino masses is usually explained through the addition of a heavy field, whose mass scale is 
typically of order of 10^° GeV. This field, which is assumed to be a Majorana spinor, can decay 
into leptons and anti-leptons, and can furthermore distinguish between them, granted that their 
couplings are complex. 

If the seesaw indeed provides mass to the light neutrinos, leptogenesis occurs, at least qualita- 
tively. 

We will see in the next chapters that under a few assumptions this mechanism can naturally 
accommodate the observed baryon asymmetry. 



Chapter 2 

Massive neutrinos 



Neutrinos have a very long history. Their existance was first postulated in 1930 by Pauli [22] 
to rescue the principle of energy-momentum conservation in P decays. Thus, neutrinos are 
even older than the Standard Model. In 1962, muon neutrinos were observed ^7}. Finaly, tau 
neutrinos were observed in 2000 by the experiment DONUT at Fermilab and so the tau 
neutrino became the last observed particle of the SM. Despite this long history, very little is 
known about neutrinos when compared to other species: neutrinos are elusive particles. 
Only left-handed neutrinos have been observed. This is because if a right-handed would 
exist, it would be a singlet under the SM gauge group. Furthermore, given the experimental 
possibilities in the 70's, when the SM was being built, (almost) no evidence of neutrino masses 
were observed. 

Therefore in the SM only the left-handed neutrino component was included. However, this 
construction faced, from 60's onward, two major problems: the solar and atmospheric neutrino 
anomalies. 

2.1 Evidence for neutrino oscillations 

The measurements of neutrino fluxes coming from different sources revealed a discrepency be- 
tween the expected and the observed signals. This anomaly was observed both in the solar 
sector [21] and in the flux of muon neutrinos coming from Earth's atmosphere [SU], and consti- 
tutes one of the first evidences of physics beyond the Standard Model. Here we just briefly discuss 
this question. We refer the reader to the review of Strumia and Vissani fST] for a quasi-exhaustive 
presentation of physics where neutrinos are involved. 

Solar neutrinos 

Historically, this was the first neutrino anomaly. Electron neutrinos are produced in the core of 
the Sun, according to the fusion reaction 

4p + 2e- ^^He + lue, (2.1) 
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and given the small neutrino cross section, they escape from the core. A flux of Ve is then 
measured on Earth. The first evidence came from the Homestake experiment j^, which uses 
Chlorine through the reaction in a water tank 

Ve + ^'^Cl ^^Ar + e" . (2.2) 

Counting the number of ^^Ar, it was possible to determine the flux of solar v^. This flux was 
measured to be only 1/3 of the expected one, hence the solar neutrino anomaly. At that time, 
it was not clear whether new physics was needed, or if solar models had to be revised. Later 
on, the SNO experiment [2S] was built being sensitive to all neutrino flavours, and also able 
to discriminate from f^^,-. The total neutrino flux was measured to be in agreement with 
solar model predictions, while the deficit in was confirmed: whereas only are expected, 
all neutrinos flavours were observed. Being confident on the fact that only the electron flavour 
is created in the Sun, the disappearance of and the appearance of muon and tau flavours 
remained to be explained. 

Atmospheric neutrinos 

Another evidence comes from the detection of neutrinos produced in the upper atmosphere. 
High energetic cosmic rays enter the atmosphere, and interact with nuclei, thus producing pions. 
These pions decay into muons, according to 

vr^ ^ /i~ + . 
The high energetic muons subsequently decay 

^ e~ + + z7e , 
^+ ^ e+ + + i^e , 

leading to about a rate for which is roughly twice that of Vg. Atmospheric neutrino experiments 
measure the value 

which compares the observed ratio to the expected one. While i? = 1 is expected, R ~ 0.65 
is observed [30] • This points towards an appearance of Ve and/or a disappearance of v^. The 
SuperKamiokande experiment [3D] settled this question: the deficit from athmospheric neutrinos 
results from disappearance. 

Interpretation 

Different solutions were advocated to explain these flavour conversions, such as neutrino oscilla- 
tions [S2|, neutrino decays (cf ref. 56 of [3T]), or decoherence of propagating u (ref. 44 of [HT]): 
the SK experiment showed that no-oscillation solutions were ruled-out at 4c7 [5H] . 
The oscillatory solution has as immediate consequence that neutrinos are massive particles. 



2.2 Neutrino mixing in the 3 flavour scheme. 
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2.2 Neutrino mixing in the 3 flavour scheme. 

A neutrino involved in EW interaction is not a mass eigenstate. Labelling the mass basis eigen- 
states by Ui, i = 1,2,3 and the flavour basis eigenstate by Va, a = e,iJ,,T, they are related 
through a unitary transformation: 

3 



(2.4) 



1=1 



The matrix U, called the PMN^ [M] 

mixing matrix, is a 3 x 3 unitary matrix which depends 
on 3 real mixing angles and on 6 CP-odd phases. Some of these phases can be rotated away by 
a redefinition of the fields. After this, for a general n x n unitary matrix, if neutrinos are pure 
Dirac spinors n(n — l)/2 mixing angles and n(n — l)/2 — (n — 1) phases remain, while n(n — l)/2 
mixing angles and n(n — 1)/2 phases if neutrinos are Majorana spinors. It is customary to express 
U as a product of three rotation matrices: 



U 




V 



CIS 





sise 

1 




Cl3 



V 



Cl2 Si2 
-Sl2 Ci2 



.diag ( e*'^^ e*'^^ 



,1 











(2.5) 

with Cij = cos {9ij) and Sij = sin (%). Oij are real mixing angles and 6 the "Dirac" CP phase. If 
neutrinos are Majorana fields, there are two supplementary phases (j)i,(j)2, the "Majorana" CP 
phases, which are absent if neutrinos are Dirac fields. 

Given this, the probability for a neutrino that propagates with a mean energy E to oscillate from 
a flavour a to a flavour (3, after having travelled a distance L, is given by 



j 



2E^ 



4 Re Up, U;^) X sin2 



i>j 



4E 



i>j 



2 



m7 



(2.6) 



where Am? • 

The determination of the mixing parameters results from the measurement of neutrino fluxes 
coming from different sources and at different distances from the emission source. Together with 
the natural sources, man-made sources such as reactors [HS] or neutrino beams are comple- 
mentary in determining U. 

In the atmospheric sector, the disappearance of is explained by an oscillation i^^ <^ Vr, whose 
detection allows to constrain ^23 and Am23, therefore called atmospheric mixing angle and (dif- 
ference of squared) mass(es). In the solar sector, the disappearance of Ve is explained by ^ i^fi 



^Pontecorvo, Maki, Nakagawa and Sakata. 
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oscillations, tlie observation of tlie latter constraining the solar sector, 612 and Amf2. Finally 
the last sector, the so-called Chooz mixing angle ^13 and squared mass difference Am^3, inter- 
venes both in solar and atmospheric oscillations and is very constrained by the CHOOZ reactor 
experiment 



We display in fig (2.1) the status of neutrino mixing in the relevant (Am?- — dij) plane for the 



different sectors. It shows the complementarity of natural and man-made sources of neutrinos, 
leading to the determination of the several parameters. 
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Figure (2.1): Values of squared masses differences and mixings angles determine by global fit of neutrino 
oscillation data. Figures taken from |37) . 



The global fit to these oscillation parameters, taken from [37j, is given in the following table: 



Parameter 


Best fit 


3a 


Amfs/lO-^eV 


7.6 


7.1-8.3 


lAmfgl/lO-^eV 


2.4 


2.0-2.8 


sin^ 6*12 


0.32 


0.26-0.40 


sin^ 6*23 


0.50 


0.34-0.67 


sin^ 6^13 


0.007 


< 0.050 



The atmospheric mixing angle is observed to be maximal, while the solar angle is large but 
non-maximal, at a very high level of confidence. The last mixing angle ^13 is still only upper- 
constrained. The precise value of this mixing angle is very important, since it is necessary for it 
to have a non-zero value if CP violation coming from the Dirac CP phase is to be experimentally 
observed. The observation of a non-zero value for ^13 could also allow to determine the neutrino 
mass ordering. As of today the sign of Amfg ~ Am23 remains unknown, leaving two possibilities 
for the neutrino spectrum: either the mass ordering is normal hierarchic, with mi < 1112 < m^., 
or it can be inverse hierarchic with < nii < m2. Considering the oscillation probability, we 



2.3 Non-oscillation experiments 
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see that only the third term 



oc 



i>j 



( 



Am? . L 



4E 



(2.7) 



can lift this degeneracy. Since the Jarlskog invariant [M] 



U*iUajUf3iU^j\ = J = Sl2 C23 S23 cfg Si3 Sin((5) 



(2.8) 



depends on sin(^i3)sin((5), we see that a non-zero value of ^13 is crucial. 

The oscillation experiments provided the first proof that neutrinos are massive. Nevertheless, 
these oscillations neither depend on the neutrino mass scale nor on the nature of neutrinos, that 
is, on having Dirac or Majorana neutrinos. In the last case, two additionnal CP violating phases 
01^2 exist, which do not enter in the oscillation probability. 

Hence, non-oscillation experiments are required in order to complete the picture. 

2.3 Non-oscillation experiments 

Cosmology 

It may not be straightforward how cosmological observations could constrain neutrino parame- 
ters. Nevertheless, robust constraints are derived, when different observations are combined [5?] . 
Were neutrinos massive enough, they would have contribute to distort the photon power spec- 
trum during recombination. Observations of the CMB alone [30] give an upper-bound on 711^: 



The effect of neutrinos is not seen on the CMB, but rather on large scale structure formation. 
Being massive, neutrinos carry a fraction of the total mass density, which has been precisely 
measured by WMAP, for instance. However this mass fraction is small. Furthermore, neutrinos 
have an important free-streaming length £p, since they are only charged under SU (2). Therefore, 
on scales that are smaller than ip, neutrinos are unable to cluster, that is, to participate in the 
formation of large-scale structures, the clustering on these scales being somewhat delayed. As the 
temperature of the Universe lowers, neutrinos become non-relativistic and their free-streaming 
length is reduced. Hence neutrinos begin to cluster, as ordinary matter. Therefore, the heavier 
the neutrinos are, the earlier do they cluster, allowing large-scale structure formation at higher 
temperature (earlier times). 

The constraint on neutrino mass from WMAP 5 year analysis, including many different obser- 
vations, reads pO] 



Erriu < 1.3 eV (95%C.L.) . 



(2.9) 



Smj.<0.61eV (95%C.L.) . 



(2.10) 
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Beta decay 

The beta decay experiment constitutes the best direct test of neutrino masses. The effect of 
neutrino masses is to distort the end-point of the beta-decay energy spectrum of the beta decay. 
Indeed, in the decay 



d —> u + e + Ve , 



(2.11) 



the energy of the electron \s Ef> = Q — Ey, which is maximal for = Q — m^. Here Q 
stands for the energy released in the f3 decay, e.g. for tritium beta decay, ^He + e + T^, 



Q = 7713 f 



18.6 keV. 



Around the end-point, the energy spectrum of the electron is oc \/{Q — Ef,Y — nT-l^ > and so 
iTT-ue / will imply a deviation from the line Q — E^. So far, the best constraint comes from 
MAINZ [?T] and TROITSK [35] experiments: 



m^^ < 2.2eV. 



(2.12) 



The future beta decay experiment, Katrin [33], which is schedules to begin data taking on 2010, 
is expected to reach a sensivity of 0.2 eV to neutrino masses. 



Neutrinoless double-beta decay 

The experiments of neutrinoless double-beta decay [3^ are very important for neutrino physics, 
since the observation of such a process would imply that neutrinos are Majorana particles. As 



can be seen in fig ]2.3[ Oz^/3/3 consists in the reaction 

{A,Z) -^{A,Z + 2)+2e- , 
which violates lepton number by two units. 



(2.13) 




Figure (2.2): Feynman diagram for neutrinoless double-beta decay. 



The decay rate for this process. 



ro,/3/3 oc \M\^ \m. 



2 

ee I ) 



(2.14) 
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depends on niee which is the ee entry of the neutrino mass matrix 
rriee = ^ rui 

i 

= cos^C^is) (mi e^'^i 0032(^12) + m2 e^'*^ sin2(0i2)) + s\n^{ei^) ■ (2.15) 

We see that these processes also depends on the Majorana CP violating phases (f)i^2- However, 
there is an uncertainty of order 0(3) on the nuclear matrix elements, which further depend on 
the nuclei used [31]. 

The status of Oz//?/3 is controversial, since the observation of a positive signal has been claimed by 
members of the Heidelberg-Moscow experiment [IS], corresponding to \mee\ < 0.1 — 0.9 eV, but 
this signal has not been confirmed by, for example NEMO or CUORICINO [S], which obtained 
the bounds \mee\ ^ 0.7 — 1.2eV and \mee\ < 0.2 — 0.9eV, respectively. Hence the importance of 
future Ov PP experiments. 



2.4 Neutrino mass models 



The observation of neutrino oscillations lead us to the conclusion that neutrinos are massive 
particles, which necessarily implies to somehow extend the Standard Model. In fact, these 
observations provided the first evidence of physics beyond the SM. 

It has been noticed by Weinberg [H] that, given the quantum numbers of lepton and Higgs fields, 
the Lagrangian can contain the following non-renormalisable dimension 5 operator 

^''=' = lfap{ira4>*) {4>^hp) , (2.16) 

which violates lepton number by two units. Actually, this operator is the only dimension 5 
operator built out of SM fields, and invariant under the SM gauge group. 

After the Higgs develops a vev, this term gives a Majorana mass to the light neutrinos rrii, ~ 
fapv"^, the smallness of rrii, coming from the smallness of /. In Weinberg's paper, the coefficient 
fafs was expected to be oc < 10~^'^ /mw, hence producing a neutrino mass m,y ~ 10~^ — 

10"^ eV. 

This operator may result from an underlying theory, where light neutrinos couple to some field, 
stabilised at some high scale M, the integration of which would give rise to the Weinberg operator 
eq.(2.16). Thus, the appeafing feature that light neutrino masses could be a window to (very) 
high energy physics. 

There is other possibifities to generate neutrino mass. For example, a natural solution is to 
introduce right-handed neutrinos, which couple to the left ones via an additional Yukawa coupling 
A, the smallness of m,^ being explained by the smallness of A. This coupfing should then be 

< 

Many other possibifities and models can account for neutrino masses. Given the SM field content, 
let us see what are the possible extensions. 

For that purpose we consider the following fields, whose quantum numbers under SU (3) x 
SU{2)l X U{1)y are given below: 
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field 


label 


G321 


lepton doublet 




(1,2,-1/2) 


charged lepton singlet 


dR 


(1,1,-1) 


neutral lepton singlet 


^R 


(1,1,0) 


Higgs doublet 


4>d 


(1,2,-1/2) 


Higgs doublet* 


4>u 


(1,2,1/2) 



In the above notation, we label the SM Higgs field (/>rf, while (pu = (j)d = i'72</'rf, cr2 being the 
Pauli matrix. 

Given the fields considered, it is straightforward to determine the possible extensions of the SM 
which would give masses to the neutrinos either at tree level or due to radiative corrections: 



Extension of the Higgs sector 


Product 


G32I 


field 


mass term 




(1,1,-1) 


singlet scalar 77+ 


Majorana (LH) mass 




(1,3,-1) 


triplet scalar A 


Majorana (LH) mass 


BR (gl SR 


(1,1,-2) 


singlet scalar k~^^ 


Majorana (RH) mass 


Extension of the fermionic sector 




(1,1,0) 


singlet i^R 


neutrino Dirac type mass 




(1,3,0) 


triplet S 


Majorana (LH) neutrino mass 



In the above table, RH and LH denote right and left handedness, respectively. In addition to the 
mass terms above, the right-handed neutrinos can also have a Majorana mass term ITJi^ vjijMij, 
appearing as a bare mass term in the SM potential. As stated above, the minimal extension of 
the SM consists in including right-handed neutrinos (RHn) which couple to left-handed ones via 
\j '^L (j^ui^ + h.c. When (p acquires a vev, this term provides a conventional Dirac type mass 
for neutrinos niy [uiu]^ + 171^^). For neutrino masses of m,j < 0.5 eV, the neutrino Yukawa 
coupling should be A < 10~^^ ~ 10~^/ie, where he is the Yukawa coupling of the electron. Such 
a model is perfectly viable; after all he ~ 10~^ ht- This model, called z^MSM, has been studied 
for example in [U], and has several cosmological implications. 

Nevertheless, one could argue (even if this cannot be seen has a proof on itself) that since the 
inclusion of right-handed neutrinos induces a violation of the lepton number, it is quite natural to 
link it to some high energy scale at which B — L is broken. Another argument is that such small 
couplings, even if possible, are non-natural. We adopt this point of view and do not consider 
this scenario. 

In the Zee |1S] or Zee-Babu models, which are based on the inclusion of the fields 77+ and 
neutrino masses emerge from radiative effects at either one or two loops, and the suppression 
of neutrino masses comes from these loop factors. We do not consider here these models, but 
examine the class of seesaw models explained hereafter, in which neutrino masses are generated 
at tree level and which up to now provides the preferred explanation for neutrino masses. 



2.5 Seesaw models 
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2.5 Seesaw models 

Common to all these models, is having a suppression of the neutrino mass which indirectly 
originates from the introduction of a new mass scale, typically far higher than the electroweak 
scale. 

2.5.1 Fermion singlet: type I seesaw 

This mechanism was originally the first of the seesaw class [SU], and originates from grand 
unification theory. Indeed, in 50(10) models, the minimal representation which is complex and 
anomaly free, is a 16 spinorial representation. This representation contains all SM fermions fields 
of one generation plus one extra field, singlet under the SM gauge group, which can be identified 
as a right-handed neutrino. 

In these GUT models, typically the RHn is not a gauge singlet. It couples to some Higgs 
representation, which acquires its vev at a high scale, of order 10^^ — 10^^ GeV. Integrating out 



this heavy field at the electroweak scale gives the effective term eq.(2.16). 

Such a GUT embedding is nevertheless not necessary at all, since the Majorana mass term can 
be put by hand as a bare mass term in the Lagrangian, the right-handed neutrino being singlet 
under G321. In such the neutrino Yukawa sector reads 

= -AL ^nl^ K-\ + h.c. , (2.17) 

— T 

where ip^ = C ij: is the charge conjugated of ip. 

After the SM Higgs acquires its vev, ((/>„) = t; ~ 174 GeV, the first term provides a Dirac mass 
for the neutrinos 

-mf^iTz'iyi, (2.18) 

with 

mD = Xv/V2. (2.19) 

Using the identity 

the above mass term can be written 

1 / I niD 




Assuming M ^ inu, the diagonalisation of the matrix above gives us two eigenstates per gener- 
ation, which are Majorana spinors. At leading order, one state is given by 
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with mass M. For now we will call tliis state the riglit-lianded neutrino. Tlie otlier state is at 
leading order 

mostly composed of the left-handed neutrino, as implied from low energy observations. The mass 
of this light neutrino is given by the seesaw relation [SU] : 

m}^ = v^\ikM-^\jk, (2.20) 

in the basis where M is diagonal. 

The light neutrino mass matrix is diagonalised by the PMNS matrix U introduced earlier: 

Dm = U^m^U, (2.21) 

where Dm = diag(mi, m2, ma). 

In addition to the nine low-energy parameters, enlarging the fermionic sector by three right- 
handed neutrinos provides nine extra parameters: three masses -|- three mixing angles -|- three 
CP violating phases. There are different ways of parametrising these unknown parameters, 
namely the bottom-up approach, where one fixes the low-energy sector to reconstruct the high- 
energy one, or the converse top-down approach in which one fixes the high-energy sector. These 
two parametrisations are useful depending on whether one considers either an effective or a full 
theory, respectively. In this thesis we use the intermediate Casas-Ibarra [JT] parametrisation, 
where one fixes both low and high energy parameters to infer the neutrino Yukawa coupling A. 
This useful parametrisation derives from the fact that since A.A^ is an hermitian matrix, it can 
be diagonalised by an orthogonal matrix i?, satisfying R.BT = 1. The matrix R is parametrised 
in terms of three complex angles. Here we use the following parametrisation 



R 



C3C2 C3S2 S3 

-C1S2 - S1S3C2 C1C2 - S1S3S2 S1C3 
S1S2 - C1S3C2 -S1C2 - C1S3S2 C1C3 



(2.22) 



where, as for the PMNS matrix, c[s)ij = cos(sin)(zjj), the ZijS being complex angles. One can 
therefore reconstruct A, according to 

Ai„ = f\/M.i?.V^.C/t) , (2.23) 

V /la 

where M = diag(Mi, M2, M3) is the heavy neutrino mass matrix. Hence, in this parametrisation, 
among the 18 parameters, 15 are free: in the low energy sector, once the solar and atmospheric 
Am^ and mixing angles are fixed, the neutrino mass scale -min(mj)- and the mass ordering 
-sign(Amf3)- remain free, as well as the mixing angle 6*13 and the three CP phases 6 — (j)i^2- In 
the high energy sector, the R matrix and the 3 right-handed neutrinos masses are unconstrained, 
hence there are 9 free parameters]^ As stated above, among these parameters, min(mj), 6*13, 



Actually there are 3 more degrees of freedom which are the CP parities of the RHns, and one should use 



R' = i?.diag(±l, ±1, ±1), with R given by eq.(2.22 1. Here we consider the case R' = R. 
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sign(Am23) and 5 are expected to be determined in a near future, whereas the 11 remaining 
parameters will probably not. This is clearly a problem if one wants to make some predictions, 
hence the will to reduce this number. There are different ways to do so. 

In the above discussion, we assumed that 3 right-handed neutrinos are added to the SM particle 
content. Actually, the minimal extension of the SM requires only 2 RHns to generate the 2 
oscillation frequencies. In this model the lightest LH neutrino is massless min(mj) = 0, 
hence the 6 yet free low-energy parameters are reduced to 4, since additionnaly the Majorana 
CP violating phase associated to min(mi), cpi, can be rotated away. Furthermore, the high- 
energy sector reduces to 4 unknowns, although there is little experimental or theoretical support 
for this hypothesis. 



2.5.2 Scalar triplets: type II seesaw 

This model consists in the addition of a scalar triplet A charged as (1,3,1) under SU{3) x 
SU{2)l X U{1)y This field is a real scalar triplet under SU{2), which lies in the adjoint 

representation. A can couple to a product of two left-handed fields. Writing 



A= h (2-24) 




and assuming that A takes a non-zero vev along its U{l)em singlet direction, (A°) = vl 0, a 
Majorana mass term for light neutrinos can be generated. Indeed, since Qy(A) = —2Qy{(l), 
we have in the Lagrangian the additional term 

^^-l k A (iaa) C^L, ^""'-^"^ -\ (m^O., ^l^CyLi ■ (2.25) 

In this so-called type II seesaw, the light- neutrino mass is given by rrii, = f vl and its suppression 
originates from the small vev of A*^ vl <^ v. This smallness can be understood from the other 
couplings of the triplet that are allowed in the Lagrangian, in particular its coupling with the 
SM Higgs doublet (pa: 

Ca 3 -/UA 0l (i^2) A 0d - Mi TV(At A) + ... (2.26) 

where the dots represent higher order couplings. When (p and A develop their vev, assuming 
Ma ^ V, one obtains a seesaw like relation between vl and v: 

-L = -^«- (2.27) 

This ensures the smallness of ml^ . This mechanism naturally occurs in GUT models, where 
scalar triplets are contained in the higher dimensional representions of Higgses used to break the 
GUT group down to the SM. In such a GUT framework, it is frequent that both type I and type 
II seesaw mechanisms occur. 
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2.5.3 Fermion triplets: type III seesaw 

Finally, a third type of seesaw mechanism can occur, where fermion triplets S are added to 
the SM particle content [Si]. These fermion triplets couple to leptons via additional couplings 
(X Xsl(f>u Ti + h.c. Since these triplets reside in the adjoint representation of SU{2) l, a Majorana 
mass term for S is allowed oc Ms S T,'^. In these scenarios, masses for the light neutrinos are 
generated according to: 



which is similar to the type I seesaw formula. 
Probing the seesaw ? 

All these mechanisms are built to account for a suppressed light neutrino mass. However, the 
scale at which B — L is violated depends on the model and on its input parameters, naturally 
implying very different low-energy signatures [55]. It would be very interesting to determine 
which scenario is the most likely to occur, since this knowledge has important implications, 
notably in the building of a high energy model consistent with low-energy data. 

Seesaw and leptogenesis 

The seesaw models share the interesting feature that the light neutrino mass is suppressed owing 
to the introduction of a intermediate scale -M, Ma or Ms- between the SM one and the GUT 
one, the seesaw scale, at which the lepton number is violated. This is a welcome by-product 
of seesaw, which is used in leptogenesis models. In chapters 3, 4 and 5, we will investigate 
leptogenesis in type I seesaw, while leptogenesis in type I-I-II seesaw is studied in chapter 6. We 
do not consider leptogenesis in the type III seesaw. In [SB], it has been found that the role played 
by S in leptogenesis is similar to the one played by fermion singlets. 




— As Mg ^ As 



(2.28) 



Chapter 3 



Leptogenesis 

in the single flavour approximation 

In the previous chapter we have seen that the seesaw mechanism, in addition to explaining neu- 
trinos masses and mixings, implies the violation of lepton number. This L violation is a welcome 
by-product of these models, and is the building-block of the Fukugita and Yanagida seminal 
paper [23], where the foundations of leptogenesis were first laid. 

In this chapter, we introduce leptogenesis in the so-called one- flavour approximation [57 } [55 1 . 
which in most cases provides good estimates of the baryon asymmetry. Here we focus our dis- 
cussion on the type I seesaw, the type II being latter addressed in chapter 6. 



3.1 Basics of leptogenesis 

In the type I seesaw model, the added right-handed neutrinos are very heavy particles, and only 
interact with leptons and the Higgs boson, via their Yukawa couplings. Being heavy, the right- 
handed neutrinos decouple early from the thermal bath. When they decouple, their decay will 
create leptons and antileptons, and if CP is violated during these decays, an excess of leptons 
over antileptons can be obtained. Once right-handed neutrinos are completely frozen out, the 
lepton-antilepton asymmetry will evolve without being affected, as all other processes are CP 
conserving. Then the fast B + L violating processes that are in- equilibrium partly convert this 
lepton asymmetry into a baryon asymmetry. This very rough picture highlights the three main 
features of leptogenesis: 

• A lepton-antilepton asymmetry is created due to CP-odd couplings/processes. 

• At the same time, leptons and antileptons are liable to inverse reactions that potentially 
wash-out the asymmetry. 

• Sphalerons partly convert the lepton asymmetry into a baryon asymmetry. 



40 



The single flavour approximation 



We will discuss the different points in detail hereafter, but let us first give a qualitative picture 
of thermal Icptogenesis. 

The right-handed neutrinos are assumed to be Majorana spinors, and so can decay into both 
lepton and antilepton, due to the Yukawa coupling 

CZi\a^ilo,(t> + \l^NiI^(t>\ (3.1) 
(.a being the lepton doublet of flavour a, while (j) is the SM Higgs doublet. 

As we saw in the previous chapter, this coupling A is in all generality a 3 x 3 complex matrix, 
which, if complex, implies that A^^'s interactions distinguish between leptons and antileptons: a 
CP asymmetry can be created in the decays of RHns. This CP asymmetry is defined by the 
difference between the decay rate into leptons and antileptons 

r(iv^£^)-r(iv^g)^ (3.2) 

In the thermal scenario of leptogenesis, a population of right-handed neutrino is first created 
from the thermal bath via inverse decays i(p ^ N , Icf) ^ N and scatterings. From these inverse 
decays and scatterings a lepton asymmetry is produced, defined as: 

r.-^^^, (3.3) 
s 

where Yx = nx/s is the comoving number density, the number density to entropy density ratio. 
During the thermalisation of RHns, Yl is proportionnal to —ecpYj^. 

Then, as the temperature of the Universe decreases, RHns begin to decouple and decay, creating 
an asymmetry oc ecpYm, that potentially cancels out with the former one. However, during 
and after thermalisation of RHns, the lepton asymmetry undergoes washout processes and so is 
partly depleted. An exact cancellation is therefore avoided. The surviving asymmetry is then 
oc ecpYn {1 — l/C^o) — ^cpf] 7^ 0, where Cyjo > 1 reflects the washout of the asymmetry. 
rj is called the efficiency factor; it reflects the competition between production and depletion 
of the asymmetry, as well as the ability of those process to thermalise the N's. In thermal 
leptogenesis, rj ~ 0(0. 1). The overall factor is the equifibrium density of RHn at the beginning 
of the leptogenesis epoch, and is roughly 10~^. A baryon asymmetry is then generated owing to 
sphaleron processes that partly convert the L as^^mmetr)' into a B as3^mmetry, with a conversion 
factor that depends on the model and on the temperature, but is roughly —1/2. 
The amount of baryon asymmetry produced by leptogenesis is then approximately 

YB--xlO-^r)ecp, (3.4) 

a value to be compared with the observed baryon asymmetry (cf appendix B for relating Yx to 
fix)- 



(3.5) 



3.2 CP violation 
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It implies that the CP asymmetry generated in RHn decays is required to be 

ecp > few X 10"^ , (3.6) 

which constrains the parameters of the model. Let us now discuss the different points in more 
detail. 



3.2 CP violation 

We have defined ecp as the difference between the decay rate into leptons and antileptons. This 
CP asymmetry emerges from the interference between the tree level decay rate and the 1-loop 
corrected one. Indeed, at tree level, for the decay of a right handed neutrino, and given that 
T{Ni — > ia4>) oc lAaip and T{Ni £a4>) oc |A*jp, trivially no CP asymmetry is possible. 
If we define the tree level decay amplitude by Atree = a, the CP conjugate process is Atree = a*, 
and their respective one-loop corrections are Ai-ioop = al3J^ and Ai-ioop = a*f3*J^, where is 
related to the loop function, and /? is some dimensionless coupling, then the CP asymmetry is 
given by 



£CP 



J[V]{\a{l + PT)\'-\a*{l + P*J^)\') 
f[D] (|a (1 + |2 + |a*(l + (3*J^) P) ' 

-2Im{/5} Im| / [V]J^] 



(3.7) 



where /[P] stands for phase-space integration. 

We see that a non- vanishing ecp requires both (3 and T to be complex, i.e. the 1-loop correc- 
tions have to develop an absorptive part, that is, on-shell particles must be running in the loop. 
Moreover, the tree level and the loop corrections have to differ by a relative phase which comes 
from p. 



Evaluation of ecp 

In this model, the CP asymmetry arises from the interference between the tree level decay 



diagram and the one-loop corrections, whose Feynman diagrams are depicted in fig (3.1 




Ni 



Ni 



p 




Figure (3.1): One-loop diagrams contributing to the asymmetry from the Ni decay. 
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There are two types of corrections to the tree level diagram (a): one from the vertex correction 
(b) and another one from the self-energy correction (c) |60| - |61) . with 

„ ^vertex , self 

£CP = £CP + ^CP ■ 

Since the decaying right-handed neutrinos are unstable particles, one cannot handle them in 
terms of asymptotic in- and out- states. One should rather deduce their properties from their on- 
shell contribution to lepton-lepton scatterings [12] • The latter being stable, the usual formalism 
applies. Then, from the 5— matrix elements of such scatterings, the CP violating properties of 
RHn decays are inferred. 

Nevertheless, it appears that in the case where the differences between RHn masses are large 
compared to their width differences, that is if \M]\f^ — M]\i^\ » \Ti — Tj\, the naive calculation 
and the careful one give the same results. We first assume such non-degeneracy of RHn masses, 
and evaluate e'"^^*^^ and e^^'^^ in the naive prescription. We then discuss the resonant regime, 
and how it affects the self-energy contribution. 

Evaluation of e'(fp^^ 

For the process Ni £aH, the tree level amplitude reads 

^■Mla = i XliUi^{pe)PRUN,{pN) , (3.8) 

where Pr(l) is the right(left) chirality projector Pr(l) = (1 + (— )75)/2- 
The vertex correction reads 

.^vertex ^ ^ ^^(i)^ A*„ A*^ A,^ (p,) tx^, (p^) , (3.9) 

i=l,2,3 13 

where 

.F = Pr f /[P]-^Pr- ^— Pl^ (3-10) 

The absorptive part of the loop-function comes from the fact that internal leptons {£f^) and 
Higgses are on-shell. Moreover, we can neglect their masses, as m£,mH <C 

In the loop function, k, q and qn are internal momenta, with pN = q + qn and k = qn + Pi- 
After a straightforward calculation, one obtains for the loop function: 




As lepton masses are neglected, the part oc ]/>£ will not contribute and we neglect it in the 
following. In the expression above we introduce the function fi{x): 

fi{x) = V^(l-{l + x)log{l + ^)]{m-{l + logx))-{l + x)L2(-^) , (3.12) 



^This assumptions holds at T = 0. When taking temperature into account, particles acquire effective masses 
due to screening of the thermal plasma. While lepton masses can still be neglected, m_ff(r) ~ OAT, and so at 
high temperatures tuh can in fact be heavier than M^i- 
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where 1/2 (a^) is the dilogarithm function. 

The CP conjugated process Ni ^ 1^ + H* reads at tree level: 



'^■M.la = '^^aiVNAPN)PLVl^ , (3.13) 



whereas the corresponding vertex correction reads: 

-i 1 



iKr"= E EW'^^"^^/3^*/3T7^irr^MS^l (3.14) 



i=l,2,3 /3 

Then the evaluation of e^^p'^^ is straightforward. The numerator simply reads 

num(e-'^*-) = -^-i.4 ( Im {/i (mJ/M^) } Im |a;,A., (AtA) . j ) ^-Mf , (3.15) 

* \i=l,2,3 J / 

where we have detailed the different contributions, and used global four-momentum conservation 
Pn = Pi + Ph- The numerator is oc A^. In this expression we write (A'I'A)^^. = \\p Xpj. 
For the denominator one considers only the leading term: 

den^r'^-) = (At a).. . (3.16) 



Thus we find the well known result concerning CP violation coming from the vertex correc- 
tion [M]: 

Notice that the sum is made over j ^ i, since the contribution from Ni running into the loop is 
the same for both CP conjugated processes. In the expression above, the loop function f{x) is 
the imaginary part of /i(x) : 

f{x) = V^(l-{l + x)log{l + ^)y (3.18) 



Evaluation of e' 



self 
CP 



Let us now turn to the evaluation of the the self-energy contribution to the CP asymmetry, still 
in the case of non-degenerate masses. 



While the tree level contribution is given by eq.(3.8|, the one loop correction reads 

i=l,2,3 

X U.M^K , \. Pr / ^—^-2 W ^PL«7V.(pAr) , (3.19) 

with obviously pj = pN-, and Qh = Q — Pn- The overall factor 2 comes from the fact that both 
components of the SU{2) doublet can run into the loop. Analogous to the vertex correction, the 
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contribution from j = i cancels out with the CP conjugated process. 
Neglecting Higgs and lepton masses, the loop function is 

T= [ ^ i \ (3 201 

and is clearly divergent. Using dimensional regularisation, and going on-shell pj^ = Mf, one 
obtains 

^ = (2 " ^ ^ ^"^'"^ " ^ ') ^ ^^'^ ' ^^-''^ 

where e = 4 — d and 7 is the Euler constant. 

The evaluation of the self-energy correction then follows the same lines than the vertex correction. 
The CP conjugate process Ni £a(t>* presents the same divergence as above, but the CP 
asymmetry resulting in the difference of the two partial decay rates is finite: 



(3.22) 




with 

9{x) = • (3.23) 

1 — X 

Had we included lepton number conserving self-energy diagrams, we would have obtained an 
additional contribution 



self _ 

Svr ^ (At A),, 




with gdx) = X ^/^(7(x), which in the hierarchical limit is suppressed when compared to the 



lepton number violating one eq.(3.22). 

Adding the two contributions, and assuming that right-handed neutrinos are non-degenerate, 
the CP asymmetry in the lepton flavour a produced by the decay of a right-handed neutrino Ni 
is finally given by [HT] : 

1 Im|A*,A„,- (AtA)..| 
Stt ^ (ATA)ii 

A remark on the calculation of Eqp 

We see that while the vertex-correction is well defined, the self-energy correction is divergent for 
both CP conjugated processes, although the difference is finite for non-degenerate RHn. 
One expects the self-energy CP asymmetry to vanish when the right-handed neutrinos are exactly 
degenerate, since then the CP violating phases of the mixing matrix can be rotated away. We 
see that in the naive estimation above, the function g{x) diverges when x ^ 1, a clear proof that 
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this treatment is incorrect. Actually, this comes from the fact that in the degenerate case, the 
widths of the RHns cannot be neglected, as will be shown below. 

One way to deduce the CP violating properties of N decays is to look at Higgs-lepton two-body 
scatterings mediated by on-shell Ni [n2]-[S3]- The unstable nature of right-handed neutrinos is 
taken into account by summing self-energy corrections in the propagator, which near the pole 
has a Breit-Wigner form: 

^(''' - ,^-M?ViM,r. ' <^-^'" 

where Fj is the width of Ni, given by 

(At A).. Mi 

Ti = ^ . (3.27) 

It has been found in [5^ that the self-energy correction reads (in the simplified 2 RHn scheme): 



with Ofc = Tk/Mi. 

Then, in the non-degenerate case, for \Mn. — Mn^] 3> iTj — Fjl, one recovers the function g{x) of 
eq.(3.23). When RHns are exactly degenerate, no divergence appears as g'^°^'{l) = 0, as expected. 
Having evaluated the CP asymmetry produced by Ni decays, let us now consider the mechanism 
that produces a lepton asymmetry. 



3.3 Leptogenesis in the single flavour approximation 

In the single flavour approximation, one assumes that only the processes involving the lightest 
RHn contribute to the production of a lepton asymmetry. In the case of a hierarchical spectrum. 
Ml <C M2,3, any asymmetry produced during A^2(-^3) leptogenesis, at T ~ M2(M3) is supposed 
be completely wiped out during the A^i leptogenesis era, at T ~ Mi. 

Another approximation concerns the decay products: in this simplified picture only the total 
lepton asymmetry is considered, without taking account the flavour content of the leptons pro- 
duced. Hence A^i decays into L in a CP violating way and the CP asymmetry is obtained by 
summing over lepton flavours: 



Im 



Assuming Mi <^ M2,3, and given the asymptotic behaviour of the loop function 

/(x)+5(x)"-^--^, (3.30) 
2 \/x 



and further using the seesaw formula for the light neutrino mass 



mf = Xo^kM^^Xlp , (3.31) 
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the CP asymmetry can be re- expressed as 



3 Ml Im{(A^.m^A)n} 
= Tr 2 ^\t\^ ■ [6.6Z) 



Davidson and Ibarra have shown [b^ that the CP asymmetry can be bounded from above by a 
function that only depends on right-handed and light neutrino masses. This bound reads 



\ei\<e =— , (3.33) 

where m™'*^(m™™) is the largest (smallest) light neutrino mass. This bound was improved in ^5]: 

kil <e'''/3(mr,mi), (3.34) 
where the function (3 is roughly given by: 

l-m, /mi ifm, « 

In the above equation, the smallest light neutrino mass is compared to the rescaled decay rate, 



The Davidson-Ibarra bound of eq.(3.33| has important consequences for leptogenesis, since it fixes 
the mass scale of right-handed neutrinos. Indeed, given that the baryon asymmetry produced 
by leptogenesis should at least match the observed one 

Yb ^ 1.38 X 10^3 ?7ei > Y§^' ~ 8.7 x lO^^^ , (3.37) 

Ml is lower bounded by: 



Ml > 6.5 X 10* GeV - x 



^ / V^b'' ) V"*™'"'' - "i™'"" /?(m™i°, mi) 
~ 6.5 X 10® GeV . (3.38) 

We will evaluate the efficiency factor r\ later; for the moment it is enough to say that ?7 < 1. We 
thus see that the RHn mass scale has to be far beyond the electroweak scale, and this has many 
important implications. 

One of them is that the seesaw mechanism is hardly testable (if even possible), since such a high 
mass scale will (presumably) never be reached in laboratory experiments, and furthermore since 
the low-energy effects of RHns are highly suppressed. 

Another caveat emerges when leptogenesis is embedded in a supersymmetric (SUSY) framework, 
and stems from the tension with gravitino over-production. We will discuss this issue in chapter 
6, when studying leptogenesis in a SUSY-GUT framework. 
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Basic framework of leptogenesis 

After having discussed its different building-blocks, let us see how leptogenesis works in the single 
flavour approximation]^ 

For the sake of illustration, we first consider a very simplified scheme, where we only include the 
leading terms at order 0{X^) that are decays and inverse-decays, and the on-shell contribution of 
AL = 2 scatterings |5B], included for consistency. We then evaluate the efficiency of leptogenesis 
in different washout regimes, and finally list the various processes which should be included in 
leptogenesis. 

The study of the evolution of the number densities nx{t) is found to be more tractable in terms 
of comoving number densities, Yx = nx/s, since both the number density Hx and the entropy 
density s scale as T^. 

The final lepton asymmetry is evaluated by solving the set of Boltzmann equations (BEs) for Yl 
and Y]\f, which drive their evolution: 



= -D{z)K,{Y^{z)-Y^^ 



= eiDiz)KiiYN{z)-Y^%z))-WUz)K,YL{z), (3.39) 

and compare the comoving number densities to their equilibrium value Y^^. In these equations, 
the evolution parameter is z = Mi/T and the decays D{z) and inverse decays Wid{z) reads: 

D(z) = z^^^^^ 



}C2{z) 



1 

where ICi{z) are the modified Bessel function of the second kind. The lepton asymmetry produc 



W^d{z) = l^^D{z) = -z']C,{z), (3.40) 



tion term is ei Ki D{z), where ei is defined in eq.(3.29), and the parameter Ki refiects how fast 
are decays compared to the Hubble expansion rate: 

_T{ml_mi 

- Him) - ' ^^-^^^ 



alternatively defined in terms of the ratio of the effective neutrino mass mi, defined in eq.(3.36) 
over the equilibrium neutrino mass m*. 



m* = ^ /256_^ ^ ^ ^P_3^^ 

Mpi V 45 ' ^ ^ 

where v ~ 174 GeV is the SM Higgs vev, Mpi ~ 1.22 x lO-'^^GeV is the Planck mass, and 
= 106.75 is the SM effective number of degrees of freedom at T Tew- For i^i » 1 or 
alternatively rh ^> m* , decays are in- equilibrium and lepton asymmetries are strongly washed 
out, while for i^i ^ 1 or m ^ m*, the asymmetries are weakly washed out, as decays are out of 



We use the labelling convention of |59j {D{z), Ki, ...), upon which this section is inspired. 
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equilibrium. 

As for the equilibrium densities, even if a Maxwell-Boltzmann distribution is used for the evalu- 
ation of the interaction rates, we correct them with a factor 3C(3)/4 ~ 0.9, ( being the Riemann 
Zeta function, to match the latter with the high energy behaviour of a Fermi-Dirac distribu- 
tion [SB], so that Y^'' reads: 



^Ni^) = o^4„ 1 ' z'}C2{z) 3.9 X 10-^ . (3.43) 

Z TV 41; 

We see how the Sakharov's conditions are fulfilled in this leptogenesis scenario: 



Baryon number violation: in this case, if L was not violated during decays, we see that 
starting from a vanishing Yl no asymmetry could be generated. 

CP violation: ei / is mandatory. 

Departure from thermal equilibrium: it is clearly seen above that, if right-handed neutrinos 
are in- equilibrium all along their evolution, the production term is null. Notice that this 
last point is only satisfied if the contribution of on-shell RHns to AL = 2 scatterings has 
been included. 



Typical solutions of the set of eqs.(3.39) are depicted in fig (3.2) for various values of Ki, while 



the CP asymmetry is arbitrarily set to ei = 10 



K,=o.m 



a:,=ioo 





Figure (3.2): Evolution of comoving number densities, solutions of eqs.(3.39l, as functions of z 

req 



Nh/T. 

The grey dashed line represents Y^', the number density of Ni in thermal equilibrium, whereas the black 
line stands for Y/v, the solution of the Boltzmann equation. In red (light grey) is depicted the lepton 
asymmetry Y^. For all these plots the CP asymmetry is taken equal to ei = 10^®, while the washout 
parameter is Ki = 0.01(100) on the left (right). 



These solutions exhibit a very different behaviour. 



3.4 Evaluation of the lepton asymmetry 
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3.4 Evaluation of the lepton asymmetry 

This evaluation proceeds mainly in two parts, which correspond to the different contributions to 
the final asymmetry. 

We work in the thermal scenario of leptogenesis, hence before decaying, right-handed neutrinos 
have to be produced. This thermalisation occurs due to scatterings/inverse decays of leptons 
present in the thermal bath. During this first stage, an (anti-) asymmetry is created. 
If Ki S> 1, since decays and inverse decays are far in- equilibrium, the thermalisation of RHns 
is fast and the latter reach thermal equilibrium at high temperatures. Conversely, if Ki <^ 1, 
thermal equilibrium occurs late, as can be seen in fig p:2)| In both cases, during this first stage 
a lepton asymmetry is produced, which is oc £i Y/v. 

Once thermal equilibrium is reached, the subsequent evolution of RHns and of the lepton asym- 
metry also depends on Ki. In the strong washout regime, Ki ^ 1, the fast decays maintain 
Y]y close to its equilibrium value, whereas in the opposite weak washout regime Ki <^ 1, these 
decays occur late, when RHns are far out of equilibrium. During this second stage, a lepton 



asymmetry is produced, which may cancel out with the former one. This can be seen in fig, (3.2) 
where the compensation between first and second stage asymmetry production is visible in the 
peak occuring at z ~ 30(1) for Ki = 0.01(100). 

Finally, as the temperature decreases, as RHns become too diluted, together with the fact that 
the different processes are freezing, the lepton asymmetry reaches its final value. 
Let us evaluate this value in the two characteristic washout regimes. 

Strong washout regime: Ki 1 

In the strong washout regime, the main contribution to the production of a lepton asymmetry 
comes from the second stage, after RHns are thermalised, because thermalisation is fast when 
inverse decays are far in- equilibrium. Thus, defining AA^ = — , the fact that the fast 
processes keep close to its equilibrium abundance implies that: 

AAf'(z) = -Ki D{z) AN{z) - Y^^'iz) ~ . (3.44) 



The formal solution of eq.(3.39| for the lepton asymmetry then reads: 

Yl{z) = ei r dxKiD{x)AN{x)e-^^^-^^''^y'^'^\ (3.45) 

= -ei r dxY^^'{x)e-^^^^^"'^y^'^y . (3.46) 
Jo 



We can rewrite the integrand of eq.(3.45) as / e with 

fix) = Ki r dy WM - HY^'Xx)) . 

J X 

This integral is evaluated using a saddle point approximation, and gets its maximum for 
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that is, for 

5 



2 ProductLog 



5 7rV5 



where ProductLog(2:) is the Lambert W{z) function|^ For Ki > 1, we find that the main 
contribution is for z > 1, as x oc \n{Ki). The lepton asymmetry in this strong washout regime 
is then accurately given by 



YL^er[^^] Y^^^z^n) , (3.47) 



and the efficiency factor is therefore 

rjs = _ (3 48) 

Weak washout regime: Ki <^ 1 



The regime of weak washout is depicted in the left panel of fig, (3.2), In this regime, since in- 
verse decays are slow compared to the Hubble expansion rate, the right-handed neutrinos reach 
thermal equilibrium at Zeq ^ 1. Once thermalised, RHns do not stay in thermal equilibrium, 
but rather decay late. 

These decays produce a lepton asymmetry that cancels out the one produced during thermalisa- 
tion. However, since leptons have more time to participate in inverse reactions, this cancellation 
is only partial, preventing Yi = 0. 
Let us evaluate the different contributions. 

Before the A^s reach their thermal equilibrium abundance, Y^^ dominates over Y/v and so 

Y^{z) ^ K,D{z)Y^''{z), 

Yiiz) ~ -eiKiD{z)Y^''{z)-Wid{z)KiYL{z). (3.49) 
The lepton asymmetry is then given by 

Y<{z) ~ -ei Ki r dxD{x) Y^\x) g-^^ '^y^^^^y^ . (3.50) 
Jo 

Defining for convenience uq = 135C(3)/87r^5f*, such that D{x)Y^''{x) = 4noWid{x), the lepton 
asymmetry produced at z < Zeq is 

y<(z) ~ -ei4no (l -e"l^^i) . (3.51) 

For temperatures below Teg, Yjy dominates over Y^^, in such a way that 

Yl^iz) = -K, D{z) Yn{z) ^ Y>{z) = Y'^\zeq) e-^(^'-^-) . (3.52) 



^The Lambert function has an accurate asymptotic expansion for 2: > 3: W(z) — ln(z) — ln(ln(z)) + ... |67) . 



3.4 Evaluation of the lepton asymmetry 
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The equilibrium temperature is defined by YN{zeq) = Yj^ (zeq)- Given that 



Y<izeg) ~ Ki I dx D{x) Y^\x) ^noKi " 







3 a + 



(3.53) 



where a = and as Y^{zeq) = no z1qK2{zeq) one immediately obtains Zeq'- 



■^eq 



-ProductLog 
2 ^ 



-22/3 r 



1/3 



2/3 



ocln(l/i^i). 



The equilibrium is reached when Y^r = Y^ ^ at Zeq ~ 10 in fig (3.2^ 
We are now able to determine the final lepton asymmetry: 



Yl{z) = ei / dxD{x) {Yn{x)-Y^\x)) e-^'^^^'^y^^'^^y\ 

r" dx D{x) Y^''{x)e-^' dyWUy) + T ^Lx Y^{x) 

10 J Zea 



£1 



£1 



4no(l-e 4^M-noaKi 



11 



£1 X l.-iKlY^^Zi 



(3.54) 



As expected, we find that the lepton asymmetry is non-zero due to the washout of the lepton 
asymmetry produced during thermalisation. The efficiency factor in the weak washout regime is 
then 



1.3 



(3.55) 



Global parametrisation 

Finally, we can give a value of rj for all washouts from a simple interpolation [58]: 



(0.8Kf2 + 2.5i^ 



-l.l6^-l 



(3.56) 



Resulting baryon asymmetry 

As the Universe cools down, at T > Afi/100, all processes are frozen and the lepton asymmetry 
remains constant, until being partly converted by sphalerons into a baryon asymmetry. In the 
single flavour approximation, one has [68J : 



Yb 



14 n„ + 9 uh 



(3.57) 



where n^ = 3 is the number of fermion generations, and uh is the number of Higgs doublets in 
the considered model (nu = 1 in the SM). Hence, 



Yb = -28/51 -1-38 x 10"^£ir/ . 



(3.58) 



This has to be compared with the observed value Y^^ ~ 8.7 x 10 ^^: assuming for example 
e\ = 10~^, K\ should lie between ^ 0.15 and ^ 15. 



52 



The single flavour approximation 



3.5 Dependence on the initial conditions 



In the thermal scenario, since right-handed neutrinos have to be produced by inverse decays, it 
is clear that a minimal amount of washout is required. 

Actually, had we assumed that initially RHn were already thermalised, or even that their abun- 
dance was the dominant one, washout processes would not be required at all. These scenarios, 
which we may call equilibrium [52] or dominant [SB], assume that before leptogenesis occurs, 
for T » Tiepto, some mechanism produces RHns. On the one hand these scenarios are model- 
dependent, but on the other hand the constraints derived for the thermal scenario are weakened. 



To illustrate this, we plot in fig. (3.3) the efficiency factor as a function of Ki, for the different 
scenarios, with Y]y{zi 



mj — in the thermal case, Yn{zv 
in the dominant scenario. 



Yjy (zin) in the equilibrium one and 




Figure (3.3): Efficiency factor of leptogenesis assuming different initial conditions. In red (solid line) 
the thermal scenario is plotted, while in long-dashed blue and short-dashed black, the equilibrium and 
dominant scenarios are respectively represented. 



We observe that in the two latter cases, the efficiency is maximal for a vanishing washout, with 
maximum values r]"^"^ ~ 1 and r/*^"™ ~ 10, while in the thermal scenario a maximal efficiency is 
reached for Ki ~ 2, for which rj^^ ~ 0.16. 

We further notice that for Ki > 10, all scenarios yield the same efficiency: in a sense, the 
strong washout regime is more robust since it does not depend on the thermal history of the 
right-handed neutrinos [M] . 

3.6 Other processes 

In the discussion above we only have included decays and inverse decays. But there are many 
other processes to take into account, that strongly modify the previous results. Let us discuss 
these terms. 



3.6 Other processes 
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©(A^) and 0{X^) processes 

Decays and inverse decays are 0{X'^), but no CP violation is possible at this order, and one 
needs to include 0{\'^) as discussed in the beginning of the chapter. At this order, in addition 
to decays and inverse decays, 2 — to — 2 scatterings mediated by exchange should be included. 

• s— channel exchange: AL = processes £a4> ^ ^pi 

• t— and u— channel A^ exchange: AL = processes 

Ct)4> ^ iai/3- 



and AL = 2 processes ia4> ^ 
(fxp ^ £aip and AL = 2 processes 



We represent in fig ](3.4) the diagrams for the lepton number violating processes. 



Figure (3.4): Feynman diagrams of AL = 2 scatterings mediated by Nj,j — 1, 2, 3. 

The s— channel receives contributions from on-shell RHns, with a corresponding rate 7^ oc 
7d(A)^. As the production term of the lepton asymmetry is ei Ki D{x) oc 7/5 (A)'^, the resonant 
s— channel contribution has to be included for consistency at lowest order, as we did before. 
Then, when including AL = 2 scatterings, this should be carefully done in order to avoid 
double counting, and the on-shell part has to be subtracted from the total AL = 2 interaction 
rates [MIES]- 

When off-shell AL = 2 processes are in- equilibrium, their late scatterings can spoil a successful 
leptogenesis by dramatically washing out the lepton asymmetry. The term that potentially 
washes out the lepton asymmetry is given at low temperature by |70| : 

0.186 / Ml \ / m ^ ^ 



W,,^,(,)^ — y-^'jy — J , (3.59) 

where = ml + -|- m^. However, for z > 15 this term is small compared to inverse decays 
Ki Wid, and can safely be neglected. Since at low temperature inverse decays read 

^1 /^,3/2 



KiWid{z:^l)^^^-z-'^'e-\ (3.60) 
this condition somehow constrains Mi and Ki [71] : 

J5I^<0.1XK,. (3.61) 

If the relation above holds, we can neglect AL = 2 scatterings; on the other hand, if this 
constraint is not satisfied, to prevent the lepton asymmetry from being washed out by these 
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processes implies to upper-constraint the light neutrino mass scale, as we will see at the end of 
the chapter. 



AL = 1 processes 0{X^h^t) and 0{X^h^) 

Given the size of the top Yukawa coupling, processes involving htQ^Hu should be included. 
At 0{\^h1), these processes consist in 2 ^ 2 Higgs mediated scatterings, both in s— channel 
^ iQs and t—,u— channels (Q^, ^ Nt and £t ^ N Qs. These processes are represented 



in fig, (3.5) 



At this order in the couplings, three body decays and inverse decays N ^ £Q^t should also 





£a 




















Q3 




t 
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Figure (3.5): Feynman diagrams of neutrino-top quark scatterings. 

be included. Actually, these processes are phase-space suppressed, and are not quantitatively 
relevant |75j, hence we will neglect them. The scattering processes violate lepton number by one 
unit, and act as source term for heavy neutrino production and as damping term for the lepton 
asymmetry. 

The interesting feature of these scatterings is that at high temperature, contrary to decays and 
inverse decays that, scaling suppressed, the scatterings involving top quarks become 

constant at high temperature. Thus, for z <C 1, thermalisation of RHns is significantly faster. 
When evaluating the lepton asymmetry in the regime i^i ^ 1, we noted that Yl 7^ only owing 
to the partial washout due to inverse decays, and this dependence was illustrated in the fact 
that Yl oc K^. Including AL = 1 scatterings strongly modifies the situation, and the lepton 
asymmetry now scales as Yl oc K. 

However, going to order 0{X'^hf), CP is violated in AL = 1 scatterings, which then also act as 
source terms for the lepton asymmetry. One then recovers the dependence of Y]^ oc in the 
weak washout regime, as source and damping terms now similarly act and cancel each other at 
order K. 

It has been shown in [TH [7H] that the CP asymmetry in scatterings is the same as in decays and 
inverse decays. For consistency, 2^3 scatterings should be included, but as for the three body 
decays they are phase-space suppressed and do not quantitatively contribute in leptogenesis. 
We illustrate the effect of top-quark scatterings in fig, (3.6)[ where we compare the lepton asym- 
metry computed with and without the AL = 1 processes (included both in washout and in source 
terms), in the regimes of weak and strong washout. When AL = 1 top scatterings are included. 



3.6 Other processes 
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the Boltzmann equations for the comoving number densities are: 



dz 
dYLjz) 
dz 



-P{z)K,{Yr,{z)-Y^\z)), 
ei P{z) Ki {Ym{z) - Y^\z)) - Ki W{z) Yl{z) , 



where the production and washout terms are defined by 



P{z) 
W{z) 



D{z)+2Ss{z) + ASt{z), 



D(z) 



YZ\z) 



(3.62) 



(3.63) 



where Sg and St represent the Higgs-mediated scatterings in the s— and t— channels. We refer 
the reader to appendix B where these scatterings are discussed in more detail. 




10"^ 10"' 




Figure (3.6): Influence of neutrino-top quark scatterings on the baryon asymmetry. The dashed- red line 
depicts the case in which scatterings are neglected, the green (light grey) line the case where scatterings 
are included only as damping term, while the blue (dark grey) curve represents the case where scatterings 
are added both as source and damping terms. 



We clearly see the influence of scatterings in the plots: when one considers only the contribution 
of scatterings to the washout, they prevent the cancellation of the lepton asymmetry produced 
during thermalisation with the one generated in decays of the RHn, so that the final Y^ is oc Ki 
instead of being oc K\. When the CP asymmetry in scatterings is included, we see that Yl gets 
suppressed as K1. 

Among AL = 1 scatterings, one should take into account scatterings involving Higgs bosons-?7(l) 
or SU{2) gauge bosons coupHngs [SH1[73]. These scatterings behave similarly to the top-quark 
scatterings, accelerating the thermalisation of right-handed neutrinos, and acting as a source and 
damping terms for the lepton asymmetry. We neglect these terms in this thesis, even if they 
are of the same order -or even larger than the top-quark scatterings-, as their inclusion does not 
drastically modify the picture here drawn. 
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3.7 Upper-bound on the light neutrino mass scale 

The requirement that leptogenesis should explain the observed baryon asymmetry of the Universe 
has led to a constraint on the light neutrino mass scale. Indeed, we have in the single flavour 



approximation the (improved) Davidson-Ibarra bound on the CP asymmetry eqs.(3.33 )-(3.34): 



|,.|<^ ^^-""-"r) ,g(„„,,^,,_ (3.64, 

Ibvr f ^ 

which decreases as the light neutrinos become degenerate in mass. Indeed, for mi ~ m2 ~ ma ~ 
m, and assuming Mi ^ M2, M3, one has: 

^1 ^: 2 — ^ P\J^v [6.bb) 

62tt mv^ 

In this quasi-degenerate case, m > matm ^ m*, and as mi < mi < ms, we are in the strong 



washout regime. Thus we can apply eq.(5.28) which tells us that Yb oc 1/K, hence 



Q Ml 

Yb ~ -1.38 X 10"^eir/(mi) oc ^ . (3.66) 

m^ 

In order to keep Yb large enough, increasing the light neutrino mass scale could be compensated 
by increasing RHn mass. On the other hand, one cannot increase Mi harmlessly. 
There is a strong constraint coming from the requirement that off-shell AL = 2 scatterings do not 
wash-out the lepton asymmetry. These scatterings are potentially dangerous at low temperature, 
when other processes relevant in leptogenesis are already frozen. Their contribution is given in 
this case by [TD] 

3x0.186/ Ml \ / m \2 
" (lO-GeV ) (lev) " ^'-''^ 

This term implies a washout of the baryon (lepton) asymmetry, which then can be expressed as 

3x0.186 f Ml N / ^ N 2 

YB = Y<e lioiuGcvjUovJ ^ (358) 



where Y^ is the would-be baryon asymmetry without AL = 2 scatterings of eq.(3.66). Max- 
imising Yb with respect to Mi and to the effective neutrino mass m, and then requiring that the 
baryon asymmetry at least matches the observed value, the authors of [TH EH ES] have derived 
an upper-bound on the neutrino masses mi,i = 1,2,3: 

mi<0.12eV. (3.69) 

This bound on the light neutrino mass is certainly one of the most important phenomenological 
implication of successful leptogenesis in the one flavour approximation. 

However, we will see in chapter 5 that this bound will no longer hold when lepton flavours are 
included. 

Finally, let us remark that this strong bound only applies in the case the heavy neutrinos are 
hierarchical. Indeed, as shown in [HS], the DI bound, which is roughly 



~ ^ X "^^-"^^ , (3.70) 
Mf m,- ' ^ ' 



3.7 Upper-bound on the light neutrino mass scale 
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gets suppressed for two reasons. On the one hand, the second factor comes from the orthogonality 
of the R matrix and explains the strong suppression of the CP asymmetry in the limit of 
degenerate light neutrinos. On the other hand, the first term of eq.(3.70 ) implies a suppression of 
the CP asymmetry when heavy neutrinos are very hierarchical. However, since we are interested 
in the regime of degenerate light neutrinos, it might be more natural to suppose that the heavy 
neutrinos are degenerate too, unless postulating accidental compensation between right-handed 
neutrino masses and neutrino Yukawa couplings. When the RHns are close to the resonance, 
both suppressions do not occur. The first one gets enhanced by a resonance factor, while the 
second suppression does not occur due to the absence of an orthogonality relation. The author 
of [15] have thus shown that in the single flavour approximation, the constraint eq.(3.69| can be 



avoided. Nevertheless, we will show in chapter 5 that assuming hierarchical RHns (even a very 
weak hierarchy), the inclusion of flavours relaxes the constraint eq.(3.69|. 



58 The single flavour approximation 



Chapter 4 



On the importance of lepton flavours 
in leptogenesis 

In the previous chapter we have studied leptogenesis in the one- flavour approximation, where one 
deals with the total lepton asymmetry, which is the sum over lepton flavours L = L^. Right- 
handed neutrinos decay into leptons and antileptons, and in this approximation one considers 
that the lepton asymmetry produced evolves coherently, until being converted into a baryon 
asymmetry through the B + L violating sphalerons. 

However, lepton doublets interact also via the Yukawa term C D —ha (t> ^% which gives mass to 
charged leptons after electroweak symmetry breaking. This interaction, if efficient, will decohere 
the lepton doublet, by projecting it onto a flavour basis, in which case the electron, muon and 
tau flavours can possibly be differentiated [TS]- [75| [7^. 

If flavours are distinguishable, we should consider the asymmetry produced in each lepton flavour 
direction, instead of the collective asymmetry, for which accidental cancellations are always 
possible. 

In this chapter, which is mostly based on the appendix of [7S|, we do not investigate the effects 
of lepton flavours in leptogenesis. This study will be led in chapter 5 and 6, in the framework 
of type I and type II seesaws, respectively. Here, we discuss how to handle flavoured lepton 
asymmetries: how should the Boltzmann equations be modified, what are the good quantities to 
consider, and especially under which conditions are flavours relevant in leptogenesis. 

4.1 A toy model for flavoured Boltzmann equations 

We want to derive the set of Boltzmann equations for lepton flavours. We could infer, given the 
BE for the lepton asymmetry, the structure of these equations. However, in this approach, the 
flavour structure, and more importantly, the relevance of flavour, may appear rather ambiguous. 
We thus derive the BEs, looking at the evolution of the flavoured leptonic number density 
operator 

paB pais rafS ft /a i\ 

/a£ =fi-f/ = - ' (4-1) 
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where aj^(a£^) and {(^ifj) ^ire the particle and antiparticle number creation (annihilation) op- 
erators, respectively. Actually, this number density operator is not truly the leptonic one, but 
rather the leptonic doublet one, since lepton number is stored both in SU{2) doublets and sin- 



glets. Moreover, as seen from eq.(|4.1|), /^^ is a matrix in flavour space [771 175]. 



A complete derivation of the EEs would require to also consider the Lorentz structure of the 
different fields. Here, we just give an heuristic derivation based on simple quantum mechanics, 
and reintroduce by hand at a later stage the Lorentzian structure. Nevertheless, such a treatment 
enables us to determine the flavour structure of the EE, and furthermore to determine explicitely 
which objects should be studied in flavoured leptogenesis. 

The operators we introduce obey the same statistic that the fields they represent, namely anti- 
commutation relations for the fermionic fields t. 



and similarly for the and A^i. The Higgs field obeys bosonic commutation relations: 







. 



(4.2) 



(4.3) 



We first discuss the Ni — £ interaction. We only consider the interaction with the lightest right- 
handed neutrino A^i, droping the indice in A'^i. The interaction Hamiltonian can be written 
as 



H 



N 



(4.4) 



The evolution of the leptonic doublet number density is given by the perturbative expansion of 
the Heisenberg equations of motion (we use here the expectation value instead of the operator): 

Off 



dt 



H 



N 



N ra/3 



Ki\\, [fN{i + m^pp - ff) - (1 - fN)uff) 
(/iv(i + /0)(5ap - /r) - (1 - /iv)^/r) ■ 



(4.5) 



In this expression, we define the RHn and Higgs distribution function as Jm = aJ^aN and 



The first and third terms of this expression represent decays of right-handed neutrinos whereas 
the second and fourth terms stand for the inverse decays. If we neglect Pauli blocking and Eose 
enhancement factors (1 it /), we can rewrite this expression as 



Off 



2/„(AAt) -/,{/,. (AAt)} 



dt " \ / al3 - ^ \ / J af3 

The evolution of the antileptonic number density is directly deduced from the previous one: 



(4.6) 



af3 



dt •' \ /a/3 



fh ( AAt 



(4.7) 



al3 



4.1 A toy model for flavoured Boltzmann equations 
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with (AA''')"'^ = (^AqiAJ^^ . For simplicity, we have taken the Higgs field to be reaj^ We therefore 
obtain the evolution of the leptonic doublet asymmetry 



dt 



Vn f (AAt) - (AAt)t 



Q/3 



fAe, 



(AAt) + (AAt)t\"'^ 



+ 2 /r, 



(AAt)-(AAt)tr''' 



where we define the "equilibrium" leptonic number density 

/•a/3 I paf3 q reqaf3 



(4.8) 



(4.9) 



The first and third terms of eq 4^ are nul at this order of expansion in the couplings. Forgetting 
about it for a while, eq j4.8| can be further simplified if we assume that the Higgs fields are in 
thermal equilibrium = /^'', and use the equilibrium condition for decays = f^'^ fe, that, 
once symmetrised, reads 



peg I fcq 



We thus obtain 

dff 



'Ae 



dt 



2 {In - fr!) ((AAt) - (AAt)t 



A 



af3 



(AAt) + (AAt)t]"^ 



(4.10) 



(4.11) 



The first term represents the production of a lepton asymmetry in the decays of N . Having only 
made the expansion to second-order, no CP asymmetry can emerge, and the production term 
is nul. However, if we introduce in the interaction Hamiltonian an effective coupling mimicking 
the one obtained when calculating ecp, cf. eq.( 3.29[ ), that is 



iJ^'^/Z ^ ^ A*^A*^Ai^ X g{xj)aNa\a\^ + ^ XjcxXj^Kn x g{xj)a^^a4,ai^ 

j,l3 



(4.12) 



we can obtain non-zero CP violation. Assuming this to be the case, and defining the washout 
factor Ki = Ylia l^aiP and the flavoured washout factors = (AA''")"^, we rewrite the evolution 
equation as: 



a/3 



(4.13) 



For the moment, we have only considered the neutrino Yukawa couplings, but we also have to 
include the charged lepton Yukawa couplings, among which we will only consider the interaction 
involving the tau Yukawa. Assuming that the process (j) ^ ^tT]^ is in- equilibrium, it contributes 
to the interaction Hamiltonian as 



holcal 



(4.14) 



^Had we considered the Higgs to be a complex scalar field, we would have obtained an asymmetry stored in 
the Higgs degrees of freedom, which could be expressed in terms of the B/3 — La one, via a conversion matrix, 

Q in EUES]. 
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In the basis where the charged Yukawa couplings are diagonal, hpr = hrSpr, their contribution 
to the lepton doublet asymmetry evolution can be written as: 

- \hr\^ {fAr^fr + fffAiY^ i^ar + Sf3r) , (4.15) 

with /at': = /r<= - ff^ and = {fr'^ + ffc)/2. 

Before moving further on, it is instructive to look at the equation for the equilibrium density 
defined eq.(4.9|. Neglecting the contribution from the interaction Hamiltonian H'^ , we find 

= 2|/..|2 5.Jpr fl' - \hr\'fj {6.r f^^ + ^r^) . (4.16) 

Here we took two flavours into account, which are labelled r and ^ for simplicity. We see that 
neither the /i^ component, nor the rr one are affected by the charged Yukawa interactions, owing 
to the equilibrium condition of the interaction rate, which is similar to eq.(4.10). On the other 
hand, for the non-diagonal entries, we have 

-^ = -\K?rrU7'\ (4.17) 

and similarly for f'p^'^. 

We see that when the charged Yukawas are in- equilibrium, "equilibrium" densities of quantum- 
correlations are exponentially damped [75], and do not affect the flavour- eigenst ate "equilibrium" 
densities. Therefore, we now write /^'^"^ = f^'^Saf^- 



4.2 Flavoured Boltzmann Equations 



Eq.(4.15) above clearly shows the flavour structure that the Boltzmann equations should have. 
Before discussing this, we make a few transformations in order to extrapolate this equation to 
the comoving number density asymmetry, through the following "rough" procedure]^ 

• we factorise the equilibrium distribution function, using for the last two terms the equilib- 
rium conditions f^'^f^ = and fff^ = f^, 

• we "sort of" integrate over the phase-space to work with number densities, 

• we reintroduce the expansion of the Universe through 

^---3i/(T), 
dt dt ^ ' ' 



^Had we kept the spinorial structure, all the steps indicated would have been rigorous. Since we are confident 
on this point, we use a simple matching procedure for our equations. Following the different steps indicated does 
allow us to identify the different terms without ambiguity. 
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and change time to z = Mi/T and number densities to comoving number densities. At 
this point we get, 



dY, 



a/3 



1 



dz 



zsH{T) 



zsH{T) 



n 



"-N 



1 



{Ya£, k} 



a/3 



n 



eq 



V^q 



Y^c 

|/irP ( >Ar<= ^aP + ^Eq ) C*^"^ + V) 



The last step consists in "sort of" redefining the couplings in order to reintroduce the usual inter- 
action rates. In an accurate calculation, the interaction rates naturally appear when integrating 
over the phase-space. 

Further using the fact that Y^c = Y^ we obtain the flavoured Boltzmann equations for the 
comoving number densities: 



dY, 



a/3 



ea/3KiD{z){YN - Y^") - Wid{z)]^{YM, 



z 7t 
sH{Mi) y/ 



^ ( 2 Y/^T-<= Sa/s + Y^^ 



a/3\ / "ar 



^ar + ^0T 



(4.18) 



The last term represents the depletion of the lepton SU{2) doublet asymmetries by the charged 
Yukawa interaction, with [^.[50]: 



Ja_ 
eq 



n 



(4.19) 



When extrapolating from eq.(4.15) to the BE for comoving number densities, the matching with 



the t'^ interaction term may seem arbitrary. However, we can derive the evolution equation for 
the r'^ asymmetry [SU], by looking at the evolution of the operator 



f — i t _ 

J At" — a^car<= — a-cO'T'' ■ 



(4.20) 



Doing so, and using the same matching procedure as before, we find 

dY/^r"^ z 



dz 



sH{Mi 



-2|/ip|^(2yA.c + y^ 



Ail 



z 7. 



sH{Mi) y 



eq (2 ^At<= + ^aL 



(4.21) 



hence justifying the quoted rate eq.(4.19). 



Given the expression above we can now deduce the temperature regime for which lepton flavours 
are relevant in leptogenesis. 

We see that the interaction projecting lepton doublets onto the flavour space is in- equilibrium 
if [7S1ES] 



z-/^>H{Mi 



(4.22) 
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which provides a bound on the decaying right-handed neutrino mass. Evaluating the equilibrium 
condition at the temperature for which the asymmetry will be mostly produced, T ~ Mi, this 
bound reads: 

Ml < 7 X W^^hl GeV , (4.23) 

Therefore: 

• r- Yukawas are in- equilibrium if Mi < 10^^ GeV, 

• muon- Yukawas are in- equilibrium for Mi < 3 x 10^ GeV, 

• electron- Yukawas are in- equilibrium for Mi < 6 x 10^ GeV. 

However, thermal equilibrium does not guarantee that the lepton doublet is to be projected onto 



flavour space [HI]. Indeed, we see that the terms involving N — £ couplings in eq.(4.18| are 
independent of the choice of the flavour basis, which is obviously not the case for the charged 
Yukawa interactions. The condition for the flavour basis to be relevant is actually more stringent 



than the equilibrium condition eq.(4.23|. The processes decohering the lepton doublets have 
to be faster than the processes which recohere lepton doublets at the time the asymmetry is 
produced, that is [5T| 15^ 



Considering only inverse decays, and evaluating them at their maxima, z ~ 2, one finds 

Ml < 5 X 10^^ GeV^ , (4.25) 

where Kaa parametrises the individual washout of the flavour a. For the flavour r, this constraint 
reads: 

M.<i5^. (4.26, 

Under this constraint, flavours are fully relevant in leptogenesis. The situation where charged 
Yukawas are in- equilibrium but not faster than recohering processes is more involved, and we 
postpone its discussion until the end of the chapter. 

4.3 Flavour structure 
4.3.1 One flavour scheme 

This case is encountered if Mi > 10^^ GeV, and in this case none of the charged Yukawas are 
in- equilibrium. It has been shown in [52] that by a rotation of lepton doublet states, one can 
recover the BE of the one flavour case studied in the previous chapter, that is 
dYN{z) 



dz 



-D{z)K,{y^,{z)-Y;,\{z)) , 



= eiD{z)Ki{YMAz)-Y'^\{z))-W,d{z)KiYL , (4.27) 
with £1 = Tr[£\. For Mi > 10-*^^ GeV, the results of the previous chapter hold. 



4.3.2 Two flavour scheme 
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4.3.2 Two flavour scheme 

Let us now consider the case where only the tau Yukawas are in equilibrium, that is 3 x 10^ GeV < 
Ml < 10^2 GeV. 

In this case, the tau flavour is distinguishable, but neither the muon nor the electron flavours 
are. Hence, the lepton produced in A'^i decays will be projected onto the tau direction, as well as 
onto an orthogonal direction to that of the tau. This second direction is composed of the muon 
and electron flavour directions that coherently mix, and which we will label L± = Lg + L^. We 



have a two flavour regime and in this case eq.(4.18) holds, with a = t and /3 =_L: 

dYN{z) 



dz 
dYl,^{z) 
dz 



-D{z)K,{YnAz)-Y'^\{z)) , 



e^LKiD{z){YN{z) - Y;,\z)) - Wu{z)K^^Y^t{^) 

- Wid{z)\ [n^rYliiz) + Ytj{z)KrX 

= e^rK,D{z) {Y^{z) - Y^^z)) - W.,{z) Y^[{z) ^^^1^ 
f Yl,Hz)+Y^}{z) \ , 

dYliiz) _ dY^liz) 

- ^^^^ 

dYlliz) _ dYttiz) 



{^^T)-Dr{2YAAz) + Yl}{z)) . (4.28) 
In these equations, the CP asymmetries e_L_L and e^r are defined by 

L = ^ee ~l~ 1 

and similar relations hold for the washout parameters k. We also define for convenience 

z -fr TxlOi^GeV o TxlQi^GeV 

Notice that, according to eq.( |4.21 ), the Boltzmann equations for the r diagonal entry can be 
written: 

dJJM ^ dY^ dYA^^ 
dz dz dz 

Since the leptonic r doublet i'^'^ interacts with right-handed neutrinos, its evolution equation 
receives contributions from the interaction Hamiltonian , corresponding to an asymmetry 
production and depletion. However, since i'^'^ also interacts with the right-handed singlet t'^, it 
receives a further contribution from H'^. 

This is quite satisfactory, since we do not want to study the asymmetry stored in lepton doublets, 
but the lepton asymmetry 

La = Lia + Leg = Lga — Le^c , (4-31) 



66 



Importance of lepton flavours 



which is stored in lepton doublets as well as in singlets 
The muon- Yukawas being out of equilibrium, the leptonic asymmetry Yg± is not influenced by 
any effect of this type: muon (or electron) isosinglet density is too far from thermal equilibrium 
to disrupt the i± propagation. 

The off-diagonal terms are however sensitive to the population, which acts as a damping term. 
Furthermore, the term oc Dt-Y^^ exponentially drives the quantum- correlations to zero: when 
the interactions involving charged Yukawa couplings of a given flavour are in- equilibrium, the 
quantum correlations between this flavour and the orthogonal one(s) are exponentially damped, 
and the flavours, being distinguishable, evolve separately. Notice that since we did not include 
the free Hamiltonian oc uj-j- a\^ai^ + w± a\^ai^, the number densities are not liable to flavour 
oscillations [75], which have been studied in detail in |52]. Actually these oscillations only affect 
the quantum correlations by inducing a faster damping 



We illustrate this damping in fig, (4.1) were we depict the evolution of the different asymmetries 



according to the set of Boltzmann equations in eqs.(4.2l|4.28 




Figure (4.1): Evolution of comoving number densities. The thick lines represent diagonal terms of the 
flavoured Boltzmann equation, with in blue (upper curve), F^J in green (light grey) and in 
red (lower curve), respectively. The thin purple lines represent the quantum-correlations, i.e. the ofF- 
diagonal terms when the charged Yukawas are in-equilibrium. The dashed line stands for Y^l whereas 
the dot-dashed one for Y^-^. The input parameters have been chosen, so that the components are 
distinguishable, with = /tj^,r = 0.1, Kt,_l = i^t,t — 20 for the washout parameters and ej^.j^ ~ 

ej^^ 0.1 X 10~^, er,_L = £r,r = 20 X 10^^ for the CP asymmetries. 



In these equations, we have chosen the washout parameters and CP asymmetries in order to 
distinguish the different components. The effect of charged Yukawa interactions in damping 
the quantum correlations is clear, hence in the temperature regime where these interactions are 
in-equilibrium the off-diagonal terms can be safely neglected. 

Furthermore, we plotted in this graph the leptonic doublet asymmetry (green-light grey) and 



4.3.3 Three flavour scheme 
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the lepton asymmetry = Y^J — Yat'= (blue- upper curve). We see that the singlet contribution 
to the leptonic asymmetry accounts only for 0(1), and thus, for now, we will neglect it on the 
evolution of the asymmetries. 

4.3.3 Three flavour scheme 

If Ml < 3 X 10^ GeV, then both muon and tau Yukawas are in equilibrium. These interactions 
decohere the lepton doublet by projecting it onto a three flavour basis (Lr, L^^, Le). The muon 
and electron flavours being now distinguishable, muonic and electronic quantum-correlations are 
exponentially damped, and one is left with the set of Boltzmann equations for the lepton number 
densities: 

-J^ = eaaK^D{z) {Ym{z) -Y^\z)) -W.,d{z) K^^Y^^iz) , (4.32) 
with a = e, n, t. 



4.4 B — L conversion and the Baryon asymmetry 

In the previous section we have derived the Boltzmann equations for the lepton asymmetries. 
However, as explained in appendix A, sphalerons conserve B/3 — La asymmetries but not 
the leptonic ones. Therefore, it is preferable to work with 

YAa = ^Yab - Yal^ , (4.33) 
computed from the leptonic doublet asymmetries by 

YAa = ^«/3^A£/3 , (4-34) 
P 

where the sum is made over all flavours that were distinguishable during leptogenesis. Details 
about the derivation of the entries of the conversion matrix Aa/s are given in the appendix B. This 
relation is obtained by expressing the chemical potentials of the B/3 — La asymmetries in terms 
of the leptonic doublet ones, and it depends on the different processes that are in- equilibrium at 
the temperature leptogenesis occurs [H^ [H^. 

Since we do not consider B violating processes, the Boltzmann equations for Yat and for the 
B/3 — La asymmetries are finally given by (dropping the redundant double indices): 

dYN{z) 



dz 

dYAa{z 



-D{z)K,{y^A^)-Y'^1{z)) 



- -Sa Ki D{z) {Yn{z) - Y'^\z)) - W{z) Ka Aafi YAfi{z) . (4.35) 

This set of equations is written here in all generality, and the expressions of the production term 
D{z) and of the washout one W{z) depend on the model, as well as the parameters which govern 
these equations, Ka, Ki and Sa- 
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Finally, B/3 — La asymmetries are conserved, until they are converted by sphalerons. The baryon 
asymmetry is then related to the individual asymmetries by: 

Yb = asph X ^ ^Aa , (4.36) 

a 

the sum being over the flavour that were distinguishable when leptogenesis occurs. The conver- 
sion factor is asph ~ 1/3, and its value also depends on the model, and on whether the sphalerons 
freeze before or after the electroweak-phase transition [H5| 15^. 



4.5 A remark on the full flavour regime 



In section 4.3 when we discussed the flavour structure of the BEs, the criterion used for the 



relevance of a given flavour in leptogenesis was presented eq.(4.23) and translates into having 
the interactions involving charged lepton Yukawa couplings in- equilibrium, Fq, > H(Miy). 



However, as explained around eqs.(4.24|-(4.26 ), thermal equilibrium does not guarantee that the 
lepton doublets will be projected onto the flavour basis. A more accurate constraint comes from 
the requirement that the decohering rates 7q are faster than any recohering interaction, either 
the production term D{z) or the washout one W{z) fST]. Considering only inverse decays as 
depletion reactions, and requiring that at least the interaction involving the charged tau Yukawa 
couplings are in- equilibrium, a bound on Mi can be derived: 

2 



IQi^GeV /h. 



hj- 



(4.37) 



If this constraint is satisfied, then the flavours are fully relevant in leptogenesis. 

Given the equilibrium condition for sphalerons. Mi < 10^^ GeV, which is roughly the same than 

the equilibrium condition for tau- Yukawas, we see that only in the strong washout regime is the 



constraint of eq.(4.37) more stringent than the equilibrium condition eq.(4.23|. 



Let us consider a lepton doublet (asymmetry) produced in decays of A^i at a temperature T ~ Mi, 
in the case where the tau- Yukawas are in- equilibrium. If during the caracteristic time for two 
i — H — t'^ interactions to happen, inverse decays have time to occur, the lepton doublet will not 
be projected onto the flavour basis, and the faster the decays/inverse decays, the more coherently 
the lepton doublet will evolve. 

In such case, the lepton asymmetry should be computed by solving the set of BEs including 
off-diagonal terms. As already stated, ref. |52] shows that by an appropriate redefinition of the 
lepton doublets, corresponding to a rotation in flavour space (allowed since the A^-interactions 
are the ^-labelling ones), one recovers the single flavour picture. 

However, through this redefinition, couplings are also modified, and so the differences between 
including or neglecting quantum correlations are less transparent. In order to study the validity 



of the full flavour regime, we solve the sets of eqs.(4.28| with and without quantum correlations. 



We depict in fig. (4.2) the effect of these correlations by showing a contour plot of the ratio 



of the lepton doublet asymmetry for a given flavour computed with and without these terms. 
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1a£*«/^A£--°"*' allowing Mi to vary around the central value Mi = 10^^ GeV and the washout 
parameter around Ki = 1. For simplicity, we have taken the CP asymmetries and the washout 
parameters of the different flavours to be equal. One expects that if the condition of eq.(4.37| is 
not satisfied, the damping of the quantum correlations will be inefficient and these off-diagonal 
terms will significantly contribute to the lepton asymmetry (in fact, given our choice of param- 
eters, in an equal amount). On the other hand, if eq.(4.37) holds, the quantum correlations will 
be damped and the ratio is expected to approach unity. 




10" 10'° 10" 10 

Ml (GeV) 



Figure (4.2): Contour plot of the ratio between the lepton doublet asymmetries computed with and 
without quantum correlations. The black line stands for Y^/ia / Y^£qq°"* = 1, whereas the blue, red, 
green and purple lines represents deviation (from left to right) of 5%, 10%, 30% and 60%, respectively. 



We see in fig, (4.2) that for Mi < 5 x 10^^ GeV, the inclusion of quantum correlations modifies 
the results up to 30%, while for Mi < 10^^ GeV, the approximate eq.(4.35) gives very accurate 
results, and so in that case quantum correlations can be safely neglected. 

To conclude on the validity of eqs.(4.35| regarding a precise determination of the lepton asym- 
metry, we can say that for Mi < 10^^ GeV, the quantum correlations are effectively damped and 
one can use eqs.(4.35|. For 5 x 10^^ GeV > Mi > 10^^ GeV, the results are vahd only up to 
10% — 30%, while for Mi > 5 x 10^^ GeV, a precise computation does indeed require the inclusion 
of quantum correlations. For heavier Mi, the one fiavour approximation is valid. 
However, we stress that the results here presented have been derived in the Standard Model (plus 
RHns). Had we considered a supersymmetric extension, the discussion would have been dramat- 
ically different. Indeed, in SUSY extensions, the charged Yukawa couplings is now corrected by 
the ratio of Higgs vevs, h'^x {\ + tan^(/3)), where tan(/?) = Vu/vd is the ratio of the SUSY Higgs 
vevs. Hence for iixr?{(3)) > 3(10), the results are vahd up to 10% for Mi < 10^^^^^) GeV, so 
that for Ml < 10^^ GeV, as required for sphalerons to be in- equilibrium, flavours are always fully 
relevant in supersymmetric leptogenesis. 
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In the following, we will neglect quantum correlations. 



Chapter 



5 



On the 



role of lepton flavours 
I leptogenesis 



in type 



In this chapter we are going to study leptogenesis in the type I seesaw mechanism, were three 
right-handed neutrinos are added to the SM model particle content in order to explain neutrino 
masses. The Lagrangian contains the extra terms 



which provide masses to neutrinos from the seesaw mechanism. These masses are not yet mea- 
sured, but rather the differences of squared masses are deduced from the oscillations of neutrinos. 
However, should we believe that leptogenesis is indeed the mechanism responsible for the ob- 
served BAU, an upper bound on the light neutrino mass can be derived [7111521 ES], rn-i < 0.12 eV, 
a bound that may be soon tested. 

On the other hand, the right-handed sector is experimentally unconstrained. This is one of the 
gaps of leptogenesis that may be never filled in. Nevertheless, the Davidson- Ibarra bound [M] 
constrains the right-handed neutrino mass scale to be above Mi > 10^ GeV . Moreover, CP 
violation in the leptonic sector -related to the CP phases of U and R- should be quantitatively 
large enough so that CP asymmetry is larger than ecp > 10"'^. 

All these constraints have been derived in the single flavour picture, where one considers that 
the lepton involved in the decays of the RHn evolves coherently after being produced. Nonethe- 
less, we saw in the previous chapter that for Mi < 10^^^^^ GeV, not only are the interactions 
involving charged lepton Yukawa couplings in- equilibrium, but their rate can also be faster than 
recohering processes, the decays and inverse decays, in which case the flavours are relevant and 
should be included for a precise computation of the baryon asymmetry obtained by leptogenesis. 
This chapter, which is mostly based on [7T] and [50], is organised as follows: firstly, we set up the 
framework of flavoured leptogenesis and we qualitatively discuss how lepton flavours affect the 
key-parameters of leptogenesis, namely the CP asymmetries and the washout factors. Secondly, 
we evaluate the resulting baryon asymmetry produced by leptogenesis, and show how it differs 
from the one-flavour picture provided in chapter 3. Finally, we consider the constraints derived 



C D -Xa^N^ -NiMiN[ + h.C. 



(5.1) 
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for the low and high-energy parameters. 



5.1 Flavoured leptogenesis 

We first recall the various definitions and the framework of Boltzmann equations. We consider a 
hierarchical spectrum for the right-handed neutrinos Mi <C M2,M^, and thus neglect contribu- 
tions from the heavier neutrinos in the lepton asymmetry production. Indeed, since in this case 
the leptons produced by leptogenesis participate via inverse decays in the thermalisation of 
A'^i, one can consider that the lepton asymmetry produced during A''2,3 leptogenesis epoch will be 
washed out. However, this assumption, if valid in the single approximation [SS]; is not strictly 
justified when flavours are taken into account IBB]) as the lepton flavour content may differ for 
the lepton involved in leptogenesis and the one involved in A''i leptogenesis. This point will 
be discussed in the following chapter; presently, we only consider the processes that involve A^i. 
The EEs for right-handed neutrinos N = Ni and = B/3 — La asymmetry comoving number 
densities are [7T]: 



dz 

In these equations, individual washout parameters are given by 



with the effective neutrino mass ifia being given by 

^2 

and Ki = J2a '^a- The other quantities have been defined in chapter 3, section 3.3. We never- 
therless recall the EE for the lepton asymmetry in the single flavour picture: 

= eiK,P{z) {Y^{z)-Y^\z)) -W{z)KiYl{z). (5.5) 
The individual CP asymmetries Sa are given by eq.( |3.25| ): 

= ^ (\u)l ^-^J (/(x,) + ,(x,)), (5.6) 



where 



fix) + g{x) = ^ ( ^ - (1 + x) log (1 + -) ) . (5.7) 

^ \ — X X 



The total CP asymmetry is the sum of the individual ones ei = 



The set of equations (5.2) governs the B — La asymmetry evolutions. As explained in the 



5.1 Flavoured leptogenesis 
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previous chapter, these B — La asymmetries are deduced from the leptonic doublet asymmetries, 
by relating the chemical potentials of the B — La asymmetry to the ones of £a according to the 
different processes that are in- equilibrium at the temperature at which leptogenesis occurs. 
We saw that if Mi < 10^ GeV, interactions involving muon and tau charged lepton Yukawa 
couplings are in- equilibrium, and are potentially faster than decays and inverse decays, thus 
projecting the lepton asymmetry onto a three flavour basis {Le, L^, Lr). Given that for T < 
10^ GeV, interactions involving charm, bottom and top quark Yukawa couplings are also in- 
equilibrium, the conversion matrix A reads |71| : 

/ -151/179 20/179 20/179 \ 
A= 25/358 -344/537 14/537 . (5.8) 
V 25/358 14/537 -344/537 / 

One then has to solve the coupled BEs for those three distinguishable flavours. The BEs are 
governed by individual CP asymmetries Sa, a = e, fi,T and individual washouts a = e, ^, r. 
If the decaying RHn is heavier, with 10^ GeV < Mi < 10^^ GeV, then only tau Yukawas are 
in- equilibrium (together with b,c,t quarks), so that the lepton asymmetry is projected onto the 
tau direction {Lr,£T, Kr) and onto an orthogonal direction in flavour space, which coherently 
mixes electron and muon flavours (Lq, Eq, Kq) = (Lg -|- L^, £e + e^, + i^fi)- 
In this two flavours scheme the ^-matrix reads: 

A = ( -417/589 120/589 \ 
y 30/589 -390/589 J ' 

For Ml > 10^^ GeV, none of the charged lepton Yukawas are in- equilibrium, thus the one flavour 
approximation holds and one solves the BE for {Yb-l, ^i, Ki). 
Finally, the baryon asymmetry produced by sphalerons conversion is 



^^ = ^E^Aa, (5.10) 

a 

where the factor 12/37 is obtained in the SM, and the sum is made over the flavours which were 
distinguishable when the asymmetries were produced. 

In this chapter, similarly to chapter 3, we consider AL = 1 two-body decays and inverse de- 



cays which act as source and damping terms, and whose expressions are given by eq.(3.40). We 
further include AL = 1 2-to-2 scatterings involving the third generation of quark doublets and 
top quark singlets, whose contribution to the washout term is given by 

where the precise expression of these scatterings is given in appendix B. Finally the total washout 
reads: 



Wscatiz)=2Ssiz) ( ] +4Stiz) , (5.11) 



W{z) = Wid{z)+Wscat{z). 



(5.12) 
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Neutrino-top scatterings also contribute to the production of the lepton asymmetry, the CP 
asymmetry in scatterings being equal to CP asymmetry in decays [731 [7S]. The production term 
is thus: 

P{z) = D{z) + 27,(z) + A^t{z) . (5.13) 

Furthermore, as we want to study the impact of lepton flavours in leptogenesis. Mi should not 
be too heavy. Mi < 10^^ GeV. Therefore, as explained in chapter 3, for values of Ki between 
0.01 and 100, we can safely neglect off-shell AL = 2 while the on-shell contribution is included. 



5.2 Effect of lepton flavours in the type I seesaw 

The key parameters for leptogenesis, which are the couplings of light neutrinos to the heavy 
states A and the mass of the right-handed neutrinos, are all we need to know to evaluate Yb- 
Going in the basis where charged leptons and RHns mass matrices are diagonal, we use the 
Casas-Ibarra parametrisation [51] for the Yukawa couplings: 

= (^/M.i?.^.C/t) , (5.14) 

where the PMNS matrix U and the matrix R are defined in chapter 2, section 2.5.1. We recall 
that the type I seesaw yet contains 15 independent degrees of freedom to fit a single observable: 
the baryon asymmetry. Even though we cannot hope to determine the would-be high-energy 
masses and mixings, we can nevertheless try to shed some light on the optimal regions of the 
parameter space. Such a work has been done in the context of flavoured leptogenesis in [57j, and 
we do not carry here such a detailled analysis. 

Using the above parametrisation, let us see how the inclusion of lepton flavours influence CP 
asymmetries and washout factors. 



5.2,1 CP asymmetry 

Assuming a strong hierarchy in the heavy neutrino sector, and using the parametrisation of 
eq.(5.14), eq.(5.6| can be recast into 

3 MiI™{EAp"^y^"^P^^^a/3^"P^i/3^ip} 



Summing over the lepton flavours, the total CP asymmetry in decays is found to be 



(5.15) 



3 MiIm{E,3"i^^?/3} 



167r Lj(3Tnp\Rm 

no longer dependent on the low-energy phases of the PMNS matrix. However, these phases appear 



with a factor oc m\^2 or ^ind so are suppressed compared to the dominant term oc 

Thus we expect them to play a subdominant role in the general case of a complex R matrix, or 



m 



3- 



5.2.1 CP asymmetry 
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if light neutrinos are very hierarchical. On the contrary, if CP is violated only in the low-energy 
sector, i.e. if R is real, we see that whereas ei vanishes, individual CP asymmetries do not [71] . 
However it has been found that relying only on low-energy phases, a viable leptogenesis requires 
a sizable amount of fine-tuning [SB] . 

In general, the CP asymmetries will receive contribution from both high and low-energy sector^ 
An important point to notice is the modification of the Davidson-Ibarra bound on the CP 
asymmetry. Indeed, in the one flavour approximation, when summing over lepton flavours and 
assuming a strong hierarchy between right-handed neutrinos, the CP asymmetry is bounded by 



3 Mi(m3 — mi 
IGtt V 



2 x/3(mi,mi). (5.17) 



Hence, for degenerate neutrinos, mi ~ ms, and thus the CP asymmetry scales as m^ ^. 

When lepton flavours are included, the bound on individual CP asymmetries is modified. Indeed, 



in the limit of hierarchical RHns, eq.(5.6| reads 



3 Ml 



Im 



{a«iA^i (mt)^^} 



"a— 2 f\f \\ • (5.18) 

Then, writting Aq,i = Aq, where A = |A| and <I> stands for the phase of A, the above expression 
reads 



3 Ml 



The washout factors are given by 



^ Im U^^p (mt) I A, \p , (5.19) 



Ml 

Using maj3 = (t/^dmC^^)^^, and since 




^rnkXplTO-i^a^pUakU^k} < max(m) A^ , (5.21) 

k,p 

we obtain that the individual CP asymmetries are bounded by |75| : 

lb TT f ^ y Ki 

This is an important result, since in the limit of degenerate light neutrinos, instead of scaling 
as m~^ like the total CP asymmetry, individual CP asymmetries scale as m hence are not 
suppressed. We study in fig )(5.1) the validity of the bound of eq.(5.70| by plotting Ea/e^^^ as 



a function of mi. As this bound is derived assuming a strong hierarchy for the right-handed 
neutrino masses, we plot Ea/e™^^ for different ratios Mi/Mj = r = 5,50. We see that even for 
r = 5, the bound on Sa can be safely used. 
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Figure (5.1): Scatter plot of £ct/e™°^ as a function of mi, for different ratios of heavy neutrino masses. 
These plots have been obtained by a random scan of the parameter space, imposing the perturbativity 
constraint A^j < \/47r. M\ ranges from 10^ GeV to 10^^ GeV, and m\ from 10~^eV to 1 eV, even if the 
upper limit is already above the cosmological bound. 
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Figure (5.2): Scatter plot of e\ (left panel) and (right panel) as a function of the mi. Mi is set at 
IQi" GeV, and we take r = 100. 



The different behaviour in the degenerate regime between the flavoured case and the single flavour 



approach can be seen in fig (5.2) , where we display e\ and function of mi. 

In the degenerate regime when mi > matm — 5 x 10~^ eV, we clearly see that while the total CP 
asymmetry is suppressed and scales as m]"^, such a suppression does not occur for individual 
CP asymmetries. On the contrary since in the degenerate regime oc m3 ~ mi, the flavoured 
asymmetries indeed increase with mi. 



^cf. e.g. IBS]. 
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5.2.2 Washout factors 

Using the Casas-Ibarra parametrisation, we can rewrite the washout factors as 

«a = (5.23) 

Then, while the total washout parameter reads 

Ki = A;i|c3C2|^ + A;2|c3S2|^ + fcalssp , ki = nii/m^ . (5-24) 

The dependence of Ki on ms being fealssp, the dependence of, for instance, the electron-flavour 
washout Ke on ms is given by sf^lsse^^]"^ and is very suppressed given the bound ^13 < 0.05. 
Then, besides the fact that the individual washouts are smaller than their sum Ki, one has in 
general an interesting non-democracy of washout strengths. Indeed, it is possible that one flavour 
is only weakly washed out while still having Ki ^ 1. Furthermore, due to the orthogonality 
of the matrix R, we see that rhi is bounded from below by mi. In general such bounds on 
rha do not exist. Then one can have rha ^ mi, which is particularly interesting in the strong 
washout /degenerate regime. 

In this sense, one somehow evades the constraint of the one-flavour picture, where the total 
washout is upper bounded by Ki <20 x (ei/5 x 10"^)°-^^. 



5.2.3 QuaHtative picture 

The rough picture of flavoured leptogenesis follows the same lines as the single flavour approxi- 
mation. First, a thermalisation period during which a lepton asymmetry is produced in hollow. 
Then, after RHns reach thermal equilibrium, decay processes wash-out the lepton asymmetry. 
When the Universe has cooled down to about Mi/lOO, the different processes are frozen, and 
the lepton asymmetry remains constant until being converted by sphalerons. 
The difference with respect to the single flavour approximation concerns how the different lepton 
flavours are coupled to the decaying RHn, that is the rate at which they are produced or washed 
out. We illustrate this difference in fig, (5.3)[ by plotting contours of the baryon asymmetry ob- 
tained by solving the set of EEs in eqs.(5.2), normalised to the observed value Y^^^ = 8.7 x 10^^^, 
in the case where flavours are included (left panel) or not (right panel). Leading to this plot 
we did not use the seesaw formula for the Yukawa couplings, since then we do not control the 
different washout parameters. Instead we uses the Kq,s as input parameters, setting the CP 
asymmetries to their maximal values. We further choose Mi = 10^^ GeV, so that we work in a 
two flavour scheme, and set for convenience mi = 10~^ eV[^ 

In this figure we clearly see the influence of having included flavours. While in the single flavour 



Recall that fhi > mi, hence for Ki — 10 , mi has to be lighter than the quoted value. 
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Figure (5.3): Contour plot of Y^"™/yj''^ for Ke+^ and Kr varying from IQ-^ to 10^. The left (right) 
panel represents the case with (without) flavours. The different contours stand for a ratio of 1, 2 and 4 
and are respectively depicted in black, long-dashed red and dashed blue. 



case the total washout is required to be 0.1 < < 10, in the flavoured case we see that a total 
washout Ki = Ke+^ + Kr as large as 10 — 50 is still possible, if one of the individual washouts is 
close to its optimal value Kq, ~ 3. 

5.3 Evaluation of the baryon asymmetry 

To begin with, we neglect the effects of the flavour conversion due to sphalerons |M|, typically 
0(0. 1), which thus have a subdominant effect. These terms will be studied later on. The EEs 
for the number density then reads: 



dz 

dYAajz) 

dz 



-Ki P{z) ^N{z) , 
-e« Ki P{z) AN{z) - W{z) Y^^{z) , (5.25) 



where AA^(x) = Yn{x) — Y^''{x) and we introduce Ra = —Aaa^a (recall that diagonal entries 
of the conversion matrix A are negatives). The formal solution of this equation for the Aa 
asymmetry is 

YAa{z) = -SaKi r dxP{x)AN{x)e-'^'^-f^'^y^'-yK (5.26) 

As in the single flavour case, the Aa asymmetry evolves differently according to the specific 
washout regime. 

All flavours in the strong washout regime 

In the case where » 1, Y]\[{z) ~ Y^^{z), so that AN{x) = —Y^'/Ki P{x) and 

lAa ^ j'^ dx Y^^'ix) e-^" '^y^^y'^ . (5.27) 



5.3 Evaluation of the baryon asymmetry 
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In this integral, since the main contribution to the asymmetry production comes from z > 1, 
the scattering term in the washout can be neglected, and therefore we recover the result of the 
one-flavour approximation, namely 

lAa ^ -Ea X f ^) Y^\zin) , (5.28) 



with Y^{zin) — 3.9 X 10 ^. The efHciency factor for the flavour a in the strong washout regime 
is thus given by: 

^« - • (5-29) 
By summing over the flavours, the baryon asymmetry is found to be 

^B = ~Y.^A^) Y^^\z,^). (5.30) 

This shows the difference with respect to the one-flavour approximation, for which the baryon 
asymmetry in the strong washout regime is 

^ ,J^"'°.i6 ■ (5-31) 

If, for example, ki = 3 and K2 = 7 (Ki = 10), then assuming equal CP asymmetries, the baryon 
asymmetry in the flavoured case is ~ 2.5 times the unflavoured one; for instance if ki = 1 and 
K2 = 14 {Ki = 15), then the enhancement is of one order of magnitude. 

All flavours in the weak washout regime 

In the case where Kq <C 1 for all distinguishable flavours, when AL = 1 scatterings are included 
in both source and damping terms, one expects the B — L asymmetry to scale as Ki Ka- Indeed, 
during the first stage of RHn thermalisation, an (anti-)asymmetry is produced, which is propor- 
tional to EaKi, and is latter washed out by decays of the RHn. However, since this asymmetry 
undergoes washout processes that are oc Kq, the surviving part is ex K-iKq^. 
The evaluation of Yao follows the same lines as in the single flavour case, where one distinguishes 
anti-asymmetry production for z < z^q during thermalisation of the A'', and the asymmetry pro- 
duction at z > Zeq that cancels out with the former one: 



Jo 



~ Ec^Iki r" dx P{x)Y^\x)e-'^'^ '^y^^y^ + r dx Y^{x)e-'^- '^^^^y^ ] (5.32) 

The first term represents the asymmetry produced in hollow during thermalisation. Its eval- 
uation is simplified using the approximate expression for the source term 



P(.)=.|^ + ., K. = ^^K, (5.33) 
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where Kg ~ O.li^i parametrises the strength of the AL = 1 scatterings. In this weak washout 

case, at high temperatures RHn are far below their thermal equilibrium, Y]\f{z < z^q) ^ 
Y^[z < Zeq), and so a similar expression can be found for the damping term 

W{z) = W,,{z) (^1 + , /? ^ + 2 . (5.34) 

The first term can be written as 

FA«-e.^4no rdx(R„T^..(.)) f^ + .)e-^"^^^^^^-(^K^+f)^^ , (5.35) 

l^a Jo \J^1 J 



where we note that 



n, = \Y^^\z^l)=^^^. (5.36) 



The expression in eq.(5.35) can be further simplified since the leading contribution comes from 
z ~ 1, such that the source term in this weak washout regime reads D{z) + S{z) ~ Ks/Ki. 
Given that 

/•oo 

/ dyy^Ki{y)=x^K2{x) (5.37) 

with x'^K2{x) ~ 2 for X < 1, and using the fact that 

dyy^K^{y) = ^-^^a , (5.38) 

one obtains the antisymmetry production during thermalisation: 

Y^^^e^^Ano^(l-e-'t^) e'^^^. (5.39) 

For z > Zeq, when RHns dominate over their equilibrium abundance, and since the washout 
terms can be neglected, the asymmetry produced in the decay can be writen as: 



Ita = r Ki P{x) Yn{x) 



Ea / Yl^{x) ~ -Sano j^O . (5.40) 



Adding the two contributions, the B — L asymmetry is given in the weak washout regime by 

YAa = ^Aa + ^Aa 

~ X (3.7i^i Ka) y^^(zi„) . (5.41) 

In this weak washout regime, the eflficiency factor is therefore found to be: 

< ~ 3.7i^i ka . (5.42) 



5.3 Evaluation of the baryon asymmetry 
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Finally, the baryon asymmetry in this case is 

= - § IZ^" ^ (^•'^^1 '^") ^Ni^in) , (5.43) 
a 

clearly differing from the baryon asymmetry evaluated in the one flavour case: 

Yb "^-(YI ^" j ^ ^1 {Y1 j K'i^in) ■ (5.44) 

In the regime of weak washout, the inclusion of lepton flavours is found to lower the baryon 
asymmetry when compared to the unflavoured case. 

Mixed regime: some flavours weakly washed out, but Ki:$> 1 

This case does not correspond to any of the cases encountered in the single flavour picture. 
We assume that one flavour, labelled P, is only weakly washed out with Kp <^1, while the total 

washout is strong since all other » 1, a 7^ /3. 

This case has the interesting feature of combining the fast RHn thermalisation, owing to strong 
inverse decays or scatterings, and which typically ensure a copious asymmetry production, with 
the fact that the flavour /3, being only weakly washed out, is somewhat protected from the strong 
washouts of the fast RHn decays. Let us see this in detail. 
During the first stage of thermalisation, one has 

Y<^{z) ^ Y^{z) - kf, W{z) Ya(3{z) . (5.45) 

This gives a contribution 

- £ ' dx y^(x) e-'^^ i-''' . (5.46) 

As <C 1 and the RHns reach thermal equilibrium at high temperatures, in this integral we 
can approximate e~'^i^ c^l — hp j W , such that 

dxYNix) J dyW{y)j , (5.47) 

where Zyj reflects the temperature below which the washout cannot be neglected with Zy] > 1. 
Therefore, in a second stage, for z > z^, one has 

Y^piz) ep Y^{z) ^ y>j -ep Ym{z^) . (5.48) 

Adding the two contributions gives: 

dxY^ix) j dyW{y)] . (5.49) 

The washouts freeze after they reach their maxima, which is evaluated to occur at Zyj 2± 2.4, and 
for which the above integrand is 2± OAY^{zin)- Therefore, the B — L asymmetry in the flavour 
/? is given by 

Y^pc^-EpQARpY;,\zin). (5.50) 
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hence rf^ ~ 0.4 K/3. For the flavour(s) a which is (are) strongly washed out, the result eq.(5.28) 
still holds. 



Global parametrisation 

Similarly to the global parametrisation in the single flavour picture, by means of a simple inter- 
polation we can infer the efficiency factor for all washout regimes 



(5.51) 



which provides a good idea of the result, even if some fine-tuning is still required in order to best 
fit the numerical results. As an illustrative example, we consider a two-fiavour case, Yao — 



We plot in fig, (5.4) the efficiency factor for fiavour /3, whose washout parameter varies from 10 



to 10^, while fixing Kq = 100. Hence, the total washout is strong, and eqs.( 5.28|5.50 ) apply. 




Figure (5.4): Individual efficiency factor 77^ as a function of the corresponding washout k^. We impose 
that the other fiavour is strongly washed-out, with Ka = 100. The black; points represent the numerical 
results, while the red line stands for the approximate formulae given in eqs.( 5.28|5.50 l. 



We observe a good agreement between the formulae derived and the numerical results. However 
if we allow the washout of fiavour a to vary, setting for example Kq, = k^, we see in fig, (5.5) 



that the formula of eq.(5.41), which is depicted in blue, is accurate only for k < 0.05. For 



K > 0.01, some fine-tuning of eq.(5.41 ) is preferable to best fit the result. In the plot of fig (5.5) 



we represent in a red dashed line the formula 

Vf3 - 



;i.3i^i47)-o.4 + (o.4K 



(5.52) 



N 1/0.4 ' 

which gives an accurate result for > 0.01. However, this is an important fine-tuning of the 
parameters, and can be thus disregarded. 



5.3.1 Thermal initial population 

The evaluation of the efficiency factor in the equilibrium scenario is very similar to the one 
fiavour case [M]. Indeed, we saw in chapter 3 that in this scenario, washout is not necessary to 



5.3.1 Thermal initial population 
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Figure (5.5): Individual efficiency factor rjfj as a function of the corresponding washout Kp, in the case 
where the washout of the other flavours also varies, with — Kjj. The black points represent the 
numerical results, while the blue line stands for the approximate formulae given in eqs.( 5.28|5.41 1. The 
dashed red line represents the fit of eq.(5.52 I. 



produce a population of RHns. Nevertheless, for the sake of illustration, we represent in fig (5.6) 
the eficiency factor for a given flavour, in the thermal case (dashed-blue) and in the case where 
the A'^s are initially in thermal equilibrium. 




Figure (5.6): Efficiency factor of a given flavour /3, in the thermal case (dashed blue) and in the case 
where the N are initially in thermal equilibrium (red). 



In the case where all flavours are weakly washed out, 

Y^a{z) = r dxY'^{x)e~~^-^^^''y^^y^ , (5.53) 
the exponential term can be neglected, and given that Yj^{zin) = y^'^(zm) one simply has 

YAa^-SaY^^Zin) , (5.54) 

resulting in a maximal efficiency factor r] = 1. 

On the other hand, when all flavours are strongly washed-out, decays maintain RHns in ther- 
mal equilibrium. Therefore, similarly to the thermal case, one can neglect high temperature 
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contribution, and so 



-^<1'(t \p-fic,f^ dyW{y) 



no 

e 



r dxKc,W{x)e-'^'^^^'^y^^yK (5.55) 



The B — L asymmetry is equal to the one derived for the thermal case: 

0.4 



^Aa ^ -£« X Y^''{zin) ■ (5.56) 



We see that the strong washout regime provides a similar dependence on the washout factor, 
regardless of the thermal history of the RHns. This, together with the fact that the atmospheric 
and solar neutrino masses point towards the strong washout regime, 

rh ~ nisoi ^ n ~ 8 , 

ffi ~ matm ^ K 45 , (5.57) 
indicates that the strong washout regime is a more robust and reliable regime. 



5.3.2 Influence of the ofl'-diagonal terms of the A matrix 

So far we have only included the dominant terms which are the diagonal entries of the matrix 

A, since the off-diagonal terms are 0(0.1). Nonetheless, it is worth studying them. 

When considering the dominant order, we saw that each flavour could be washed out very 

differently. As the off-diagonal terms couple flavours among themselves, it is likely that a flavour, 

being suppressed by either a very weak or a very strong washout, may nontheless provide a 

non-negligible contribution to the final baryon asymmetry, having been pushed up by another 

flavour. 

The evaluation of this off-diagonal term has been done in where we showed that while 
preventing a complete washout of flavours in particular cases, the effect of the diagonal terms on 
the baryon asymmetry is less than 20%. Here we illustrate these results. 
The off-diagonal term contributes to the B — L asymmetry as: 

Y^i{z) = K^Y.^o.p f dxW{x)Y^p{x)e-~^-^^''y'^^y^ , (5.58) 

and so does not depend on the initial condition of the right-handed neutrinos. The evaluation 
of the above integral is made difficult by the fact that it implies knowing the dynamic of the 
individual flavours. 

We thus make a few approximations. First, we neglect the off-diagonal terms entering in the BE 
for the flavour (3. This is quite justified, since its effect on the flavour a is A'^g <^ 1. Then we 
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assume that Ya/3{z) ~ Ya^{oo) so that Ya/3 is factorised out of the integral. We then obtain, 
dropping out the sum on /?: 



rod 

' Aa 



Yap (1 - e^-^-) (5.59) 



A 



where 



a 



dxW{x)^\.'i. (5.60) 



Had we only included inverse decays, we would have obtained a = Svr/S. 

The relative factors Ao^pjAaa strongly differ from one flavour to the other. For example, in a 
two-flavour scheme, the relative factors are 

^Aep. - -CSFa, (1 - e-^^''^) 

Yli ~ -O.OSyAe;. (1 - e-^-°) . (5.61) 
If ^ 1; the off-diagonal term is about ~ {Aa/s / A^a) while if <C 1, this term is roughly 

The B — L flavour conversion is expected to have a significant impact in the case where one 
flavour is very strongly or very weakly washed-out. Indeed, flavour conversion acts as a source 
term for the asymmetry. Therefore in the former case, it partly compensates the washout from 
RHns decays, and so we expect a small increase in the asymmetry. In the weak washout case, 
this extra source term further cancels with the asymmetry produced during thermalisation. So, 
in the weak washout case, one expects a depletion of the asymmetry. 



We illustrate in fig (5.7) the effect of including the off-diagonal terms on the efficiency factor. 
We consider a two-flavour regime, with washout varying from 10-2 to 10^ On the left hand 
side, the CP asymmetries of the different flavours are chosen equal and arbitrary, while on the 



right hand side the CP asymmetries are set to their maximum values, cf. eq.(5.70|. 
Looking only at washout effects (left panels), we see that, as expected, the flavour conversion 
affects Aq asymmetries mostly when » 1 or <C 1. The fact that the flavour Ae/i is (much) 
more affected comes from the relative size of the matrix elements of A: A(,+i_i,t ~ 0.2 = AAr^e+^L- 
When realistic values are used for the CP asymmetries (right panels), we see that the global 
behaviour remains unchanged, albeit this further restricts the influence of the off-diagonal terms 
to the case where ^ KajKx is not strongly suppressed. 
What about the baryon asymmetry? 

Individual asymmetries are affected by ±50%, and given that Yb = 12/37^^ Yao, the baryon 
asymmetry is affected up to 20%, Yb being enhanced in the strong washout regime, and depleted 



in the weak washout regime. This can be seen in fig, (5.8 



We see that the flavour conversion does not significantly modify the baryon asymmetry. However 



^For simplicity, we set Kg = so that the washout of Yae+zj is 
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Figure (5.7): Influence of the flavour conversion on the asymmetry in lepton flavours. The influence on 
the flavour r (e + /i) is displayed in the up (down) panel. The contours represent different ratios of 
^otai j ^^chag ^ with either equal CP asymmetries Sr = ee^^ (left panel) or else CP asymmetries set to their 
maxima — (right panel) . An increase in the efficiency is represented by solid lines, whereas the 
dashed lines represent a reduction of 77, with a colour code: black (±5%), red (±20%) and blue (±50%). 

we notice that for some specific configurations, it prevents Yb from vanishing: in the case where 
Ee^j + = 0, and = Kt, without flavour conversion the baryon asymmetry is zero, while with 
flavour conversion one has 

~ 7xl0-2y^^, (5.62) 

which is non zero, even if compatibility with the observed baryon asymmetry may be hardly 
obtained. 

5.4 Constraints from a successful leptogenesis 

We saw in chapter 3 that requiring a successful leptogenesis in the single flavour picture somewhat 
constrains light and heavy neutrinos. It is mandatory that the decaying RHn is heavy enough. 
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Figure (5.8): Influence of the flavour conversion on the baryon asymmetry. The contours represent 
different values of the ratio Y^°*"'/Y^'°^, the enhancement being depicted by solid lines and the reduction 
by dashed Hues, with the colour code: black (±5%), red (±10%), blue (±20%). 



Ml > 2 X 10 GeV, and that light neutrino masses are lighter than ~ 0.12 eV. Concerning 
the different mixing angles and CP violating phases, the only generic constraint derived is 
the necessity that the matrix R contains CP-odd phases, otherwise no CP asymmetry can be 
generated. 

Given the importance of this constraint, and the fact that the inclusion of lepton flavours strongly 
modifies the key parameters of leptogenesis, it is interesting to investigate the possible influence 
of lepton flavours on the above constraints. 



5.4,1 Lower bound on heavy neutrino mass 

Concerning the lower bound on RHn masses, the inclusion of lepton flavours is a mixed blessing. 
On the one hand, when light neutrinos are (very) hierarchical, since the upper bound on the 
CP asymmetry is roughly the same, we do not expect any important modifications to the lower 
bound on Mi. Actually, since individual CP asymmetries are roughly 1/3 of the total one, we 
even expect in general a slight increase of M™". On the other hand, assuming hierarchical 
RHns, when light neutrinos are (quasi) degenerate, we see that while Si becomes suppressed, 
the CP asymmetries in each lepton flavour are in fact enhanced. In this degenerate limit, we 
expect the lower bound on Mi to decrease. Let us see this in detail, for both the thermal and 
the equilibrium RHn abundance scenarios. 
From the bound on each individual CP asymmetry |75| . 



< 



3MimmE 
IGvrw^ 




(5.63) 
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one has 

\Yb\ ~ 1.26xlO"^Vec 



from which a lower bound on Mj\f^ is derived, 

Mn, > 7.1 X 10« GeV f^^) 



1b 



So V K 



(5.65) 



Since the lower bound on Mi is inversely proportional to the efficiency rja, it will therefore depend 
on the thermal history of the decaying right-handed neutrino. In the case where A^i are produced 
by scatterings, the efficiency is maximised for a washout ~ 1, where 770, — 0.2. In the case 
where A'^i are non-thermally produced, the efficiency reaches its mamximal value of 1 for a very 
weak washout Kq <C 1. The lower bound also depends on the alignment of flavours, and in the 
case of democratic washouts one has jHU] : 

^ I 4.1 X 10*^ GeV in the thermal case 
Mn, > < (5.66) 
I 2.5 X 10 GeV in the dynamical case . 

This bound is close to the one derived in the one-flavour approximation, where M^^ > 4.2 x 

10^ GeV in the thermal case and M^^ > 2.1 x 10^ GeV in the dynamical one [55] . 

The parameter space Mi — mi compatible with successful leptogenesis is shown in fig, (5.9)[ where 



we represent the thermal scenario (up panels) and the case where the A^s are in thermal equi- 
librium at high temperature (down panels), both in the flavoured case (right panels) and in the 
single flavour approximation (left panels). 

These scatter plots have been obtained by numerically solving the set of EEs, using the Casas- 
Ibarra parametrisation for the neutrino Yukawa couplings. We choose a normal mass ordering for 
both low and high energy sectors, setting the solar and atmospheric mass differences, together 
with the corresponding mixing angles, to their best fit values, while ^13 < 13°. For the light 
neutrino mass, we impose the constraint mi < leV. We use as perturbative limit A < 1, hence 
the upper-bound on Mi. The hierarchy for heavy neutrinos is imposed M3 = rM2 = r^Mi, with 
r = 10. The low-energy CP violating phases, as well as those of R, are not constrained. 
We scan over 33000 random sets, and select the solutions satisfying 8.7 x 10^^^ ^ < 
2 X 8.7 X 10^^^. When scanning over the different parameters, we tried to target the lower 
part of the curve, with unequal success, hence the lower bound numerically differs from analyti- 
cal estimates, by a factor 2 — 3. 

We clearly see in these plots the different behaviour regarding the washout. 
In both scenarios, the transition between strong and weak washout rapidly occurs in the one 
flavour case at mi ~ m^:, while this transition occurs at slightly higher mi in the flavoured case, 
due to the ^ < 1 matrix elements. Furthermore, we clearly see that washouts are necessary in 
the thermal scenario, with m > 10"^ eV or equivalently Ki > 0.1, while this is not the case in 
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10 ' 10 * 10 10 - 10 ' 1 10 ' 10 * 10 ' 10 - 10 ' 1 

mi(eV) mi(eV) 



Figure (5.9): Mi — rhi parameter-space compatible with successful leptogenesis: bound on Afi, in the 
case of a zero (thermal) initial A^i abundance for the up (down) panels. The right (left) panels show the 
allowed (Mi-mi) parameter space in the case the lepton flavours are included (neglected). The vertical 
lines represent the equilibrium mass m* ~ 1.08 x 10~^eV (dashed-black), the solar mass -y/Aro^ ~ 
V7.6 X 10-5 eV (in green) and the atmospheric mass yj^ml^^ ~ V2.4 X 10-3 eV (in blue). 



the equilibrium scenario. 

Moreover, we see the important impact of lepton flavours in the limit of degenerate light neutri- 
nos, with a re-opening of the parameter space. Accordingly, the lower bound on Mi is turned 
down in this regime, with for example a diminution of one order of magnitude when rhi ~ 0.1 eV. 

Influence of AL = 2 and the "full flavour regime" 

An important point to notice is that the above plot have been obtained neglecting AL = 2 
scatterings, and assuming that the flavours were fully relevant in leptogenesis. These assumptions 
are justified as long as 

Mi/lO^^GeV < ciO.lKa (5.67) 
Mi/lO^^ < ^2 ^-1 ^ (5 68) 



respectively. Here ci^2 are numerical factors that are 0(1). 

If satisfied, eq.(5.68| implies that lepton flavours are fully relevant in leptogenesis as explained 
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in chapter 4, while eq.(5.67) tells us that AL = 2 lepton violating scatterings can be safely 
neglected. 

We wonder about the robustness of the above discussion regarding these two constraints. On 
the one hand, strengthening the constraint of eq.(5.67| reduces the Mi — fhi parameter-space 
by increasing the minimal allowed washout, whereas on the other hand strengthening eq.(5.68) 
lowers the upper bound on Mi. 




Figure (5.10): Influence of AL = 2 and of recohering processes, in the flavoured case (left) and unfiavoured 
case (right). In light red is depicted the case c = 10, whereas in dark red, light blue and dark blue are 
represented the case ci = C2 = c, with c = 1, c— 0.6 and c = 0.1, respectively. 



In order to quantify this point, we illustrate in fig, (5.10) the Mi — rhi parameter-space which is 
allowed by leptogenesis when the above constraints are included. These constraints are strength- 
ened/weakened by modifying the factor ci^2'- we thus vary ci and C2 for the flavoured case, and 
only ci in the single flavour case. For simplicity, we took ci = C2 = c, with c = 10 depicted 
in light red, which would correspond to the case eqs.(5.67 5.681 are satisfied. In dark red is 
depicted the case c = 1, in light blue c = 0.6 and in dark blue c = 0.1. We see how the allowed 
parameter-space dramatically reduces when we strengthen the constraint. 

Since the plots of figs. (5.9)|(5.11) have been obtained assuming c = 0.8, we can say that the 
conclusion of the previous (and following) section applies. 



5.4,2 Upper bound on light neutrino mass 

The reason why the light neutrino mass is bounded from above in the single flavour approximation 
is twofold. First, in the one flavour approximation, when summing over lepton flavours and 
assuming a strong hierarchy for the right-handed neutrinos, the CP asymmetry is bounded by 

.r^ = ^^^^i^^^V^></5(-i^-i)- M9) 

167r 

Hence, for degenerate neutrinos, mi ~ ms and the CP asymmetry scales as m^f^ and is therefore 
suppressed. This suppression could be compensated by increasing the heavy neutrino mass scale. 
However, Mi cannot be increased harmlessly: AL = 2 washouts that are oc Mi must be out of 
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equilibrium, which means that Mi is upper-constrained and so is mmin- 

When lepton flavours are included in leptogenesis, the situation is drastically different, since 
the flavoured CP asymmetries are no longer suppressed. Indeed, the upper-bound on the CP 
asymmetry in a lepton flavour a reads 

3 Ml nis Fk^ 
IQtt y Ki 

no longer suppressed for degenerate light neutrinos. Therefore, the bound 

mi < 0.12 eV 



(5.70) 



(5.71) 



no longer holds. We can see this in the scatter plots of fig (5.11) which have been obtained in a 



similar way as the plots of fig (5.9) It is clear how the inclusion of lepton flavours re-opens the 



parameter space in the degenerate/strong- washout regime. 
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Figure (5.11): Successful leptogenesis: Mi-mi parameter-space, in the case of a vanishing (equilibrium) 
initial RHn abundance for the up (down) panels. The right (left) panels show the allowed (Mi-mi) 
parameter space for the flavoured (unflavoured) case. The vertical lines represent the equilibrium mass 
~ 1.08 X 10^'^ eV (dashed-black), the solar mass v/Am2^i ~ V7.6 X 10-5 eV (in green) and the 
atmospheric mass -y/Arof^ ~ y/2A x 10"^ eV (in blue) . 



Working in the flavoured framework, the authors of have derived an upper-bound on mi, 
under the constraint that flavours are fully relevant in leptogenesis. However this bound roughly 
reads m < 2 eV, which is already above the cosmological constraint on neutrino mass. 
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5.4,3 Constraints on mixing angles and CP violating phases 

Constraining the seesaw parameters from leptogenesis is an ambitious task that requires a sys- 
tematic analysis of the different high energy phases. Such work has been attempted, both in the 
the single flavour approximation e.g. in [HT], as well as in flavoured picture [571 ESI IS2]- Due to 
the degeneracy among the parameters, even if favourable configurations have been identified, it 
has not been possible to clearly constraint the parameters. Hence the necessity to tie the con- 
straints coming from leptogenesis with other (non-)observations, such as lepton-flavour violating 
processes [M]. But even then, high energy parameters would remain unconstrained, and only 
most favourable parameter space regions could be determined. 

5.5 Summary 

In this chapter, we have studied the effects of lepton flavours in the type I seesaw model of 
leptogenesis. Actually, the conclusions to which we are led are not particular to the type I 
seesaw, but rather they are generic conclusions which may apply in all leptogenesis models, 
under the constraint that at the time the asymmetry is created, the interactions involving the 
charged lepton Yukawa couplings are in-equilibrium/faster than the coherent processes. 
However, we made this study in the type I seesaw, and so let us summarise the results we derived 
in this case. 

In the flavoured treatment of leptogenesis, the individual asymmetries Yao are independently 
produced and washed out. For example, the washout factor oc depends only on the coupling of 
the lepton flavour to the right-handed neutrino. This coupling being in general different for the 
distinct flavours, the latter are differently washed out. This non-alignment of lepton flavours is 
particularly interesting in the strong washout regime Ki ^ 1, since it is possible to have flavours 
which are only weakly washed out, allowing to somewhat escape the usual constraint that Ki 
should not be larger than ~ 15. A similar effect is found for the CP asymmetries. 
We showed that in the relevant temperature regime, the CP asymmetries in each lepton flavour 
have to be considered, instead of the total CP asymmetry (summed over the flavours). This 
has mainly two consequences. The first one comes from the fact that the e^s depend on the 
low-energy CP violating phases of the PMNS mixing matrix. Therefore, one could ask, relying 
only on this low-energy phases, whether or not leptogenesis could work. Such a study has been 
done in [SH], where it was shown that while leptogenesis is indeed possible in that case, there is 
however no clear constraint on the seesaw parameters, owing to the unknown i?-matrix elements. 
It has been further shown, in [51], that leptogenesis is actually insensitive to these low-energy 
phases, in the sense that for any values of the PMNS matrix, one can always find a configuration 
of the i?-matrix for which leptogenesis works. 

The other consequences, which have a direct implication for the low energy phenomenology, stem 
from the modification of the upper-bound on the e^s for hierarchical heavy neutrinos. Since the 
CP asymmetries in lepton flavours are no longer suppressed in the limit of (quasi) degenerate 
light-neutrinos, the upper-bound on this mass scale, which was derived in the single flavour 
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picture, docs not hold in the flavoured context. 

We studied the lower-bound on the lightest right-handed neutrino mass, and showed that the 
bound derived in the single flavour case is only slightly modified, and still sets the scale of 
leptogenesis to be Mi > (2 - 3) x 10^ GeV. 

Finally, we can say that while the inclusion of lepton flavours somewhat weakens the constraints 
derived in the single flavour picture, it only significantly affects leptogenesis in the strong wash- 
out regime. Furthermore, the inclusion of flavours does not provide any significant constraints on 
either the high-energy CP violating sector, or the low-energy one. This comes from our complete 
ignorance of the neutrino Yukawa couplings, which are not constrained in the SM type I seesaw. 
However, this is not always the case, as we will see in the next chapter, where we consider a 
grand-unified framework, where the latter couplings are theoretically predicted. 
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Chapter 6 

Leptogenesis in the type II seesaw 



We saw in the previous chapter that leptogenesis in the type I seesaw can explain the observed 
baryon asymmetry under the main constraint that the lightest right-handed neutrino is heavier 
than ~ 4 X 10^ GeV. However, we also observed that while the low-energy sector is somehow 
constrained by low-energy observations, we have failed to constrain the high-energy sector. This 
lack of predictivity originates from our ignorance of the neutrino Yukawa couplings A and the 
masses and mixings of right-handed neutrinos. Indeed, in the SM, these masses are bare mass 
terms, which are thus unconstrained. 

In 5O(10)-based grand unified theories, the RHns are no longer gauge singlets, and thus the 
Major ana mass term can only arise from the dynamical breaking of an underlying symmetry. 
Moreover, due to the unifying picture, the neutrino Yukawa couplings are related to the Yukawa 
couplings of other field(s), and therefore one could think of a predictive GUT framework in which 
to embed the type I seesaw. 

However in 5O(10)-GUT scenarios, the typical mass of the lightest RHn predicted by a type I 
seesaw is Mi < 10^ GeV, and so hardly compatible with thermal leptogenesis 
Nevertheless, we saw in chapter 2 that different realisations of the seesaw mechanism are possible. 
In the (pure) type II seesaw where one advocates the exchange of scalar SU{2)l triplets, a 
mass term 

rriu ~ vlJl (6.1) 

is obtained, the smallness of being explained by the small vev of an SU(2)l triplet, vl <C v. 
Such a small vev could appear quite unnatural. On the contrary, for GUT models that are based 
on SU{2)l X SU{2)r [MIE^; such a suppression appears very naturally, due to a seesaw-like 
relation among the vevs of the different Higgs representations that the model contains. 
In the type II seesa'wQ light neutrino mass emerges from the neutrino couplings with fermion 
singlets (the RHn) and SU{2)l triplets. Since these models contain two sources of lepton number 
violation, it is interesting to look at the feasibility of leptogenesis. 

^Here we label pure type II seesaw the case where only scalar triplets are added, and type II seesaw the case 
where fermion singlets and scalar triplets are added. 
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In this chapter, we investigate the viability of leptogenesis in a supersymmetric GUT model, 
based on a Left-Right symmetric type II seesaw model. This study is based on where 
we complete the previous study of [51] by the inclusion of lepton flavours and heavier neutrino 
effects. We first introduce the GUT framework in which this study is conducted, untill the type 
II seesaw relation is obtained. Inverting this relation will allow us in a second time to infer the 
high-energy parameters (right-handed neutrino masses and couplings...), using a reconstruction 
procedure developped in [M]. The knowledge of neutrino Yukawa couplings and RHn masses 
enables us to discuss the thermal scenario of leptogenesis in this SUSY-GUT framework. This 
scenario has a priori important differences with respect to the type I scenario discussed in the 
previous chapter, mainly coming from the gauging of SU{2)ji, the inclusion of a scalar triplet and 
the fact that we now work in a supersymmetric extension of the SM, with its inherent problems. 
We finally conduct a numerical study of leptogenesis in this framework, discussing in particular 
the influence of the key parameters of our model, and if the requirement of a viable leptogenesis 
yields constraints for the parameter-space. 

6.1 Framework 

This supersymmetric GUT model is based on ^3221 = SU{3) x SU{2)l x SU{2)r x U{1)b~l- 
Our model consists of three 16 representation, which account for all SM fermions plus 3 right- 
handed neutrinos. This 16, which is spinorial, is the smallest complex representation of 50(10). 
Given the product decomposition under 50(10), 



fermion bilinears can be formed by coupling this product to the appropriate Higgs representations, 
either the symmetric lOg and 126s or the antisymmetric 120a. Another solution consists in 
using non-renormalisable operators, which couple fermion bilinears to products of Higgses, whose 
decompositions contain one of the representations above. 

We consider here only symmetric representations, which consist in two fundamental lO's and a 
rank 5 (complex anti-self dual) 126. The inclusion of a (self-dual) 126 is then necessary in order 
to prevent supersymmetry breaking at the GUT scale by a non-vanishing D-term. Therefore, 
the superpotential contains the following Yukawa couplings: 



Due to the symmetry of the chosen Higgs representation, the couplings V^'^ and / are complex 
symmetric matrices. Assuming that only the SU{2)l doublets of the 10s acquire a vev, the 
relation above tells us that at the GUT scale 



16 X 16 = 10s + 120a + 126; 



(6.2) 



Wy = Y^^j l&i 16j 10„ + Yf^ 16i 16j lOrf + ^ 16^ 16^ 126 . 



(6.3) 



md = me . 



(6.4) 



The GUT framework therefore constrains the neutrino Dirac type mass m/) to be equal to the up 
quark type mass m„. On the other hand, the relation = nie among down-quark and charged 
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lepton masses is clearly not satisfactory |100| . and we will correct it later on. 

Under the Pati-Salam group SU{4) x SU{2)l x SU{2)r [SB], the l26 decomposes as follows 



126 = (6,1,1) + (10,3,1) + (10,1,3) + (15,2,2). (6.5) 

In the type II seesaw, we assume that the two scalar triplets of 126, Al = (1,3,1,1) and 
Ar = (1, 1,3, —1) under G3221, respectively couple to light and heavy neutrinos, according to: 

= ]^ri Al U L, + ^ fl AnNf + y«L,if„ iV/ . (6.6) 

A non-zero vev of the SU{2)r triplet then implies a Majorana mass for right-handed neutrinos, 
whereas a vev for SU{2)l triplet implies a Majorana mass for the left-handed neutrinos. 
We suppose here that the breaking of G3221 occurs due to Ar, whose vev vr lies at an inter- 
mediate scale between the GUT and the electroweak scales. This breaking leads to a Majorana 
mass term of RHns, which violates B — Lhy two units. 

After electroweak symmetry breaking, the light neutrino masses receive two contributions, one 
coming from the type I seesaw [5U] : 

^iij = -vlkkM^^Xlj = --^Xik Xjk , (6.7) 

where Vu is the vev of the up-type Higgs doublet, Vu = vsm{P), and one coming from the pure 
type II seesaw [S3J: 

mi\,=VLf/. (6.8) 

The smallness of the type I mass is obvious, since in this model the breaking of B — L occurs 
far above the electroweak scale, vr » Vu- Conversely, the smallness of pure type II contribution 
comes from the suppression of the SU{2) l triplet vev. 

For instance, this seesaw- like relation for vl can be obtained through the addition of a 54 dimen- 



sional Higgs that couples both to the 10s and to the 126. The Higgs sector of the superpotential 
contains 

Wh D /9i54 126 126 + P2IO 10 54 + M126T26 126 + (6.9) 

where pi^2 are model-dependent couplings, the knowledge of which is beyond the scope of this 
chapter. Since under G3221, the 54 contains a bi-triplet 54 D A54 = (1,3,3,0), it can couple 



to the SU{2)l,r triplets of 126 and to the the product of two) Higgs bidoublets 10 D $ = 
(1,2,2,0), 10 X 10 D (1,3,3,0), according to 

D piA54A/jAl + P2^^ A54 +MaAlAr. (6.10) 

Minimising the scalar potential with respect to vl, one obtains the seesaw relation for the SU{2)l 
triplet vev: 

PlP2V^VR 
VL^ ' 



^For the different decomposition, see, for example, ref. jlOI 
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which ensures that the (pure) type II contribution to light neutrino mass is indeed suppressed. 
In general, the superpotential contains many other terms (as indicated by the dots), which can 
contribute to the above relation. Therefore, in general vl depends on many other couplings. 
Nevertheless, we always have 

Ma^ being defined as a function of the relevant couplings. 

Moreover, we assume here that G3221 is Left-Right (LR) symmetric [H], so that the Lagrangian 
invariant under G3221 is also invariant under the SU{2)l <-> SU{2)r exchange 





Al^Ar. (6.13) 

This LR symmetry enforces fi = fR = f, leading to two important consequences: 

• Firstly, the seesaw relation now reads [53] 

m. = -^/-^AriA. (6.14) 
Ma vr 

Given the low-energy constraints on m^, and the GUT relation between the neutrino and 
the up-quark Yukawa couplings, we can extract /. 

• Secondly, the right-handed neutrino masses simply read 

M = fvR, (6.15) 

and so the knowledge of / from the seesaw relation in turn predicts the RHn spectrum, 
a clear gain when compared to the non-GUT seesaw, where the RHn masses are put by 
hand. 



6.2 Inferring the unknown parameters 
6.2.1 Inverting the seesaw 



Starting from the seesaw relation of eq.(6.14|, the authors of [5U |[TU^ have developed a procedure 
which allows to determine the properties of the complex symmetric matrix /, assuming the 
knowledge of rrii, and A. Here we summarise their results, which are the building-blocks of our 
study ^E\- Given the symmetry of A, it is possible to find A^/^ such that A = A^/^.A^/^-^. Hence 

Y = aX + f3X-'^ (6.16) 



eq.(6.14) can be written 



where 
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and 



Ma vr 

The complex symmetric matrix Y is then diagonalised by an orthogonal matrix Oy, which also 
diagonalises X: 

diag(yi, 2/2,^3) = OyY-Ov 

diag(xi,X2,X3) = O^.X.Oy. (6.19) 



The equation (6.16) then leads to a second degree equation for the eigenvalues yi = axi + I5x^ ^, 



which is easily solved: 



= ^ (yi ±sign(Re(y,))^y,' +4a/3^ • (6.20) 
The couplings / are then given by: 

/ = \^I\X. (a-V2)^ = AV2.o^.diag(xi,X2,X3).Ol,. (a-^^)^ . (6.21) 

Finally, the masses of the RHns are obtained after the diagonalisation of /, through an unitary 
matrix C/j 

M, = f,VR, fi= (uj.f.U})^^ , (6.22) 

where the fiS are positive, and by convention /i < /2 < fs- The unitary matrix Uj relates the 
basis where the right-handed neutrino mass matrix is diagonal to the basis where the coupling A is 
symmetric. We place ourselves in the former basis, which implies a redefinition of A"*^*" = UpX°'''^. 
Since there are two possible choices for each eigenvalue Xi, there are 8 different solutions for the 
matrix /, that is for the right-handed neutrino spectrum. These 8 solutions (2" for n generations) 
constitute the "eight-fold" ambiguity [SS] of the LR symmetric seesaw mechanism: for one low- 
energy spectrum, we have eight possible solutions, which can be distinguished through their 
high- energy effects. 

Behaviour of the difFerent solutions 

In the above relations, the contributions from the different seesaws are clear: /3 stands for the 
type I contribution whereas a denotes the pure type II seesaw. 

It is instructive to consider the asymptotic value 4a /3 <C since in this case the ± solutions 
exhibit the different seesaw contributions: 

^ (6.23) 

The " -I- " solution corresponds to the dominance of the pure type II while " — " solutions 
correspond to a dominance of the type I. 
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The different solutions for / are labelled according to the constitutive : for example (+, +, — ) 
refers to the solution built from (x^, .t^, .T3 ). For the /jS, in general one cannot clearly assess 
which seesaw is the dominant one, apart from the /(+'+'+) which clearly tends, in the 4a /3 <C \yi\^ 
limit, towards a pure type II solution. Similarly, in this limit the f^~'~'~^ solution approaches 
the type I case. 

However, we notice that from the asymptotic behaviour of the XiS we can infer the evolution of 

right-handed neutrino masses as a function of the B — L breaking scale vr. 

Indeed, for the " + " solution, as xf ~ yi/a oc vr, one expects at high vr that the RHn mass 

will continuously increase with vr. For the " — " solutions, since x~ ~ — oc I/vr, VRfi ~ 

constant, and one expects that at high B — L breaking scale the " — " RHn masses will also 

become constant. The value of this constant is of major importance for leptogenesis. 

In the opposite asymptotic limit, 4a/3 S> l^/jp, we have 

xf ~ ±sign(Re(yi))vW« (6-24) 

which shows that the type I and the pure type II somehow cancel each other in rrii,. Indeed, in 
this case, one roughly has / ~ A \fWfoL and 

rn^J ~ ^J~c^\ - mi . (6.25) 



In the intermediate regime, type I and II contributions are of the same order. 

We are interested in the dependency of the different parameters on the B — L breaking scale. 

Therefore, in the following, we fix the ratio 

^ = ^ = ^.0.1. (6.26) 

a VR VL VR 

In fact, depending on the different couplings, (3/a{M^) could be either larger or smaller than 
1{vr). Its value roughly sets the relative importance of the different seesaw contributions to the 
light neutrino masses. Modifying this value only accounts for a translation of the solution on the 
Vr axis. Wc choose (3 /a = 0.1 due to numerical instabilities encountered with (3/a = 1 in the 
high Vr regime. 



Fixing parameters 

In the type I seesaw, we saw that there are 18 parameters among which only 4 are known. Let 

us sec how many parameters wc have in our type II framework. 

Wc place ourselves in the basis where charged lepton (and down-tA'pe quark) and RHn masses 
(/) are diagonal. Then, using the GUT relation mo = m„, we have 

A = Uq Amg{yu,yc,yt)-Uq (6.27) 

rriu = J7/.diag(mi, 777-2, m3).t/j , (6.28) 
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where yq is the Yukawa coupling for the g-quark, and rm are light neutrino masses. The 
unitary matrices Uq and Ui are related to the CKM |103| and PMNS [M] mixing matrices. They 
respectively read: 

Uq = Pu-VcKM.Pd, (6.29) 
Ue = Pe.UpMNS-Pu. (6.30) 

At the GUT scale, since quarks and leptons belong to the same 16 representation, they cannot 
be independently rotated, hence the matrices P^ {k = u, i^, e, d). 

However, at low energy Pe^^u can be rotated away, while the matrix Py contains the usual 
Majorana CP violating phases. All these matrices contain 3 phases, among which two global 
phases can be re-absorbed. We parametrise these matrices as follows: 

Pk = diag(e''^^ , €^'^2 , ei<^3 ) ^k = u,v (6.31) 
Pi = diag(e''^i , e^'^a , 1) , /fe = e, d . (6.32) 



Hence, at the GUT scale, we have 7 additional CP violating phases when compared to the type 
I seesaw. This is a welcome feature regarding leptogenesis, even if these extra phases induce a 
loss of the predictivity for our model. 

The input values for the quark masses and for the CKM matrix at Mz have been taken from 
Refs. |104| and |105| . respectively, and subsequently evolved to the GUT scale using the Mathe- 
matica package REAP [TUS] with an effective SUSY threshold MsusY = 1 TeV and tan (3 = 10. 
The light neutrino mass spectrum has been assumed to be hierarchical, with mi = 10^^ eV, and 
the oscillation parameters have been set to the best fit values of |107| . while keeping ^13 only 
upper-constrained ^13 < 13°. Renormalisation group effects induce a multiplicative factor for 
the different Yukawa couplings: for the first two generations of quarks, renormalisation induces a 
factor ~ 0.41, the top quark Yukawa factor being ~ 0.56, while the neutrino sector is multiplied 
by a factor ~ 1.08. 

All in all, we are left with 17 "free" parameters: the light neutrino mass scale, the ordering of 
the neutrino spectrum and the mixing angle ^13; we also have 10 free CP phases. The B — L 
breaking scale can vary from ~ 10^^ GeV to 10^^ GeV; however perturbativity of the couplings / 
restricts this range: for the value a/ fi = 0.1 that we choose, from ~ 10^^ GeV to ~ 5 x 10^'* GeV. 
Finally, there is another constraint coming from the reheating temperature, which we will discuss 
later and which gives us another parameter, Trh which is allowed to lie between 10^ GeV and 
10^^ GeV. Working in a supersymmetric framework, we have to fix the ratio of the low-energy 
Higgs doublet vev Vu/vd = tan(/3). As already said, we choose here tan(/?) = 10 as a represen- 
tative value. 

Considering the large number of parameters, and our limited computation ability, it is not feasi- 
ble to scan over the full parameter range, but we instead fixed some of the parameters. We list 



^In general, the Yukawa couplings we use are linear combinations of those of eq.(6.3 1, and their precise relation 
is beyond the scope of this chapter. 
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in the following table the usual values for the inputs that we used, distinguishing the parameters 
which we usually vary from those whose value is fixed (NMO stands for neutrino mass ordering, 
which is chosen normal ordered). 



vr{ GeV) 


TRHiGeV) 




^13 




10^2 _ 5 X 1014 


5 X 10^ - 10^1 GeV 


10-5 - 10-1 


- 13° 


0,7r/4 



(6.33) 



NMO 


(3/a 


tan(/3) 


't>l,3 




,e,d 


NH 


0.1 


10 












(6.34) 



Correcting the GUT relation mg = 

Since the only source for the Dirac mass term is the SU{2)l doublet of the two 10 Higgses, we 
have ruu = rnr) and me = rud- The latter relation is, strictly speaking, not a problem since the 
Dirac-type mass of light neutrinos is unknown. However, the relation rUe = rud is clearly wrong, 
and must be improved. 

A usual correction is to assume that the SU{2)l doublets of the 126 also acquire a vev at the 
electroweak scale. Doing this, one thus obtains |1UU| : 

rud = vf^X'^ + vUf, 
me = vf^X" -?,vi2^f , 

= <oA"-3<26/, (6-35) 
such that the GUT relations for the masses are 

mD = ruu-^ f 126 / , 

me = nid- 'ivf2Qf ■ (6.36) 

Inserting the first relation in the seesaw formula, one could in principle invert it in the same way 
as previously explained. However, the fact that rue oc / renders the extraction of / a really hard 
task. 

We could assume that an anti-symmetric 120 does the job, but our model does not contain such 
a representation. The remaining possibility is to invoke adjoint 45 representations which couple 
to the 10 and create an effective 120, as can be seen from the following decomposition: 

10 x 45 = 10 x 120 x 320. (6.37) 

Neglecting the higher dimensional representation, we see that both the 10 and the 120 can 
contribute to the fermion masses. However, the 10 equally contributes to leptons and quarks, 
and so we neglect it (should it take a vev, a redefinition of A"^ would absorb it). Therefore, we 
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need the 120 to acquire a vev. Among the different possibilities, the only products which leave 
GsM invariant at the GUT scale after the 45 develops a vev are 

£ D ^(16 X 16)|(i,2,2) X [(1,2, 2)io X (1,1,3)45] 

+ ^(16 X 16)1(15,2,2) X [(1,2,2)10 X (15,1,1)45] , (6.38) 

where we wrote the decomposition under the Pati-Salam group 5C/(4) x SU{2)l x SU{2)r. We 
see that the two terms contribute differently: the first one is oc (1, 1,3)45, and so a vev in this 
direction does not distinguish quarks from leptons. On the contrary, the second term, which is 
oc (15, 1, 1)45, clearly does. 

Considering that only the latter term acquires a vev, at the GUT scale the relation rUe = 
can be corrected as: 

me = Vd {\d - ^^2) , (6.39) 
md = Vd {\d + "^^2) , (6.40) 

where vi^Vd = ((15,1,1)45). Once corrected, the masses of charged leptons and down type 
quarks at the GUT scale are found to be: 

nid = 0.94 MeV rus = 17 MeV = 0.98 GeV , 

me = 0.346 MeV = 73.0 MeV = 1.25 GeV . 

However, for a consistent correction of the relation nie = m^, since we only include one 45, we 
should obtain that — ~ m^, which is not the case with the values above. 
Nevertheless, supersymmetric threshold effects arising from the SUSY breaking sector (see, e.g. 
|108| and references therein) allow to solve this problem. Indeed, at the tree level, the down 
quarks couple only to one Higgs doublet, Hd, which is not the case for their scalar partner, given 
that W D fiHu-Hd + YdQ.HdD'^. Once SUSY breaking terms are included, this will provide a 
correction to the down quark masses, according to |108) : 

rub ^ nihil + 6m,tMP)) , (6-42) 

where 

(6.43) 

which is evaluated to be Sm ^ 2%. Once these corrections are added, the bottom quark mass 
is mb{MGUT) — 1.17 GeV, which is just enough to fit me and md with the procedure discussed 
before. 

The two mass matrices differ me / md, and so we introduce a unitary matrix Um, the matrix 
rotating from the basis where me is diagonal to the basis where md is diagonal: 

me = The ^ md = U^.md.Um ■ (6.44) 



104 



Leptogenesis in the type II seesaw 



The precise expression of Um is found after fitting the charged lepton and down quark masses, 
and is given by 



I 



\ 




m (XiTi am ^^r^^^ 
125 t^is, ^23, 



V 








(6.45) 



y is a CKM-like mixing matrix, with three real angles and one complex phase. The different 
values we will use for the phases of Vm are listed below, with the corresponding non-zero phases 
of the matrices P^. 



Sets 


nm 
^12 


am 
^13 


^23 




0r 


4>r 


-/-^ 


4>t 




7r/4 


1 


1.07 


0.22 


0.21 


5.80 


3.21 


4.37 


5.86 


0.87 


6.16 




2 


1.07 


0.22 


0.21 


5.80 


3.21 


4.37 


5.86 


0.87 


6.16 


^2 


3 


0.17 


0.066 


0.29 


0.23 


3.14 


0.54 


0.015 


6.27 


0.0032 




4 


0.28 


0.089 


0.37 


0.062 


3.15 


3.12 


6.03 


2.94 


6.19 





Perturbativity of the couplings 

Now that we have been able to reconstruct the couplings /, let us see how they depend on the 
B — L breaking scale. 

First of all, given the asymptotic behaviour of the XjS, we can say that among the different 
solutions, the case (+,+,-1-) will reach the highest values of /, and so of M, while on the 
contrary the case (—,—,—) will reach the lowest ones. 



We thus plot in fig, (6.1) the evolution of the coupling fi for these two solutions, as a function of 
the B — L breaking scale, also considering different values of the smallest light neutrino mass. 
We assume here a normal mass ordering and put all phases (j)^ and ^13 to zero; furthermore, we 
neglect Um for the moment. 

We observe that the coupling / becomes non-perturbative for vr > 10^^ GeV. From now, we 
therefore restrict ourselves to vji < 10^^ GeV. 



6.2,2 Right-handed neutrino spectrum 

Let us study the right-handed neutrino spectra for the different solutions. 



In fig, (6.2) we plot the masses Mj,i = 1,2,3 for the 8 solutions (±,ib,ib) as a function of 
the B — L breaking scale vr. We choose mi = 10~^ eV, the other parameters being fixed as 
before. Qualitatively, we see that the different solutions behave as expected: (-I-, -|-, -|-) masses 
continuously increase with vr, whereas (— , — , — ) masses must reach a plateau. 
Moreover, a " — " solution for Xi entails that one of the masses reaches such a plateau. The 
precise relation between x~ and the mass Mj which stabilises can be found in the appendix B 
of [M]. 

Finally, we can classify the different solutions using the plateau argument, and group them in 
four pairs, (±, -|-, -|-), (it, — , -|-), (it, -|-, — ) and (it, — , — ) as displayed in figj(6.2) 
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Figure (6.1): Evolution of the coupling / as a function of the B—L breaking scale vr for the case (+, +, +) 
(up panels) and (—,—,—) (down panels), for different light-neutrino masses mi, with mi — 0.1 eV, and 
10"^ eV for the left and right panels, respectively. In green (upper-curves) we depict /3, in blue /2 and 
in red /i. 



• In the first pair (it, -|-, -|-), right-handed neutrino masses keep growing with increasing vr-, 
at least until / ~ 1 is reached. For vr < 10^^ GeV, we have Mi < M2 <C M3. A priori, 
one could think of A^i and N2 as potential candidate for leptogenesis. 

• In the second case (±,—,-1-), the lowest RHn mass rapidly reaches a plateau, for vr > 
10^^ GeV. The constant value is Mi ~ 10^ GeV, closed to but still slightly below the type 
I leptogenesis bound Mi < Mf "'^ ~ (2 - 3) x 10^ GeV. We notice that Mi ~ M2 at low 
Vr. Here also, A''i and are "leptogenesis-friendly", albeit Mi seems too light. 

• The third pair (it, -|-, — ) exhibits a very different mass pattern: above the B — L breaking 
scale considered. Mi ~ 10^ GeV, while M2 varies from 10^ GeV to GeV for lO^^ Q^y < 
Vr < 10^^ GeV. Thus A''i cannot be of any help for leptogenesis, while N2 clearly is. 

• In the fourth type of solutions, (it,—,—) cases, both Mi and M2 reach a plateau. Mi 
stabilises around 10^ GeV, and thus A''i is too light for leptogenesis, while A''2, reaching a 
constant mass value M2 ~ M™™'^ provides an interesting possibility. 

Therefore, in the following, we will focus only on four typical solutions: -|-, -|-) which cor- 
responds to the pure type II in the high vr limit; (—,—,—) which conversely corresponds to 
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Case {+,+,+) 
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Figure (6.2): Spectra of right-handed neutrinos resulting from the different solutions for /, as a function 



of vr. The colour code is similar to fig, (6.1) The vertical dashed line represents the perturbative limit 
for /, while the horizontal gray band denotes the lower bound on M, as derived in the type I seesaw (see 
previous chapter). 



the type I solution in the high vr limit; solutions (+, +, — ) and (+, — , +), which are interesting 
mixed type I-pure type II solutions. 
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6.3 Leptogenesis 

Now let us proceed to the analysis of leptogenesis in this scenario, which has already been widely 
studied in the single flavour approximation [M], |109| - |115| . 

The sources of lepton number violation consist in interactions involving right-handed neutrinos 
and scalar triplets. The former are found similar to the type I seesaw, except for the fact that 
the RHns are no longer singlets, since they belong to an SU{2)ji doublet. Thus they interact 
with the additional gauge bosons, which may have consequences for leptogenesis. 
Furthermore, the inclusion of SU{2)l Al triplets induces additional CP asymmetries, as we will 
see below. 

Before examining these points, a remark is in order, regarding the different right-handed neutrino 
spectra discussed above. 

In leptogenesis, it is customary to assume that the contributions of the heavier neutrinos are 
negligible, since the asymmetries produced during decays will be erased by the inverse 
decays which produce A''i. In our case, since M3 » M2,Mi, it is true that the asymmetry 
produced during decays will be washed-out by A''i^2 processes. However, we see that in the 
(±, lb, — ) cases. Ml is too small for leptogenesis, while M2 has potentially the good order of 
magnitude. Hence when evaluating the produced baryon asymmetry, we must consider both A^i 
and processes, while N3 ones can be disregarded. 

6.3.1 Right-handed neutrinos as part of an SU{2)ji doublet 

Having gauged SU (2)/j x U{1) b-l, our model contains four extra gauge bosons ~ B' ,W^, VF^'^, 
which acquire mass when the SU{2)r Higgs triplet Ar develops a vev. Right-handed neutrinos 
are charged under G3221, N'^ = (1, 1/2, 1, 1), and so they interact with the extra gauge bosons 
through different processes |1L)9| I112j : 

• Decays of into Wr, either with on-shell Wr via two-body decays when Mi > Mwji, or 
into off-shell Wr via three-body decays, when Mi < Mwn- 



• Gauge boson mediated scatterings: either (a) right-handed quark-lepton Wr mediated 
scatterings , through s— channel A^e^j, urcIr, and t—, u— channels N dR ^ ^rur 
and Nur <-> crcIr, or (b) lepton-lepton scatterings involving two A^s, via Z' mediated 
scatterings in the s— channel N N"^ ^ grcr, and Wr t—,u— channels crWr ^ N N and 
CR Ni ^ Nj eR. 




Figure (6.3): Two-body RHn decays. 



108 



Leptogenesis in the type II seesaw 




The decay channel of N into Wr may affect leptogenesis by diluting the CP asymmetry generated 
in RHn decays into leptons |116| . However, as Ma ^ Mi we can neglect this dilution. The 
potential danger associated to these new interactions comes from scatterings, which tend to 
keep right-handed neutrinos in- equilibrium with the thermal bath. Two cases may happen: 
either Mi > Mwb, or Mi < Mwj^. It has been shown in |109) that in the former case no lepton 
asymmetry can be generated, unless Mj > 10^^ GeV when the scatterings cr ^ N ^ crW^ 
become out of equilibrium. In the considered case, under the constraint vr < 10^^ GeV, we do not 
reach such high values for the Mj. However, since Mwn — 9rvr/2 ~ vr/3 and Mj = fiVR, this 
case may be only encountered in the -|-, -|-) solution for the heaviest right-handed neutrino, 
for VR ~ lO^^ GeV. 

Thus the danger may only arise from the other case Mj < Mw^, when scatterings mediated 
by the extra gauge bosons force the RHns to be in thermal equilibrium, and so prevent the 
third Sakharov condition from being satisfied. The condition for these scatterings to be out of 
equilibrium is |109| 

Mi<9 X 10^2 GeV x {vr/10^^ GeV)^^^ , (6.46) 
while the condition for the scatterings to be slower than AL = 1 scatterings is 

M. < 2 X 10" K, X (j^)' X {«M (H)' + l.is) , (6.47) 

where Ki is the washout factor of the RHn Ni (see below). 

These two conditions are fulfilled for both Ni and A^2, for the B — L breaking scale we considered 
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vr 10 GeV and for light neutrino masses mi < 0.3 eV. Actually, this constraint may only 

apply for heavy RHns, Mi > 10^^ GeV, for which one encounters an even stronger constraint 
which stems from the gravitino problem, as we will see hereafter. For that reason, we can safely 
neglect the effect of the SU{2)ji x U{\)b-l gauge bosons. 



6.3.2 Scalar triplets in leptogenesis 
Additional CP asymmetries 

The first thing to notice is that contributes to the leptonic CP asymmetry in two different 
ways Since A^ couples to leptons and to the Higgs doublet, a CP asymmetry can be 

generated during A^, decays through the interference between the tree level A^ — > LjLj and the 
vertex correction involving right-handed neutrinos, depicted in fig ](6.5) 



Nk 



Figure (6.5): Vertex correction for decays. 

These decays violate lepton number by two units. However, since we have Ma^ — vr/S, the 
triplet is much heavier than the lightest RHn, Ma^ » Mi. Therefore the asymmetry produced 
during these decays will be washed out by the subsequent processes involving N2 and Ni, and 
so we can neglect all processes with on-shell A^s. 

The other possible contribution comes from the vertex correction of the right-handed neutrino 
decay diagram, with an off-shell A^ running in the loop: 



H /\ H* 




Figure (6.6): Additional CP asymmetry for the type II seesaw. 



As the couplings of A/, to lepton and Higgses are both complex, they induce a new interference 
term with the tree level diagram N — > The corresponding CP asymmetry in lepton flavours 
reads [rnifTTT] : 

// ^ _3_ >p ImjA^fc {mil i) kg \*^] Mj jj 
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where the loop function is given in the MSSM by I113| : 



f'{y) = xlog('i±^V-^l- (6-49) 



6.3.3 Supersymmetric leptogenesis 

The supersymmetric version of thermal leptogenesis is not very different from the SM one |118^ 
[55] . Indeed, given the energy scale considered, ~ Mj » ITeV, the soft SUSY breaking terms are 
negligible, and the masses of fermions and sfermions are equal Mjy = Mj^. Furthermore, from 
the Yukawa couplings in the superpotential 

W D Xi^Lo^.HuNf , (6.50) 

T{Ni ^ = T{Ni ^ l^Hl) = T{N, ^ l^Hu) = r(iV, ^ = ^ |A,„|2 , (6.51) 

lovr 

T{Nf ^ I^&u) = r(iVf ^ i^Hu) = ^ \Ka? . (6.52) 
Supersymmetry implies the equality of the total decay widths, Fat^ = = Mj(AA^^)jj/47r. 

i 

Hence a lepton asymmetry is created through decays of right-handed neutrinos and (RH) 
sneutrinos N in an amount equal to that of the slepton asymmetry 

£Ni,La = ^JV=i,Lc« ^ ^NiXc ^ ^N'^i,La ■ (6.53) 

The CP asymmetry generated in A', decays is [M] : 

Since we have an additional scalar trilinear coupling LHuN'^, the vertex correction is modified 
when compared to the SM well as the self-energy correction, since twice more particles 

are running in the loop. The type I correction thus reads in the MSSM |61| : 

/.(-) = V^f^-logf^lV^-^. (6.55) 



1 — X \ X I I Jx 



As we saw in chapter 3, when RH neutrinos are almost degenerate, the self-energy correction 
should be treated with care. As we can see in fig |(6.2) for example in the (-I-, — , -|-) at low f j?. 



we have such a (partial) degeneracy M\ ~ M2. We therefore use for the CP asymmetry the 
following formula |73| : 

d = ^7>^ El™ {^-^*" [(AAt)i,(C^- + 2C%) + 2 (AAt),, C^^J } , (6.56) 



where C^^ is the vertex correction 



\^x\ r- Mj 



Cl^ = -V^log{ — ],V^=j^, (6.57) 
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the other term standing for the self-energy correction 

In the limit of hierarchical RHns, 2; » 1 and the terms Cg b can be neglected, so that we recover 



the non-resonant formula eq.(6.54). 

Another difference in the supersymmetric version of leptogenesis, is the relation among chemical 
potentials, from which we deduce the conversion between the lepton doublet asymmetry and the 
B — La asymmetries, 

YLM = Y,AapY^P ■ (6.59) 

For Ml < 10^ GeV x (1 -|-tan/?^) , tau- and muon- Yukawa couplings are in equilibrium, thus 
the three flavoured asymmetries lAe,^A/x,^AT are distinguishable. In the MSSM, the A matrix 
is given by [55] : 

/ -93/110 6/55 6/55 \ 
A= 3/40 -19/30 1/30 . (6.60) 
V 3/40 1/30 -19/30 / 

For Ml between 10^ GeV x {I + tan/?^) and 10^^ QgV x (1 tan/?^), only the tau- Yukawas 
are in equilibrium, thus the lepton asymmetry is projected onto a 2 flavour-space (Y^e+n = 
^Ae-i-A/x, ^Ar)- The conversion L ^ B — L now reads: 

A = ( ] . (6.61) 

y 46/761 -494/761 J 

Finally, if Mi is above 10^^ x (1 + tan/?^) GeV, none of the interactions involving charged 
lepton Yukawa couplings are in equilibrium, and we recover the flavour-independent treatment 
of leptogenesis with A = — diag(l, 1, 1). 

Finally, with two light Higgs doublets, the baryon asymmetry resulting from the fast B + L 
violating processes is [S3]: 

^^ = ^E^A. . (6.62) 

a 

With our choice of tan/3 = 10, the 3 flavour regime will generally apply, since in most of the 
solutions ((— , — , — ), (-I-, -I-), etc..) the RH neutrino masses are below 10-*^^ GeV. 
Let us turn to the discussion of the major problem a supersymmetric scenario of leptogenesis 
faces: the constraint on the reheating temperature of the Universe, Trh- 



6.3.4 Reheating temperature and the gravitino problem 

The standard scenario of cosmology relies on the existence of an inflationary period [BdU], which 
explains, among other issues, the observed homogeneity and isotropy of the Universe. The in- 
flationary epoch lasts until the field responsible for driving inflation, the inflaton, reaches the 
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minimum of its potential. At this time, the energy density of the Universe is dominated by 
the vacuum energy of the inflaton. As this field oscillates around its minimum, it reheats the 
Universe by a release of its potential energy. Then, when T ~ Trh-, the inflaton field decays, 
creating degrees of freedom, and so the temperature drops down. 

Therefore all species existing before the reheating of the Universe will be completely diluted by the 
entropy production. Similarly, any pre-existing lepton asymmetry will be washed-out, so that it is 
mandatory that leptogenesis takes place after reheating Tiepto ^ Trh- As Tiepto ~ ^'h ~ 10^ GeV 
thermal leptogenesis requires Trh ^10^ GeV. 

On the other hand, such high temperatures can be problematic when one comes to local super- 
symmetry |119| . Indeed, in such a case, the SUSY partner of the graviton field, the spin 3/2 
gravitino, imposes stringent constraints on the reheating temperature |120| . The gravitino only 
interacts gravitationaly, thus very weakly. Essentially two situations can occur. The lifetime of 
the gravitinos can be larger than the age of the Universe, in which case they might either over- 
close the Universe or their relic density may give an overabundant dark matter density |121| . In 
the opposite case, gravitinos can be unstable, and their decays can either spoil the success of Big 
Bang Nucleosynthesis |122) . or their decay products, if stable, can also provide an overabundant 
dark matter relic density. 

The constraints on the reheating temperature come from the fact that the thermal production 
of gravitinos, which happens during or just after the reheating stage, is more efficient for high 
reheating temperatures. Indeed, the Boltzmann equation for the gravitino number density 713/2, 
in a bath dominated by the inflaton field and radiation, reads [123^ I124j : 

^ +3F(r)n3/2=7(r), 



where pR = n'^g* is the radiation energy density, and pip the energy density of the infiaton 

scalar field. The decay width of the infiaton field is related to the reheating temperature by: 



where in the MSSM the effective number of degrees of freedom is = 228.75. The thermal 
production rate of gravitino is roughly |124| 7(r) ~ 10 T^/M^^. 

Since gravitinos interact only gravitationnaly, one can consider that the density after reheating 
remains unchanged until gravitinos decay (if ever) , such that one has 

Clearly, thermal production of gravitinos is more efficient for high reheating temperatures. 
Now, if gravitinos are stable, they might either overdose the Universe, which leads to the con- 
straint TfiH ^ 2 X 10^'^ GeV, or if they provide the dominant contribution to dark matter, their 
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relic density might be greater than the inferred density ^dmh^ = 0.1143 ±0.0034 [10]. Assuming 



that gravitinos provide the sole contribution to dark matter, one has 

"^3/2 



Orf^/i^ = 1^3/2/1^ ~ 2.78 X 10^° 



3/2 



100 GeV 



(6.65) 



This relation strongly constraints the reheating temperature, as can be seen in fig, (6.7) where we 



plot the relic density of gravitinos (assumed to correspond to the total dark matter abundance). 
In this plot, we further assume that gravitinos are only produced via thermal scatterings (l24j . 
which may not always be the case. Indeed, assuming that gravitinos are the lightest supersym- 
metric particles (LSP), decays of the next-to-LSP may give rise to a non-thermal production of 
gravitinos. However, we see that a stable gravitino implies Trh < 10^° GeV, which is slightly 
above the leptogenesis bound. 



ndmA^= 0.1143+0.0034 




10 10= 

'"3/2 (GcV) 



Figure (6.7): Constraint on the reheating temperature from gravitino abundance, assuming ^Idm — 
The plot has been obtained assuming universal boundary conditions for the gaugino masses, with TOi/2 — 1 
TeV (green), 500 GeV (blue) and mi/2 — 150 GeV (red). For each case, upper and lower curves 

represent — 0.11435 + 0.0034 and ^^/2 = 0.11435 — 0.0034 respectively, while the central value is 
indicated by a grey line. Plot adapted from refs. |124) . 



As already stated, it is also possible to have unstable gravitinos. Nonetheless they can be long- 
Hved particles and decay late. If these decays happen during or after BBN, T < Tbbn, the 
entropy released might jeopardise the successful predictions of BBN |122| . As the decay rate of 
gravitinos is roughly 

3 

r3/2 ^ ^ , (6.66) 

gravitinos decay after BBN if ^ H{Tbbn)- Assuming this is the case, the gravitino abun- 
dance is required to be small enough to not affect BBN, implying an upper-bound on the reheating 
temperature. For 7713/2 ^ 1 TeV, one requires Trh ^ 10^~^ GeV. This value is clearly incompat- 
ible with thermal leptogenesis. If one has 1 TeV < 7^13/2 < 50 TeV, then the constraint roughly 
reads Trh < 10^ GeV. 
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Assuming even heavier gravitinos, one evades the BBN constraint, but then gravitino decays will 
result in a non-thermal production of the LSP, which is assumed to be stable and the dominant 
(if not the only) contribution to the observed dark matter abundance. Considering only this 
non-thermal production for the LSP, its number density is related to the gravitino one by 

^iLSph' 2.78 X lO^o X Yy, ( ^J^^^ ) , (6-67) 

where Xi is the lightest neutralino (in this case the LSP). Under the constraint that ^lsp ^ ^dm, 
reheating temperatures below 3 x 10^*^ GeV are required. 

To summarise this discussion, we can say that even if the gravitino abundance puts stringent 
constraints on the reheating temperature, leading in some models to a clear incompatibility with 
the scenario of thermal leptogenesis, there are however possibilities in which the upper-bound 
on the reheating temperature is sufficiently high to allow for a successful leptogenesis. 
Therefore, even if aware of this caveat, in the following we study the thermal scenario of lepto- 
genesis, assuming that Tr/^ = 10^^ GeV. This high reheating temperature will already allow us 
to distinguish the solutions which provide enough baryon asymmetry from those which do not. 
Then, in a later stage, we will investigate more realistic values for Trh- 

6.3.5 Boltzmann equations 

As before, the evolution of the comoving number densities is obtained by solving the set of 
Boltzmann equations. For the RH (s)neutrino one has: 

Yl^iz) = -2i^, (A(^) + S,(z)) (y^^(z)-y^^(z)) , (6.68) 

where z = Mi/T. In this equation, Yat. is compared to the number density in thermal equilib- 
rium 1^ , which is given by the following equation: 

r-(.) = l|gi«|.^A-.(fl..) '■ii'- - 1.8 X 10- . (6.69) 

where we have introduced Ri = Mi /Mi. 

In the MSSM, the number of degrees of freedom in the thermal bath is g^, = 228.75 (without RH 



''Assuming a Maxwell-Boltzmann (MB) distribution for the equilibrium abundance gives F^' = Y^''. Taking 
Fermi-Dirac (FD) or Bose-Einstein (BE) distributions, the abundances differ at high temperatures by a spin 
factor j66] 

n^^ (T » M) = I n^^ (T » M) = ^ nA„(T » M) , 



and are equal at low temperatures. So even if we assume a MB statistic for the distribution functions, we 

./■eq _ 3-yeq 



nevertheless correct the abundances by the high temperature expansion, and therefore we take Y^'' = jY^'', with 



Y^'' given by eq. (6.69 I. 
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neutrinos or triplets). The individual B — L asymmetries Yaq are driven by: 



i=l,2 

+ 2 J] J];T^i(z)yl„^yA^(z) . (6.70) 

i=l,2 13 



The washout parameters 

^ia ell's ctS usual given by: 



H{Mi) 

but with a slight modification of rh and m*: 



(6.71) 



m, = ^ ~ 1.56 X 10-W . (6.72) 

3^ Mpi 

The total washout parameter Ki is obtained by summing over flavour indices: 

Ki = yKi^ = ^c^ ^ . (6.73) 

^ 1.56 X 10-3ey ^ ^ 

The thermally averaged decay rate Di are given by: 

Di{z) = zRi . (6.74) 

The AL = 1 scatterings are Higgs-mediated processes involving top quarks and antiquarks, and 
receive contributions from the s— and t— channels: 

S^{z) = 2-fl{z)+4f,{z). (6.75) 

The washout term Wi{z) = W-'^{z) + Wl{z) results from the contribution from inverse decays 

Wt{z) = ^ Rf z^ lCi{Riz) , (6.76) 

and AL = 1 scatterings 

A careful study of the reheating problem should be done by considering the influence of the in- 



flaton field on the expansion rate of the Universe, similarly to eqs.(6.63), as done in [SF t lHT tlT^ . 
We do not carry such an analysis, but take the reheating into account in an effective way. We 
consider that the inflaton field decays instantaneously when the temperature drops below Trh, 
and thus we consider that the temperature at which leptogenesis starts Tin coincides with Trh- 
For masses Mi < Trh, this is a rather good approximation. On the other hand, for heavier 
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masses, thermalisation processes of Ni are Boltzmann suppressed and for M; > (2 — 3) x T^ih 
the number density of the created RHns is negligible. 

Had we properly included the contribution of the inflaton field, such a suppression of the thermal- 
isation processes would not have occurred; however, as for T > Trh the Universe is dominated 
by the inflaton energy density, the large Hubble expansion rate oc VP</> + Pr + Pn ~ would 
have diluted processes involving RHns. The difference between the naive reheating we consider, 
and the proper treatment similar to eq.(6.63| has been done, for example, in [124] for the grav- 
itino abundance, where it has been found that the naive prescription, albeit twice larger, is of 
the correct order of magnitude. 

As already stated, in a first stage we will use Tin = 10^^ GeV to discriminate among solutions, 
and then we will lower this value to more acceptable values. 



6.4 Results for the four typical solutions 

Now that the framework of our model has been set up, we can examine the viability of thermal 
leptogenesis for the four characteristic solutions we will focus on: 

• Solution (+,+,-1-) for which the pure type H seesaw dominates at large vr, and where 
RHn masses continuously grow. 

• Solution (+, — , +), an interesting type I+H mixed solution where, contrary to M2 (which 
increases with vr), Mi reaches a plateau close to the type I bound. 

• Solution (+,+,—), which is also a mixed solution, with Mi reaching a constant value, but 
with Ml ~ 10^ GeV, far below the type I bound. 

• Solution (— , — , — ), which tends to the type I in the high vr limit. As in the previous case. 
Ml < 10^ GeV. 

From the above spectra, we can infer that for (+, +, +), the lightest RH neutrino will contribute 
dominantly to the lepton asymmetry production, while for (+, — , +), N2 might play a role since 
Ml could be excessively light. On the other hand, for the remaining solutions (±,ib,— ), it is 
obvious that is the (ultra) dominant contribution. 

We first assume that = mg, with a high reheating temperature Trh = 10^^ GeV, and then 
relax these assumptions in order to study their influence. In the next section, we will study the 
dependence on low-energy parameters that are the light neutrino mass, the CP violating phase 
SpMNS and the mixing angle ^13, but for the moment, we fix mi = 10~^ eV, ^13 = 0, put Spmns 
to zero as well as all the other phases, except for 02 = 7r/4. 



6.4.1 Influence of flavours 



We first investigate the influence of flavours on our four characteristic solutions by considering 
fig, (6.8) , where we plot the baryon asymmetry as a function of the B — L breaking scale, including 
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Case (+,+,+) Case (+,-,+) 




10'^ 10''' 10'" 10'^ 10'^ 10 

Vr (GeV) Vr (GeV) 



Figure (6.8): Evolution of the baryon asymmetry as a function of the B — L breaking scale w/j, for the 
four characteristic solutions. In these plots we compared the flavoured case, depicted in solid-black lines, 
and the unflavoured case, in dashed-red lines. 



lepton flavours (solid black line) or not (red dashed line). 

The different behaviours are more easily understood by looking at fig (6.9) where we plot the 
CP asymmetries and washout factors involving the dominant RHn for the various solutions (see 
the caption for the colour code). 

We first notice that for the present choice of parameters, in the flavoured case (black lines) 
(+, +, +) can easilly saturate the observed baryon asymmetry, while (+, +, — ) is only marginally 
compatible. This is easily understood in fig (6.9) These two solutions have in common that the 
dominant RHn strongly couples to the r flavour (in green). Hence an important CP asymmetry 
is created in this direction, even if somehow compensated by the strong washout. For (+, +, +), 
we see that around uj? ~ (2 — 3) x 10^^ GeV, ei^r ^ 10^^ while at the same time the washout 
decreases, dropping from ~ 10^ down to ~ 5. Hence the baryon asymmetry is enhanced, as can 
be seen from the bump in Yb around u/j ~ (2 — 3) x 10^^ GeV. 
Why is then (+, + , — ) comparatively so small? 

It is because the r flavour, even if copiously produced, is also strongly washed-out: for vr > 
(2 — 3) X 10^^ GeV, K2,t — 20. Hence, a washout four times larger than for the -|-, -|-) case 
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Case (+,+,+) Case (+,+,+) 




10'" 10" 10'-* 10'" 10" lO'* 

Vr (GeV) Vs (GeV) 



Figure (6.9): Evolution of individual CP asymmetries „ and washout factors as a function of for 
the four characteristic solutions. Working in a three flavour scheme, we depicted in green the tau flavour, 
in blue the muon flavour and in red the electron one. We only plot the ei_Q,s and Ki_Q,s corresponding to 
the RHn Ni that provides the dominant contribution to leptogenesis. 



requires a CP asymmetry four times greater, and such a value for the CP asymmetry can only 
be reached for higher vr. 



6.4.1 Influence of flavours 
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Finally, as vr increases, Mi or M2 approach Tm = 10 GeV, and so A'"! or A^2 are no longer 
efficiently produced by inverse decays, thus the dramatic fall of Yb- 

We observe a huge difference when we consider the unflavoured case (dashed red line): while 
for (+, +, +) the unflavoured case is close to the flavoured one, particularly in the high vr-£cp 
regime, the (+,+,—) is highly suppressecj^ In the next section, we will concentrate on this 
feature, which corresponds to the case discussed in [5 ^ I126| . 

Let us examine the remaining (it,— ,ib) solutions. For the parameters chosen, the observed 
value yj''* is not reached. Still, the qualitative description carried for the previous case holds, 
up to the different behaviours of the CP asymmetries. Notice also the strong suppression of the 
unflavoured Yg in the (— , — , — ) case, which also nicely illustrates the survival of the asymmetry 
produced by due to flavour effects [5^ that we now discuss in greater detail. 

Survival of leptogenesis due to lepton flavour effects 

As we just saw, for solutions (+,+,—) and (—,—,—) a huge difference exists between the 
flavoured and the unflavoured picture. 

Assuming that all high energy phases and angles are set to zero, except for $2 = '^/"Q ^^'ve 
for example the following CP asymmetries and washout parameters at vr ~ 10^^ GeV: 



(-,-,-) 


Ni,e 






N2,e 




N2,T 




1.1 X 10-^6 


9.6 X 10-15 


5.8 X lO-i'' 


-1.2 X 10-^ 


-6.4 X 10-® 


-3.4 X 10-"^ 




0.04 


17.2 


16.2 


2.3 
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{ + ,+,-) 


Ni,e 




Nl,T 


N2,e 


N2,l^ 


N2,T 




1.2 X 10-16 


9.7 X 10-15 


5.7 X 10-14 


7.0 X lO-'^ 


2.0 X lO-'^ 


2.6 X 10-6 




0.04 


17.2 


16.2 


0.5 


0.2 


3.5 



We see that in both solutions the washout parameters and the CP asymmetries exhibit similar 
patterns. In both cases, the CP asymmetry from A''! decays is very small, and one can neglect the 
lepton asymmetry produced by A'^i. Moreover, in both cases A'^i interactions present important 
differences regarding the washout of lepton flavours: while the electron flavour is very weakly 
washed out [nie ~ 4 x 10-^), the — t flavours are strongly washed out [nifj^^r — 15). 
On the other hand, e, ^ and r couplings to N2 are of the same order. 

Hence, once N2 has completely decayed, the three asymmetries (L^, L^, L-r) are comparable, un- 
til 2; ~ 1, when and Lt are strongly suppressed by A''! washouts. However, since the electron 
flavour is not affected by these washouts, the baryon asymmetry generated by A^2 processes is 

■^Actually, in this case, as in the (—,—,—) one, the unflavoured picture gives Yb — IQ-^**; but due to the 
numerical precision required to generate these plots, some instabilities arise, so we preferred to impose a cut 
around Yb ~ 10"^^. See the next section for a plot with an increased precision to avoid the instabilities in 
(±, ±, — ) solutions. 

^The case considered here corresponds to a non-trivial Um- For convenience, we rather choose the set 1 of Um, 
deflned in section 6.2.1. This choice only affects the quantitative result, and not the discussion. For example, 
while we obtain Yb ~ 1.25 x 10"^°, had we taken Um = 1 we would have obtained Yb — 7.8 x 10 slightly 
below the observed value. 
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not erased by A'^i leptogenesis, opposed to what occurs in the single flavour approximation. The 



effects of {Le.,Ni) orthogonality can be seen without ambiguity in figs, (6.10) Notice also the 
dramatic effect of flavour conversion on the muon and tau flavours, albeit this effect is only 
sub-leading for the baryon asymmetry. 
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Figure (6.10): Evolution of the asymmetries in the (+ + function of z = Mi/T. The 

thick black and grey lines represent the baryon asymmetry computed with and without flavour effect, 
respectively, highlighting the very effect of lepton flavours in leptogenesis. The thin lines represent the 
contribution to the baryon asymmetry of the different flavours: electron (purple), muon (red) and tau 
(green). On the left panel we include off-diagonal terms when solving the BE, whereas in comparison 
those terms are not included in the right panel. This clearly illustrates the muon and tau survival due 
to their non-diagonal couplings to the electron flavour. 



Following [5H], we evaluate those effects in the illustrative (-|- H — ) solution. Given that right- 
handed neutrino masses are very hierarchical, we can address and A^i leptogenesis indepen- 
dently. Neglecting off-diagonal terms of the matrix A, after A^2 leptogenesis the asymmetries are 
given by: 



-4 X 10-1° , 1a/.|^2 ~ -4 X 10-11 



-1 X 10-9 , Yl, 



6 X 10 



-10 



As the Universe cools down, the comoving number densities remain constant until z ~ 1, where 
A'^i washout processes come into equilibrium, damping the asymmetries according to: 



(6.78) 



and yielding the formal solution 



(6.79) 



6.4.2 Influence of mass corrections 
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The first term corresponds to the depletion of the Yaci from A'"! washouts, whereas the second 
term represents the effect of non-diagonal couplings between flavours. The contribution of the 
former to the final asymmetry turns out to be: 

^i«-l^A«,, xe^^""'^-, (6.80) 



and the second term, evaluated in [HU], gives a contribution which is: 



<^E>-j.xA„, ^^ (6.81) 



The above formulae show that only the electron flavour survives A^i washouts, while muon and 
tau flavours are suppressed by a factor ~ 10"^, illustrated on the right panel of fig (6.10) The 
baryon asymmetry, produced from the electron flavour, is: 

^ ^E^A" =^ ^l-19^Ae|^2 ^ -1-25 X 10-^° , (6.82) 

a 

which is slightly above the observed value (notice the excellent agreement with the numerical 
result Yb ~ -1.45 x 10^1°). 

In the one flavour approximation, A^i washouts exponentially suppress the total lepton asymme- 
try, as seen in fig ](6.10) 



Yl^Yl^^^ xe-'-r-K (6.83) 



We further notice the effect of the off-diagonal couplings which prevent the total depletion of the 
/i, r asymmetries: 

^AMr) - 0-12 Yac ^ 0.12yAe|jv2 e^^"''!- ~ -4.4 x 10"" . (6.84) 

This effect is spectacular for the individual lepton asymmetries, but only marginal for the baryon 
asymmetry. 

This example clearly illustrates how the orthogonality of the leptonic directions and the decay 
of A''i and N2 in flavour space allow N2 leptogenesis to survive the washout of A^i. 
In the following, we do not pursue the discussion of the unflavoured picture. 

6.4.2 Influence of mass corrections 

Let us now examine the influence of the correction to the relation nie = rrid. This correction 
entails a redefinition of the down-quark Yukawa coupling — > ula through the matrix 

Um which diagonalises the charged lepton masses. Going to the basis where down-quark masses 
are diagonal, the neutrino Yukawa couplings are affected by Um, as is the right-handed neutrino 
spectrum, as well as the CP asymmetries and washout factors. Using the four sets of Um defined 
in section 6.2.1, we plot in fig ](6.11)" the baryon asymmetry as a function of vr. We do not 
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Case (+,+,+) Case (+,-,+) 




Vr (GeV) Vr (GeV) 

Figure (6.11): Evolution of the baryon asymmetry as a function of the B—L breaking scale vr, for the four 
characteristic solutions, and for the four sets of matrix Um defined in section 6.2.1. The different colours 
represent the different sets, with sets 1,2,3 and 4 depicted in red, blue, green and purple respectively. 



analyse the different Um separately, but just remark their main features. 

The first thing to notice is that for solutions (+, +, ±), the qualitative picture remains unchanged: 
Yb increases with vr, until the right-handed neutrino becomes too heavy and thus its production 
is suppressed, accounting for the fall of Yb at high vr. If we look more carefully, we see that 
these two solutions are enhanced, especially for the set 1 (in red), for which both two solutions 
reach the observed value. 

Interestingly, for the (+,+,-1-) case, A''i and N2 become partially degenerate around ~ 4 x 
10^^ GeV, leading to a resonance of the baryon asymmetry, as can be seen from the peak in Yb 
for this solution. This resonance is shown in fig ](6.12)" where we plot the baryon asymmetry as 
a function of vr and the corresponding right-handed neutrino masses. 

Nevertheless, this degeneracy is accidental and is not a generic feature of this model. Moreover, 
relying on resonances to obtain a sufficiently large baryon asymmetry, while Mi ~ 10^ GeV may 
appear as extremely fine-tuned, and so we no longer consider this possibility. 
A more robust and welcome consequence of having included the mass correction stems from 
the success of the (-1-,— ,-|-) solution, for which the observed baryon asymmetry is reached for 
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Figure (6.12): Zoom on the quasi-degenerate region of the (+, +, +) solution depicted in red in fig, (6.11 



On the right-hand side, the plot shows the evolution of Mi and M2, which become quasi-degenerate 
around ~ 3.8 x 10^^ GeV. This quasi-degeneracy implies a resonance of the baryon asymmetry, as 
shown on the left-hand side. 



10^^ GeV <VR< 10'^ GeV. 

Consider for example the solution corresponding to the first set of Um (in red) of fig (6.11) , which 
has the same high-energy phases (f)"^ = ir/A of that in fig, (6.8) For vr ~ 10^^^^'* GeV, we see 



that the baryon asymmetry is enhanced by two orders of magnitude when the mass correction 
is included. 

This comes from the very different behaviours of the CP asymmetries and washout factors, as 



can be seen by comparing figs (6.9) and (6.13) In the former case the maximum CP asymmetry 
is reached for ~ 2 x 10^^ GeV, and is eir — 10"'' with a washout ki^- ~ 20. In the latter, when 
the mass correction is included, the CP asymmetry is comparatively enhanced by one order of 
magnitude, while the washout factor is reduced by one order of magnitude. 



10-2 „ 



Case (+,-,+): U,„ i= I3 



Case (+,-,+): U„, i= I3 




10'3 10'-* 10'^ 

Vr (GeV) 



10'^ 10" 10'^ 
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Figure (6.13): CP asymmetries and washout factors for the (-I-, — , +) solution depicted in red fig. 
Colour code as in fig (6.9) 



(6.11) 
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The enhancement of the baryon asymmetry coming from the inclusion of the mass correction 
can be analytically explained from the influence of 9^ of the matrix Um- Indeed, we roughly 
have that 



am \ 1 2 
^12 



Ml oc VR{fi\cos{eT3)cos{9T2)\^ + f2\cos{eT3)sm 
M2 oc t;,j(/i|cos(^2"^3)sin(0r2)l' + /2|cos(0"3)cos(er2)P + -) ' (6-85) 

where the dot stands for terms of the same order for Mi and M2. Hence, all other parameters 
being fixed, varying 6Y2 modifies the relative (/i — /2) weights, and so potentially increases Mi, 
which was previously slightly below the type I bound (cf. fig ](6.2) ). 



In the following, we will only consider the favourable case which is depicted in red in fig (6.11) 



6.4.3 Dependence on the reheating temperature 

Let us now discuss the important question in thermal leptogenesis of the lower bound on Mj, or 
equivalently, the lowest possible reheating temperature. In the type I seesaw, a bound is derived 
on Ml, namely Mi > 2 x 10^ GeV, which follows from the constraint of having a large enough 
CP asymmetry. This bound implies a lower bound on Tjih, which must be > 2 x 10^ GeV. In 
the type II seesaw, assuming strong hierarchy among right-handed neutrinos, the individual CP 
asymmetries are bounded by |113j : 



3 Ml 



max 



According to the mass pattern for the RHn we have encountered, this relation is not likely to 
apply, since either the neutrino that dominantly contributes to the asymmetry production is N2, 
with Ml <C M2 <C M3, or else the main contribution comes from A'^i, but with a mild hierarchy 
between A^i and N2. 



Nevertheless, we partially answer this question with fig, (6.14) where, for the four solutions 
represented in the Tjih — vr parameter-space, we display the limits above which the computed 
baryon asymmetry at least equals the observed one (except for the (—,—,—) case, cf. figure 
caption). In order to derive this plot, we fix all other parameters, with mi = 10~^eV, and 
choose the first set for with $2 = the other phases as well as ^13 being set to zero. 
The first thing to notice is the failure of (— , — , — ) to produce enough baryon asymmetry. 
Solution (-|-,-|-, — ) works for Trh > 5 x 10^'^ GeV; for lower temperatures, N2 becomes heavier 
than Trh and its production is Boltzmann suppressed. 

Particularly interesting are the solutions (-I-, ±, -|-), for which the reheating temperature can 
be brought down to more reasonable values: the lower bound derived in these cases is Trh > 
5 X 10^ GeV. 

However, we again want to stress that we do not scan over the entire parameter space, but 
restrict ourselves to the case where only one high energy CP violating phase is non-zero. As a 
consequence, the lower bound on Trh should be seen only as a rough approximation of the true 
lower bound. 



6.5 Dependence on the low energy parameters 



125 




10'' 10" 10'5 

V« (GeV) 

Figure (6.14): Contour lines of the ratio of the computed baryon asymmetry over the observed one as 
function of vji and Tin, for the different characteristic solutions: solution (+ ++) is depicted in thick black 

line, solution (+ H — ) is in thick long-dashed blue, whereas solution (H h) is in thick dashed red. For 

these three solutions, the contours represent a unit ratio, and inside the contour we have y^"^™- /Y^^ > 1. 

We have depicted in a thin black line the solution ( ) for which the observed bound is not reached, 

and so for this solution we plotted the ratio y^"^™ /Y^^ = 0.1. 



6.5 Dependence on the low energy parameters 

As stated many times before, the conclusion of this study crucially depends on the different input 
parameters. The high energy parameter-space, which is 10 dimensional, exceeds our computation 
ability, and furthermore is not implied in low-energy experiments. On the other side, m^, 613 
and 6) are very important low-energy parameters which are involved, for example, in a possible 
violation of CP in neutrino oscillations. In this section we numerically investigate their influence. 



6.5.1 Dependence on the light-neutrino mass 

We first examine the influence of the light neutrino mass scale, assuming a normal mass ordering. 



As can be seen from fig (6.1), where we plot the fiS for mi = 0.1 eV and 10~^eV, increasing 
nil raises the right-handed neutrino masses corresponding to a " -|- " solution, while conversely 
decreasing its value lowers masses which correspond to a " — " solution. This behaviour can 



be seen in the contour plot of fig, (6.15) where we show for our characteristic solutions the 
region in the {m^-^ = mi^vn) plane where enough baryon asymmetry is created. We choose 
Trh = 10^^ GeV, and the second set for Um - ^2- 
We notice the strong influence of mi on the different solutions. 

Furthermore, even if solution (— , — , — ) fails to produce enough baryon asymmetry, its behaviour 
is similar to the one of (-|-,-|-,— ), as expected. These two solutions rely on the decays of N2, 
but also strongly on the washouts from inverse decays producing A^i, as we saw in the previous 
section. The fact is that when mi increases, ki^qS also increase, hence the asymmetry produced 
by leptogenesis will be erased by subsequent A''i processes. For these solutions, we have an 
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Figure (6.15): Contour lines of the ratio of the computed baryon asymmetry over the observed one, as 
functions of and Hght neutrino mass scale mi. The ratio y™'''^ /Y^^ = 1 is fulfilled for the three 

solutions: (+ + +) depicted in thick black line, (+ H — ) in thick long-dashed blue, and (H h) in thick 

dashed red. Inside the contours we have y^^"^ /Y^^^ > 1. In the case of solution ( ), for which the 

observed bound cannot be reached, we have depicted in thin black line the ratio yg^'m jy^^s _ q -|^_ 



upper-bound on mi, namely mi < 2 x 10~^eV. 

In the (-I-, — , -I-) case, when mi increases. Mi, which corresponds to the " — " solution decreases, 
dropping below 10^ GeV. This could be somewhat compensated with the increase of M2, but 
since M2 > 10^^, its contribution is suppressed. For this solution, a similar bound on mi is 
derived, namely mi < 5 X 10-2 eV, with vr > lO^^QgV. 

In the (-I-, -I-) case, increasing mi increases Mi, and consequently the efficiency of leptogenesis, 
until Ml > (2 — 3) X Trh and Yb drops. To avoid this decrease in Yb, lower values of Mi are 
required, and so lower values of vr. 



6.5.2 Dependence on low-energy mixing angles and phases. 

Let us now discuss the influence of the low-energy CP violating phase 6 and mixing angle ^13. 
These two quantities are crucial for low energy neutrino physics, since ^13 = 0° and 5 = 0, vr 
would imply that CP is conserved in neutrino oscillations |127) . 

In order to derive constraints, we choose Trh = 7 x 10^ GeV for (-I-, it, -|-), and Trh = 



5 X 10 GeV for the remaining solutions, according to fig, (6.14) We further take mi = 10 eV 



since this value suits all solutions. 

We first consider in fig, (6.16)| the infiuence of ^13 on Yb, when vr varies, with a fixed CP phase 



5 = 0. 

We see that the infiuence of ^13 on Yb strongly depends on the type of solutions. 
For all the solutions, we notice that for 5 = 0, ^13 = gives the best results. For the (it, it, — ) 
solutions, we observe a huge dependence on ^13: Yg^"'^ decreases by two orders of magnitude 
when 6*13 goes from 0° to its upper-bound, 13°. We can consequently say that 6 = favours 



6.5.2 Dependence on low-energy mixing angles and phases. 
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Case (+++) Case (h — i-) 




Vr (GeV) V« (GeV) 

Figure (6.16): Baryon asymmetry as a function of vn^ for the four typical solutions. In these plots are 
displayed Yg for different values of 6113, taking (5 = 0. We represent in black, purple, blue, red and green 
the values 6*13 = 0°, 2°, 5°, 9°, 13°, respectively. For the top panel, we choose Trh = 7 x 10^ GeV, while 
for the bottom panel Trh = 5 x 10^" GeV. 



013 = 0°. What about non- vanishing values of 5? 

To study the influence of 5, we fix f = 6 x 10^^ GeV for (+, ±, +) and vr = h x 10^^ GeV for 
(lb, lb, — ), keeping the values of Trh quoted above; we then plot in fig, (6.17)| in the ^13 — 5 
plane. 

We see that changing these parameters does not improve the situation for (—,—,—), which still 
fails to produce enough baryon asymmetry. However, we notice that its behaviour is very similar 
to the (+, +, — ) case, for which the observed bound can be reached: in both cases, the baryon 
asymmetry is larger for small values of ^13, with a factor ~ 100 between ^13 = 0° and ^13 = 13°. 
Moreover, we clearly see a preference of these solutions for 5 ~ vr, for which higher values of ^13 
are possible. Actually 5 = tt is the most favourable case for (+, +, — ), while (— , — , — ) seems to 
favour 5 ~ vr — 7r/8. 

This preference for 5 ~ vr is also observed for the (+, ib, +) solutions, but in these cases non-zero 
values of ^13 maximise . 

We can conclude that, for the choice of parameters made, our model points towards non-zero 
values of 5, in fact values close to vr. Hence it implies small CP violation in oscillation experi- 
ments [T?7] . 

However, we stress once again that this conclusion is not general, being based on a partial scan 
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Figure (6.17): Contour lines of the ratio of the computeci baryon asymmetry over observed one, as a 
function of the low-energy parameters 5 and ^13. We depicted in red lines the unit ratio, if reached. 



of the 17-dimensionnal parameter space. 

6.6 Summary 

In this chapter we have studied leptogenesis in the type II seesaw. This model has the interesting 
feature of containing two sources of lepton number violation, originating from the couplings with 
the scalar SU{2)l triplet and with the right-handed neutrinos. 

Moreover, our model is embedded in a SUSY-GUT framework: it has the advantage of constrain- 
ing the neutrino Yukawa couplings, the ignorance of which is one of the weakness of the model we 
studied in the previous chapter. Furthermore, the fact that we assume a Left-Right symmetry 
enables us to invert the seesaw relation, and so to characterise right-handed neutrinos in terms 
of light neutrino and up-type quark properties. For the RHn spectrum, 8 solutions are found, 
which can be grouped in four pairs according to their different properties. We focus on one 
solution per pair, which exhibits either the dominance of one (pure) type of seesaw, for (-I-, -|-, -|-) 
and (— , — , — ), or else an interesting mixing between type I and pure type II, for (-I-, — , -|-) and 
(-I-, -|-, — ). We thus discriminate among this eightfold degeneracy with the requirement that lep- 
togenesis is viable. 

The leptogenesis criterion has already been used by the authors of [5^, but as they already 
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stressed, the study was far from being definitive, mainly for two reasons. Firstly, the treatment 
of leptogenesis was not correct, since it was based on the single flavour picture. Secondly, the fact 
that only a pair of 10 Higgs representations had been included to give mass to light particles, 
imposed the relation rrie = rud at the GUT scale, which clearly does not hold. Therefore in this 
chapter, which is based on [HH], we improved the study regarding these two points. 
Thus, we have included the effects of lepton flavours, inclusion which is not only required [75] , 
but instrumental [7T]. This is especially true in the case we considered here, where usually the 
total washout Ki is strong but the individual ones are not. We further included heavy neutrino 
flavours, at least partially, since the heaviest RHn is usually much heavier than the other two 
RHns. This inclusion follows the fact that in the (it, it, — ) cases, while A''i is too light to be of 
use for leptogenesis, N2 is of the good order of magnitude. 

We found that indeed, adding both light and heavy flavours, improves the situation. The most 
spectacular example is the (+, +, — ) solution, where the asymmetry is produced due to N2 decays, 
and the inclusion of lepton flavours prevents this asymmetry from being completely washed-out 
by A''i processes, as it would be the case in the unflavoured picture. 

Moreover, by adding in an effective manner a non-renormalisable term involving a 45 Higgs 
representation, which is thought to be present in typical S'O(10)-based GUT, the correction 
of the GUT relation rUe = constitutes a clear improvement of the model towards realis- 
tic fermion masses. This correction implies a redefinition of the lepton right handed neutrino 
Yukawa couplings A, which turns out to modify the right-handed spectra, individual washouts 
and CP asymmetries, the relevant quantities for leptogenesis. Owing to some freedom in this 
redefinition, we found that correcting fermion masses greatly improves the results. A typical 
example is the success of the (-I-, — , -I-) case. 

Requiring that leptogenesis successfully predicts the observed amount of baryon asymmetry, we 
thus found that 3 of our 4 characteristic solutions do work. However, these results might be 
pointless, since these might be a conflict with the reheating temperature. Examining this prob- 
lem, we found that the solutions -|-, -|-) and -|-, — ) work in a large part of the parameter 
space with Trh ^ 10^'^GeV, whereas (-|-,-|-, — ) would require Trh ^ 5 x 10^'^GeV. Under the 
constraint that Trh ^ 2 x lO^'^ GeV we thus found that two of our four characteristic solutions 
satisfy the leptogenesis criterion. Hence, we can say that four of the eight solutions are compat- 
ible with leptogenesis. Nevertheless, our model suffers from its lack of predictivity, since the 10 
dimensional high-energy parameter-space remains unconstrained. Therefore, it would be inter- 
esting to link the constraints coming from leptogenesis to other constraints coming, for example, 
from lepton-flavour violation. It could be also interesting to build a complete SO{10) GUT 
model which provides the LR symmetric seesaw, together with the required proton stability, 
doublet-triplet splitting and realistic fermion masses and mixings. 
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Conclusions 



We have seen that under the assumption that neutrinos get their masses through the seesaw 
mechanism, leptogenesis qualitatively occurs; in chapter 3, we discussed the requirements under 
which the single flavour picture can quantitatively work. These conditions are: 

• The matrix R must be complex, and furthermore CP violation has to be large enough in 
order to have ecp ^ 5 x 10^''. 

• The right-handed neutrino masses are bounded: 2 x 10^ GeV < Mi < 10^^ GeV. The 
lower bound comes from the requirement that ecp is large enough, while the upper-bound 
represents the temperature below which electroweak sphalerons come in-equilibrium. 

• Assuming hierarchical heavy neutrinos, the light neutrino masses are upper-consrained: 
mi < O.lSeV. 

This is for the single flavour picture. 

However, we showed [7S] that lepton flavours have to be included in leptogenesis, since for Mi < 
10^^ GeV, the interactions involving the charged lepton Yukawa couplings are in-equilibrium. 
Since the different lepton flavours have different Yukawa couplings, the lepton doublet interacting 
with the RHns has to be decomposed into the different flavours. Consequently, one has to consider 
lepton flavour asymmetries. 

In chapter 4, we show how to include the evolution of flavoured asymmetries in the Boltzmann 
equations, and what are the constraints under which the flavours are relevant. We showed in 
this chapter that, for instance, when the tau Yukawas are in-equilibrium, then the quantum 
correlations between the tau flavour and the other flavours are exponentially damped, which 
is not the case for the correlations between the muon and the electron flavours. Regardless of 
whether leptogenesis is efficient or not, we identified different temperature regimes: 

• If 10^ GeV < Ml < 10^^, only the interactions involving the tau Yukawa couplings are in- 
equilibrium, and therefore 2 flavours are distinguishable: the tau flavour and the coherent 
sum of e -|- ^ flavours. 

• If Ml < 10^ then the muon- Yukawa interactions are also in-equilibrium, and the lepton 
asymmetry is projected onto a three-flavour space, (r, ;U,e). 



132 



Conclusions 



The electron- Yukawa interactions are in- equilibrium for Mi < 10 GeV, that is for masses far 
below what is required by leptogenesis. 

We subsequently applied these flavoured recipes to the type I seesaw in the SM model, and to 
the type II seesaw in a GUT model, respectively in chapter 5 and 6. 

The chapter 5 is based on [751 [711 ED] , and provides the main results obtained during my thesis, 
which can be compared with the constraints of the single-flavour picture above. 

• In the flavoured context, the i?-matrix can be real, with non-zero lepton flavoured CP 
asymmetries, since the latter now depend on the CP phases of the lepton mixing matrix 
UpMNS- However, the lower bound on the CP asymmetries is still e > 5 x 10"'^. A 
noticeable modification, which may provide one of the most important effects of lepton 
flavours, is the modification of the upper-bound on the CP asymmetry in the regime of 
degenerate light neutrinos. Indeed, while in the single flavour picture scp oc and is 
thus suppressed, when flavours are included the upper-bound rather reads ~ m and is 
even enhanced in this limit. 

• Regarding the lower bound on right-handed neutrinos, lepton flavours do not significantly 
modify the picture, with Mi > 2.5 x 10^ GeV. 

• On the contrary, since individual CP asymmetries are no-longer suppressed in the regime 
of degenerate light neutrinos, the allowed parameter-space is enlarged, and light neutrino 
masses are no longer upper-constrained by leptogenesis when right-handed neutrinos are 
hierarchical. 

The main effect of lepton flavours lies in the degenerate regime, which corresponds to a regime of 
strong washout. This is simply because the couplings of different flavours with the right-handed 
neutrino are not equal, and so the flavours interact at different rates. 

Another aspect of leptogenesis, which perhaps constitutes one of its weakness, is the lack of 
predictivity regarding the different CP violating phases and especially for the high-energy sec- 
tor. In the single flavour picture and in the flavoured one as well, only optimal regions in the 
parameter-space can be inferred. 

In order to remedy to this problem, we studied in chapter 6 supersymmetric leptogenesis in 
a Left-Right symmetric type II seesaw model. In this model, the GUT embedding yields the 
interesting result that the Dirac type mass of light neutrinos is no longer free, but related to the 
masses of the other fermions. This is a clear gain compared to a non-GUT scenario. 
In this scenario, the high-energy sector was reconstructed from an inversion of the seesaw formula, 
but nevertheless suffers from an "eightfold" ambiguity in the solutions. We therefore investigated 
if leptogenesis enabled us to discriminate among these solutions. Furthermore, given the spectra 
we encountered, for half of the solutions it was clear that only the second lightest right-handed 
neutrino was relevant for leptogenesis. We therefore studied leptogenesis, including both light 
and heavy flavour effects. In this case, we found that 2 of our 8 solutions survived, namely the 
solution which correspond to a pure type II seesaw in the high B — L breaking scale limit. 
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We further refined our model, correcting the GUT relation among charged lepton and down t^'pe 
quark masses, through the inclusion of non-renormalisablc operators. We found that this step 
towards more realistic fermion masses strongly affects leptogenesis, increasing to 6 the number 
of allowed solutions. We further found a nice illustration of the importance of both light and 
heavy flavours for the case where the lepton asymmetry is generated by the decays of A^2, and 
the inclusion of lepton flavours prevents this asymmetry from being washed-out by A^i processes. 
However, since we were working in a supersymmetric extension, we had to have in mind reheat- 
ing temperature constraints. When investigating whether or not leptogenesis could work with 
Trh ^ 10''^° GeV, we found that 4 solutions were compatible, in a large part of the parameter 
space. We then consider whether the requirement of successful leptogenesis enabled us to make 
some predictions for the seesaw parameters. Unfortunatlj^, the "cure seems worse than the dis- 
ease", since our GUT model contains more free parameters than the non-GUT one. 
This lack of predictivit)' is a common plague to high-energ}^ seesaw models, and thus to the ther- 
mal scenario of leptogenesis. One then one asks: is there a remedy? There are at least several 
ways out: the first one would be to link the different observables which involve seesaw couplings, 
another could be to study leptogenesis in a complete and consistent GUT scenario, where the 
number of free parameters is dramatically reduced or even over-determined. Finally, lowering the 
seesaw scale could provide additionnal observable effects at the LHC or other future experiments 
and is therefore an interesting solution. Nevertheless, lowering the seesaw scale enters in conflict 
with the scenario of thermal leptogenesis that we discussed in this thesis. Hence the need to 
low-energy leptogenesis models. 
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Appendix A 

A few words on particle physics 



A.l The Standard Model 

The Standard Model of particles relies on the local invariance under SU{3)c x SU{2)l x U{1)y 
symmetry, the matter sector being composed of 3 generations of 15 fundamental fermion fields 
which are given in the following table: 



field 


G32 


Qrs 


Qy 


B 


L 


UL 


(3,2) 


1/2 


1/3 


1/3 





di 


(3,2) 


-1/2 


1/3 


1/3 





UR 


(3,1) 





4/3 


1/3 





d-R 


(3,1) 





-2/3 


1/3 





l^L 


(1,2) 


1/2 


-1 





1 


SL 


(1,2) 


-1/2 


-1 





1 


eR 


(1,1) 





-2 





1 



(A-1) 



The SM successfully describes the strong and electroweak interactions thanks to the exchange of 
vector bosons, which are in the adjoint representation of the different groups: the 8 gluons for 
SU{3), the 3 W bosons of SU{2) and the vector of C/(l)y. 

Fermion mass terms emerge from the coupling of fermion bilinears with the scalar sector, which 
consists, in the SM, in a scalar field, doublet under SU{2) 



(A-2) 



The neutral component of (f) is singlet under SU{'S)c x [/"(l)em) and so it can take a vev (0°) = 
■u/v^ ~ 246 GeV, leaving the low-energy group invariant and in the meanwhile fermions aquire 
masses, thanks to the Yukawa couplings: 



V V V 

Ae -j=eLeR + Xd-j^ dL + Ad ul ur + h.c. 



(A-3) 



The SM does not contain a right-handed neutrino, hence m-u = 0. Furthermore, all interaction 
terms are renormalisable. This renormabilisity, together with the field content given above. 
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imply that besides the invariance under the local symmetry G321, 4 global U{l)s are conserved: 
the baryon number B, and the lepton flavour number Le^^^r, and consequently the total lepton 
number L. As a consequence of this, the neutrino remains massless to all orders of perturbation. 

A. 2 Sphalerons 

Baryon and Lepton numbers are accidental symmetries of the Standard Model. This means 
that they do not reflect a higher symmetry but are satisfied given the particle content and the 
renormalisable couplings of the model. B and L are conserved at tree level and to all order in 
perturbation theory. 

However, 'T Hooft showed [13] that non-perturbative effects, called instantons, can lead to the 
violation of B + L, while conserving the orthogonal B — L. This property is associated with the 
topological structure of any SU{N) gauge group. The ground state of the theory is not unique, 
but degenerate vacua exist which are topologically inequivalent. Going from one vacuum to 
another vacuum can be done by tunnelling through field configurations called instantons, but 
typically the probability will be highly suppressed. Another possibility arise from the existence of 
static but unstable field configurations that help the transition to occur: these are the sphalerons. 
In concrete terms, the baryonic and leptonic number currents, 

= ^^La^,iLi-eli''Jeli , 

i 

•^/^ = ^ *5l tlfiQi i - ul i^jul i - dl i^-fdl i , (A-4) 

i 

have a non-zero divergence from the ABJ triangle anomaly |128| : 

d'^J^ = d^jj: = ^ (g'^W^'W,, - g'^B^'^B-^,) , (A-5) 

where W^y and B^y are SU{2) and C/(l) field-strengths respectively, and Ng the number of 
fermions generations. Therefore, we see that -B — L is conserved, as dii{jj^ — jj^) = 0. The 
orthogonal combinaison, i3 -|- L is violated: 

dfi{Jll + Jl^) = 2Ngd^K'^ . (A-6) 
The interesting point is that JC^ is related to the topological structure of the vacuum, by: 

NcsiU) - Ncsito) = r dt I d^'xd^K^' = n , (A-7) 

J to J 

where n is an integer, and Nqs are Chern-Simons number. Therefore we see that if during a 
time At a transition between two vacua with distinct topological charges occurs, this will induce 
a baryon-|-lepton number violation, as 



A(B + L) = Ng^Ncs 



(A-8) 



A. 2 Sphalerons 
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At zero temperature however, the transition rate between two vacua is exponentially suppressed, 
with 

r ^ e-8-'/f' e~^'^ , (A-9) 

and effectively no transition occurs. The picture changes however at high temperature, T > Tewi 
as showed by Kuzmin, Rubakov and Schapovnikov[n] : instead of tunnelling through the barrier, 
the available thermal energy allows to step over it, thanks to the Higgs and gauge boson field 
configurations that lie on the top of the energy barrier separating two topologically inequivalent 
vacua. The situation is schematically depicted in fig jA-l[ 




Figure A-1: Energy of gauge field configuration versus Chern-Simons number 



The sphaleron transition amplitude is roughly e ^^p'"/^, with Egph the height of the energy 
barrier 

^ Sit V . . . „N 

Esph ^ — ■ (A-10) 

Hence, for T » Egph, the transition amplitude gets unsuppressed and rapidly occurs. 

In the Standard Model, we have the strong sphalerons, which are related to the SU{3) gauge 

structure, and whose density rate is given by |129| : 

Tqcd/V:^2504T^. (A-11) 

As y ~ T^, the QC-D-sphalerons enter in- equilibrium when Tqcd{T) < H{T), which happens 
at r < 10^'^ GeV, and then induce an effective interaction 

Oqcd = I^i=i,2,3iqLiqLiU%i d%i) , (A-12) 

which relate left and right-handed quarks. 

The SM also contains electroweak instantons, which are related to the SU{2) gauge structure, 
with a density rate |130) : 

Tew/V c^25a^wT^- (A-13) 
If in-equilibrium, these sphalerons induce an effective interaction 

OeW ='n^i=l,2,3{qLiqLiqLdLi) , (A-14) 
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which consist in a B ^ L exchange. The SU{2) sphalerons are faster than the Hubble expansion 
rate for T < 10^^ GcV. Hence, for Tew < 10^"^ GcV, any B + L number will be driven to zero 
by the fast B + L violating sphalerons. On the other hand, B~L asymmetries are not affected by 
sphalerons transition. The early GUT baryogenesis model were based on SU{5), which conserves 
B — L, hence the possible B number creation was wiped out as soon as electroweak sphalerons 
enter in-equilibrium and this models cannot work. On the contrary of leptogenesis models, were 
the asymmetry is produced in the L direction, which is then reprocessed into the —B direction. 



Appendix B 
Thermodynamics 



The evolution of the Universe is described by the Einstein equation which relate the geometry 
of the Universe to its content: 

- lg^^un = ^-kGnT^v + A^M^ . (B-1) 

In this equation TZ is the Ricci scalar, TZ^y the Ricci tensor, T^j^ the stress-energy tensor, the 
space-time metric and A is the cosmological constant. Gjq is the Newton constant, 

Mpl ~ 1.221 X lO^'^GeV . (B-3) 

Assuming that the Universe content is a perfect fluid, we can write the stress-energy tensor as: 

T^v = -pQfM^ + {p + p)u^Uv , (B-4) 

where p is the pressure and p the energy density of the perfect fluid. The velocity vector of the 
fluid, is given in the rest frame of the plasma b)' n = (1, 0, 0, 0). 

Assuming further homogeneity and isotropy of the Universe, it can be described through a 
Robertson- Walker metric, 

ds^ = dt^ - a{tfdf^ , (B-5) 
so that the (0-0) component of the Einstein equation can be written: 

a ) 3 3 

where A: = 0, -|-1, —1 is the curvature of the Universe (resp. flat, open and closed geometry), and 
a(t) is the scale factor of the expanding Universe. The Hubbler expansion rate R is defined as: 
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and tells us how fast the Universe is expanding. 

The species which compose the Universe are liable to interaction. In an expanding Universe, 
if this rate is faster than the Hubble expansion rate, then theses particles will be maintained 
causally connected with the thermal bath: the species are in- equilibrium and in that case they 
all have the same temperature. However, Uinverse is cooling. Then it is likely that at some 
temperature Tj, interactions will not be able to keep species in equilibrium as they freeze-out: 
r(r) < H{T) for T < Tf. 

When out-of-equilibrium, the species decouple from the thermal bath and evolves independently 
of it. The criterion T{Tf) ~ H{Tf) gives a rough picture of the decoupling, but a more accurate 
description require to study microscopic evolution of the particle number densities, by solving 
the Boltzmann equations. This appendix briefly introduce equilibrium and out-of-equilibrium 
thermodynamics, which is needed when studying leptogenesis. We refer the reader to refs.[I].[M] 
for more detailled presentations of section 1 and 3. 



B.l Equilibrium thermodynamics 

Species in- equilibrium can be described by ideal Fermi-Dirac (FD) or Bose-Einstein (BE) fluid, 
whose distribution function is: 

/r(p, ^„ T) = , (B-8) 

exp((Si - /Xij/T) ± 1 



where gi denotes the number of degree of freedom of the species i, Ei = y -|- mf and /Xj is 
the chemical potential of the species. The -|- (— ) sign refers to FD (BE) statistic. The classical 
approximation of Maxwell-Boltzmann (MB) statistic, which will extensively use in the network 
of Boltzmann equations, is: 

/f" = Qi exp{{Ei - fi,)/T) . (B-9) 

Given these distribution functions, one defines the equilibrium number density rij and the equi- 
librium energy density pi are given by: 



pT(T) = I ^,friP,f^uT)E{p) = ^l dE—i- (B-11) 



Defining Zi = rrii/T^ the number density for a massive particle is given by: 

27r2 



r^\p,Tr^ = '^e~'^z^K,(z,), (B-12) 



for Maxwell-Boltzmann statistic, and: 

<"(/., r)f g = 1^ E l+i ^'^2 ((fc + 1)^0 , (B-13) 

fc=0 



B.l Equilibrium thermodynamics 
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for fermions and bosons. Here K„ {x) is the modified Bessel function of the second kind. 
The asymptotic expansion of Bessel functions is well known and one can therefore deduce asymp- 
totic behaviour of the number densities. 

• For high temperatures T S> mj, neglecting chemical potential: 



7r2 V 4 ' 



where C(3) ~ 1.202. 

The low-temperature expansion reads: 



^eq^j^^MB ^ ^if^j^^BE ^ nl^Tf ~ gi ^'^ e-"^^/^ . (B-17) 

The usual assumption is to disregard quantum statistic and handle fermions and bosons in terms 
of MB particle, with either 



<=%j4i^2(^x), (B-18) 



for a massive particle, or 



= (B-19) 



for a massless one. Similarly, one obtains the asymptotic behaviour of the energy density: 
• for the high temperature expansion T 3> mj one has: 

^2 



p^^(r)^^ '^giT\ (B-20) 

Pfm^^ - y^9iT'. (B-21) 

For the low temperature expansion T <^ rrii, 

pf{T) = rrii n1\T) oc e""^*/^ . (B-22) 



Therefore, at a given temperature T, the energy density of non-rclativistic species will be expo- 
nentially suppressed, thus being safely neglected in the total energy density of the Universe 



^2 



Ptot = ^g^T" , (B-23) 



142 



Thermodynamics 



where counts the effective number of degrees of freedom: 

i=BE ^ ^ i=FD ^ ^ 

Here, Tj is the temperature of the species i: Ti oc T with equality when the species i are 
in- equilibrium. For the typical temperatures we consider, T » 100 GeV, g^, = 106.75 in the 
Standard Model, while g^ = 228.75 in the MSSM. 



Therefore, according to eq.(B-6), assuming flatness and negligible cosmological constant, that is 
valid in the very early Universe, one deduces 



'2 



In the radiation dominated epoch a ~ t^/^ and therefore one relates time and temperature by: 



-2 



Neglecting the chemical potential we have for the entropy density: 



+ (B.27) 



T 45 
where 

i=BE ^ ' i=FD ^ ' 

If all relativistic species are in-equilibrium, then Ti = T and q^: = g^,. Furthermore, since 
T oc a~i and s oc a~^, as long as the number of particles remains constant, the comobile number 
density Yi = rij/s is constant. On the other hand, when degrees of freedom freeze and become 
non-relativistic, the conservation of entropy implies that qf = q^ and thus: 

r> = r< (jPj ' . (B-29) 

That is, when dof decouple from the thermal bath, their entropy is transfered to species that are 
in thermal equilibrium, and their temperature therefore increase: entropy creation (re)heats the 
Universe. 

Given the above definitions, we further have that 

s = 45*^^(3^ ^7 - l-8g*ra^, (B-30) 
so that the actual entropy and photon number density are related according to 

So ^ 1.8 (7*0 — 7.04 n^o ■ (B-31) 



B.2 Relating B and L numbers 
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B.2 Relating B and L numbers 

Let us see how to relate the baryon number of the Universe to the lepton number that is produced 
during leptogenesis [HS]. Since that particles and antiparticles have an opposite chemical potential 
IJ-x = —fJ'x, the net number density of a particle x can be written, for small Hx/T 

9xT^ / fJ-x/T for a fermion /tj oo^ 

o I Zjix/1 for a boson 

Therefore, since B is stored both in SU{2)l doublet and singlet components, the baryon number 
is 



6 

i 

where the sum is made over the different asymmetries that are populated. 
Similarly, the lepton number is 

L = Y^Li = (^/x^^ + fi^y + /Xg.^) . (B-34) 

i i 

The different species undergo reactions, which if in- equilibrium enforce algebraic relations among 
the chemical potentials involved, since foiSii + j-^k + l interaction, the equilibrium condition 
implies 

Mi + = /"fe + W • (B-35) 

Then, according to the different chemical equilibriums that hold at a given temperature, B and 
Li can be related. Actually in the SM, B and L are only related through the effective interaction 
induced by the electroweak sphalerons. Since this interaction is in- equilibrium for T < lO^^ Q^y^ 
no baryon asymmetry can be generated via leptogenesis for T > 10^^ GeV. Furthermore, since 
these interactions conserve B — L but violate B + L, which they set to zero, one could naively 
think that 

B + L B-L^B-L 

2 2 2 ^ ' 

when electroweak sphalerons are in-equilibrium. Actually, since the latter only involve left- 
handed fields, while the lepton and baryon number are stored both in SU{2) doublets and 
singlets, the situation is more involved, and relating B to the Lj require to determine which 
interactions are in-equilibrium when the lepton asymmetry is created. Since leptogenesis occurs 
at T ~ Ml, the different temperature regime then correspond to range of Mi. 

We can list the different interactions and the temperature at which they become in- equilibrium [53] . |76| . 

• Gauge interactions are always in-equilibrium. Consequently, gluons and B^ have a null 
chemical potential, and furthermore the different quark colours have a same fi. Interactions 
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with the enforce 

fJ-w = ^^ij)- + fJ'4>o , 

fiw = -fJ-UL + I^BL , (B-37) 
where nw = fJ'W- — ~l^w+^ denotes the Higgs doublet. 

Moreover, since SU{2) is gauged, the sum of third weak-isospin component should be zero: 
^ Yj {j'ui - /^dl) + ^ {j'l^l ~ ~ ^ (^'t'° + - = , (B-38) 

i i 

where m is the number of Higgs doublet. This relation implies that /j^w = and so 
all SU{2) doublet components have an equal chemical potential fiuL = Md^ = 

• The requirement that the total hypercharge is null enforces: 

J2 {iJ'qi + 2 - l^d},) - Yl (^4 + l^^l) + 2 /^^ = . (B-39) 

i i 

• For T < 10^3 QgY^ sphalcrons are in- equilibrium: 

i 

For T < 10^^ GeV, the electroweak sphalerons are in- equilibrium and so B and L doublets 
are related according to: 

^(3/Xgj^+/x^.J =0 . (B-41) 

i 

Finally, the interactions involving charged fermion Yukawas, when in- equilibrium, yield the fol- 
lowing relations: 

for up-type quarks, down-type quarks and charged leptons, respectively. 

Wc arc interested in relating the = B/3 — asymmetries to the leptonic doublet ones iia, 
and this relation depends on the temperature regime. 

• If Ml > 10^^ GeV, none of the charged lepton Yukawas are in-equilibrium. Therefore, 
the SU(2)l singlets do not store any as^^mmetry since only the left-handed part interact 
with the gauge bosons. Consequently the lepton number and the leptonic doublet number 
are equal. Since no baryon number is generated, as electroweak sphalerons are out-of- 
equilibrium, we have Yao = = —Yl- 



B.3 Out-of-equilibrium thermodynamics: Boltzmann equations 



145 



If 10^ < Ml < 10^2 QgY^ tijgn the tau- Yukawas are in- equilibrium, as well as electroweak 
sphalerons. A non-zero asymmetry can thus be stored in tr. Since t, b and c quarks are 
in- equilibrium, the conversion matrix reads: 



( ] 


T3 


1 

9 




~ 6s 




\ y^r ) 




I - 



i ^ 








I- 











(B-43) 



However, the e and ^ flavours are indistinguishable, therefore one should rather consider 
the asymmetry produced m. e + ^ direction, that is, one should sum the chemical potentials 
/Xe^ + fif,^ , SO that 



Since Yi 




(B-44) 



■YAa. 



(B-45) 



If Ml < 10^ GeV, the muon Yukawas are in- equilibrium, as well as the s-ones. Thus, the 
three lepton flavours are distinguishable, and one has 



( \ 


rp3 ( 








V 



that is, 



Yi^ 



1 
179 



V 



9 






1-2 


11 


4 


9 


15 


9 


4 


4 


4 


9 


15 


9 


151 


20 


20 


25 


344 


14 


2 


3 


3 


25 


14 


344 


2 


3 


3 




(B-46) 



■Yac 



(B-47) 



Finally, the baryon asymmetry is related to the B/3 — La ones by the relation 

32 + 4m 



Yb 



8 + 13m 

OL 

which gives Yb = 12/37 Yaq in the SM with (m = 1) one Higgs doublet. 



(B-48) 



B.3 Out-of-equilibrium thermodynamics: Boltzmann equations 

Let us see how to characterise species which undergo non-elastic collisions. The evolution of the 
distribution function is given by the Boltzmann equations (BE) 
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where C is the collision term and L is the Liouville operator, 

which is conveniently express in a homogeneous and isotopric space-time as : 

L[f] = E^-Hit)p^. (B-49) 

Integrating over momentum, one recover the usual BE: 

rii + 3H{t)ni = d . (B-50) 

All the work will then to evaluate the collision term. For a generic 2-to-2 scattering i+j ^ k + l, 
it can be written as: 



d = {2Trf J S{pi + pj -Pk- Pi) 



X \M\^x{{l±m±fj)hh-hfj{l±fk){l±fi)) , (B-51) 

where + (-) stands for FD (BE) statistics, and M. is the invariant amplitude. We can simplify 
this expression by using MB statistics and neglecting Pauli blocking and stimulating emission 
factors (i.e., quantum statistic), so that the collision term reads: 



Ci = (27r)'^ J d{pi + Pj -pk- Pi) 



2{2TT)^Ei"'2{2TT)^Ei 
X \M\''xifkfi-fifj) . (B-52) 

Since that for what concern us, the elastic scatterings, which conserve particle densities but 
not their distribution, are much faster than inelastic scatterings, we can approximate fx = 
f^nx/rfx' so that the collision term can be further simplified: 

'~ r,^ r,^1'^J e<l eq'kl ' y° 

where we define the reaction density 7^]: 

7g = (2.r/i(« + »-„-„)A...._«^ X lA^Px/r/f. (B-54) 



Given this, the BE reads: 



rii 



+ 3 Hit) n. = E I S ^^7^' - ^ ^7^] 1 . (B-55) 



B.4 Boltzmann equations in leptogenesis 
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B.4 Boltzmann equations in leptogenesis 

It is customary in leptogenesis to parametrise number density evolution in function of the variable 

z = Mi/T , (B-56) 

instead of the time t, where M. The relation between those two is 

dz Mn^ dT 

-dt=-^Tt='''^'^- ^^-^^^ 

It is also customary to use the comoving number density Yx that are defined by 

^ — , (B-58) 
s 

where s is the entropy density. Such transformation is done because both Ux and s scale as T^, 
and therefore Y is constant during the cooling of the Universe, up to non-elastic processes. 
With these two redefinitions, the BE eq.(B-50| translates to 

X,j \ k I i j / 

This is the general form of the Boltzmann equations for comoving number densities. It is written 
here for a two-to-two scattering, but the collision term should obviously include all relevant 
processes of the model under consideration. In this manuscript, we only focus on decays X ij 
and on 2-to-2 scatterings. The density rate for a two body decay is given by: 



ifj (z.) = nl^{zx) X r(X - ij) , (B-61) 



where Ki^2 are the modified Bessel function of the second kind, and T{X ij) is the decay rate 
evaluated in the rest-frame of the decaying particle. For a two body scattering, one has: 

^ki'^(') = lUr-A— / dx^K^{^z^)a{x) , (B-62) 
647r4 Zm Ji 

where Smin = max ((mx + rnY)"^, {rrik + Tni)"^), Zm = \/Smin/T and d{x) is the reduced cross- 
section, the cross section summed over initial phase space: 

^ 2A(s,m|,m2,) 2 (s - (mx + my)^) (s - (mx - my)^) 

a[s) = —a[s) = — ^ — '^{s) ■ (B-63j 

s s 

Decays/inverse decays and top-scatterings 

Let us discuss the different terms that we consider in this thesis. 
At first order, we included decays and inverse, that are given by: 

Ki{z) 

N 



'^^D{z) = =KiZ 

sH{Mi) K2{z 



2 y/«(z) "^^^^ " 4 



i^iWUz) = D{z) = ^ z' K^iz) . (B-64) 
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where ki (denoted Ki in the thesis), is the wash-out parameter, 



Given the high temperature expansion 



2 . 2 

(B-66) 



z z^ 



we have that D{z ^ 1) ~ 2Wid{z <^ \)z^ jl <^ 1: decays and inverse decays are highly 
suppressed. In the opposite low-temperature regime, since 

we obtains 

Z)(z>l) ~ z-^, 

Hence, while decays increase with z, the inverse decays get Boltzmann suppressed at low tem- 
perature. The latter reach their maximum for z ~ 2.4, at Wid{z) ~ 0.29. 

We also include AL = 1 neutrino-top scatterings, which are mediated by a Higgs field in the s— 
and t— channels, as seen in chapter 3. 

The reduced cross-section for these scatterings, involving a right-handed neutrino Ni, is given 
by [S7]: 

a,(x) = ^^,^^, (B-69) 
2 

where x = s/Mf is the rescaled energy in the centre-of-mass frame, ai = Mf /Mf and at = h^/47r. 
The t— channel contribution reads |57| : 

^ , . o X - ai f X -2ai + 2ah ai-2ah , f x - at + ah\\ ^r, r^n\ 

as{x) = 3 Of Ki ■ \ In , (B-70) 

X \ X - ai + ah X - ah \ an J J 

where ah correspond to the rescaled Higgs mass, which is used to regularise an infrared divergence, 
ah = (a^/Mi)2 ~ 10-16. 

At high temperatures, the reaction density scales as T^, and one has [5^ : 

z ^s,t z<i 9 mj 



and D(z) + S{z) ~ S{z) ~ 0.1, z <Cl. 

In the opposite regime of low temperatures, D{z) + S{z) ~ D{z) ~ z. 
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